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We survey the basic principles of atmospheric dynamics relevant to explaining existing 
and future observations of exoplanets, both gas giant and terrestrial. Given the paucity of data 
on exoplanet atmospheres, our approach is to emphasize fundamental principles and insights 
gained from Solar-System studies that are likely to be generalizable to exoplanets. We begin 
by presenting the hierarchy of basic equations used in atmospheric dynamics, including the 
Navier-Stokes, primitive, shallow-water, and two-dimensional nondivergent models. We then 
survey key concepts in atmospheric dynamics, including the importance of planetary rotation, 
the concept of balance, and simple scaling arguments to show how turbulent interactions 
generally produce large-scale east-west banding on rotating planets. We next turn to issues 
specific to giant planets, including their expected interior and atmospheric thermal structures, 
the implications for their wind patterns, and mechanisms to pump their east-west jets. Hot 
Jupiter atmospheric dynamics are given particular attention, as these close-in planets have been 
the subject of most of the concrete developments in the study of exoplanetary atmospheres. We 
then turn to the basic elements of circulation on terrestrial planets as inferred from Solar-System 
studies, including Hadley cells, jet streams, processes that govern the large-scale horizontal 
temperature contrasts, and climate, and we discuss how these insights may apply to terrestrial 
exoplanets. Although exoplanets surely possess a greater diversity of circulation regimes than 
seen on the planets in our Solar System, our guiding philosophy is that the multi-decade study 
of Solar-System planets reviewed here provides a foundation upon which our understanding of 
more exotic exoplanetary meteorology must build. 



1. EVTRODUCTION 

The study of atmospheric circulation and climate began 
hundreds of years ago with attempts to understand the pro- 
cesses that determine the distribution of surface winds on 
the Earth (e.g., Hadley 1735). As theories of Earth's gen- 
eral circulation became more sophisticated (e.g., Lorenz 
1967), the characterization of Mars, Venus, Jupiter, and 
other Solar-System planets by spacecraft starting in the 
1960s demonstrated that the climate and circulation of other 
atmospheres differ, sometimes radically, from that of Earth. 
Exoplanets, occupying a far greater range of physical and 
orbital characteristics than planets in our Solar System, like- 
wise plausibly span an even greater diversity of circulation 
and climate regimes. This diversity provides a motivation 
for extending the theory of atmospheric circulation beyond 
our terrestrial experience. Despite continuing questions, 
our understanding of the circulation of the modem Earth 
atmosphere is now well developed (see, e.g.. Held 2000; 
Schneider 2006; Vallis 2006), but attempts to unravel the 
atmospheric dynamics of Venus, Jupiter, and other Solar- 
System planets remain ongoing, and the study of atmo- 



spheric circulation of exoplanets is in its infancy. 

For exoplanets, driving questions fall into several over- 
lapping categories. First, we wish to understand and explain 
new observations constraining atmospheric structure, such 
as light curves, photometry, and spectra obtained with the 
Spitzer, Hubble, or James Webb Space Telescopes (JWST), 
thus helping to characterize specific exoplanets as remote 
worlds. Second, we wish to extend the theory of atmo- 
spheric circulation to the wide range of planetary param- 
eters encompassed by exoplanets. Existing theory was pri- 
marily developed for conditions relevant to Earth, and our 
understanding of how atmospheric circulation depends on 
atmospheric mass, composition, stellar flux, planetary rota- 
tion rate, orbital eccentricity, and other parameters remains 
rudimentary. Significant progress is possible with theoret- 
ical, numerical, and laboratory investigations that span a 
wider range of planetary parameters. Third, we wish to un- 
derstand the conditions under which planets are habitable, 
and answering this question requires addressing the inter- 
twined issues of atmospheric circulation and climate. 

What drives atmospheric circulation? Horizontal tem- 
perature contrasts imply the existence of horizontal pres- 
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Fig. 1 . — Atmospheric circulation results from a coupled interaction between radiation and hydrodynamics: horizontal temperature and 
pressure contrasts generate winds, which drive the atmosphere away from local radiative equilibrium. This in turn allows the spatially 
variable thermodynamic (radiative) heating and cooling that maintains the horizontal temperature and pressure contrasts. 



sure contrasts, which drive winds. The winds in turn push 
the atmosphere away from radiative equilibrium by trans- 
porting heat from hot regions to cold regions (e.g., from the 
equator to the poles on Earth). This deviation from radia- 
tive equilibrium allows net radiative heating and cooling to 
occur, thus helping to maintain the horizontal temperature 
and pressure contrasts that drive the winds (see Fig. 1). Spa- 
tial contrasts in thermodynamic heating/cooling thus funda- 
mentally drive the circulation, yet it is the existence of the 
circulation that allows these heating/cooling patterns to ex- 
ist. (In the absence of a circulation, the atmosphere would 
relax into a radiative-equilibrium state with a net heating 
rate of zero.) The atmospheric circulation is thus a coupled 
radiation-hydrodynamics problem. On the Earth, for exam- 
ple (see Fig. 2), the equator and poles are not in radiative 
equilibrium. The equator is subject to net heating, the poles 
to net cooling, and it is the mean latitudinal heat transport 
that is both responsible for and driven by these net imbal- 
ances. 

The mean climate (e.g., the global-mean surface temper- 
ature of a planet) depends foremost on the absorbed stel- 
lar flux and the atmosphere's need to reradiate that energy 
to space. Yet even the global-mean climate is strongly af- 
fected by the atmospheric mass, composition, and circula- 
tion. On a terrestrial planet, for example, the circulation 
helps to control the distribution of clouds and surface ice, 
which in turn determine the planetary albedo and the mean 
surface temperature. In some cases, a planetary climate can 
have multiple equilibria (e.g., a warm, ice-free state or a 
cold, ice-covered "snowball Earth" state), and in such cases 
the circulation plays an important role in determining the 
relative stability of these equilibria. 

Understanding the atmosphere/climate system is chal- 
lenging because of its nonlinearity, which involves multiple 
positive and negative feedbacks between radiation, clouds, 
dynamics, surface processes, planetary interior, and life (if 
any). The inherent nonlinearity of fluid motion further im- 
plies that even atmospheric-circulation models neglecting 
the radiative, cloud, and surface/interior components can 



exhibit a large variety of behaviors. 

From the perspective of studying the atmospheric circu- 
lation, transiting exoplanets are particularly intriguing be- 
cause they allow constraints on key planetary attributes that 
are a prerequisite to characterizing an atmosphere's circu- 
lation regime. When combined with Doppler velocity data, 
transit observations permit a direct measurement of the ex- 
oplanet's radius, mass and thus surface gravity'. With the 
additional expectation that close-in exoplanets are tidally 
locked if on a circular orbit, or pseudo-synchronized^ if on 
an eccentric orbit, the planetary rotation rate is thus indi- 
rectly known as well. Knowledge of the radius, surface 
gravity, rotation rate and external irradiation conditions for 
several exoplanets, together with the availability of direct 
observational constraints on their emission, absorption and 
reflection properties, opens the way for the development of 
comparative atmospheric science beyond the reach of our 
own Solar System. 

The need to interpret these astronomical data reliably, 
by accounting for the effects of atmospheric circulation and 
understanding its consequences for the resulting planetary 
emission, absorption and reflection properties, is the cen- 
tral theme of this chapter. Tidally locked close-in exoplan- 
ets, for example, are subject to an unusual situation of per- 
manent day/night radiative forcing, which does not exist in 
our Solar System^. To address the new regimes of forcings 
and responses of these exoplanetary atmospheres, a discus- 
sion of fundamental principles of atmospheric fluid dynam- 
ics and how they are implemented in multi-dimensional, 
coupled radiation-hydrodynamics numerical models of the 
GCM (General Circulation Model) type is required. 

Contemplating the wide diversity of exoplanets raises 



^Combining Doppler velocity and transit measurements lifts the mass- 
inclination degeneracy. 

■^Pseudo-synchronization refers to a state of tidal synchronization achieved 
only at periastron passage (=closest approach), as expected from the strong 
dependence of tides with orbital separation. 

Venus may provide a partial analogy, which has not been fully exploited 
yet. 
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Fig. 2. — Earth's energy balance. The Earth absorbs more sunlight at the equator than the poles (blue curve, denoted "shortwave"). 
The Earth also radiates more infrared energy to space at the equator than the poles (red curve, denoted "longwave"). However, because 
the atmospheric/oceanic circulation act to mute the latitudinal temperature contrasts (relative to radiative equilibrium), the longwave 
radiation exhibits less latitudinal variation than the shortwave absorption. Thus, the circulation leads to net heating at the equator and net 
cooling at the poles, which in turn drives the circulation. Data are an annual-average for 1987 obtained from the NASA Earth Radiation 
Budget Experiment (ERBE) project. Copyright M. Pidwirny, www.physicalgeography.net, used with permission. 



a number of fundamental questions. What determines the 
mean wind speeds, direction, and 3D flow geometry in at- 
mospheres? What controls the equator-to-pole and day- 
night temperature differences? What controls the frequen- 
cies and spatial scales of temporal variability? What role 
does the circulation play in controlling the mean climate 
(e.g., global-mean surface temperature, composition) of an 
atmosphere? How do these answers depend on parame- 
ters such as the planetary rotation rate, gravity, atmospheric 
mass and composition, and stellar flux? And, finally, what 
are the implications for observations and habitability of ex- 
oplanets? 

At present, only partial answers to these questions exist 
(see reviews by Showman et al. 2008b; Cho 2008). With up- 
coming observations of exoplanets, constraints from Solar- 
System atmospheres, and careful theoretical work, signif- 
icant progress is possible over the next decade. While a 
rich variety of atmospheric flow behaviors is realized in the 
Solar System alone — and an even wider diversity is pos- 
sible on exoplanets — the fundamental physical principles 
obeyed by all planetary atmospheres are nonetheless uni- 
versal. With this unifying notion in mind, this chapter pro- 
vides a basic description of atmospheric circulation princi- 
ples developed on the basis of extensive Solar-System stud- 
ies and discusses the prospects for using these principles to 
better understand physical conditions in the atmospheres of 
remote worlds. 

The plan of this chapter is as follows. In Section 2, we 
introduce several of the equation sets that are used to in- 
vestigate atmospheric circulation at varying levels of com- 
plexity. This is followed (Section 3) by a tutorial on ba- 
sic ideas in atmospheric dynamics, including atmospheric 
energetics, timescale arguments, force balances relevant to 
the large-scale circulation, the important role of rotation in 
generating east-west banding, and the role of waves and 



eddies in shaping the circulation. In Section 4, we survey 
the atmospheric dynamics of giant planets, beginning with 
generic arguments to constrain the thermal and dynamical 
structure and proceeding to specific models for understand- 
ing the circulation of our "local" giant planets (Jupiter, Sat- 
urn, Uranus, Neptune) as well as hot Jupiters and hot Nep- 
tunes.'* In Section 5, we turn to the climate and circulation 
of terrestrial exoplanets. Observational constraints in this 
area do not yet exist, and so our goal is simply to summa- 
rize basic concepts that we expect to become relevant as 
this field expands over the next decade. This includes a 
description of climate feedbacks (Section 5.1), global cir- 
culation regimes (Section 5.2), Hadley-cell dynamics (Sec- 
tion 5.3), the dynamics of the so-called midlatitude "baro- 
clinic" zones where baroclinic instabilities dominate (Sec- 
tion 5.4), the slowly rotating regime relevant to Venus and 
Titan (Section 5.5), and finally a survey of how the circu- 
lation responds to the unusual forcing associated with syn- 
chronous rotation, extreme obliquities, or extreme orbital 
eccentricities (Section 5.6). The latter topics, while perhaps 
the most relevant, are the least understood theoretically. In 
Section 6 we summarize recent highlights, both observa- 
tional and theoretical, and in Section 7 we finish with a sur- 
vey of future prospects. 

2. EQUATIONS GOVERNING ATMOSPHERIC CIR- 
CULATION 

A wide range of dynamical models has been developed 
to explore atmospheric circulation and climate. Such mod- 



''The terms hot Jupiter and hot Neptune refer to giant exoplanets with 
masses comparable to those of Jupiter and Neptune, respectively, with or- 
bital semi-major axes less than 0.1 AU, leading to high temperatures. 
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els are used to explain observations, understand mecha- 
nisms that govern the circulation/climate system, determine 
the sensitivity of a planet's circulation/climate to changes 
in parameters, test hypotheses about how the system works, 
and make predictions. 

Developing a good understanding of atmospheric circu- 
lation requires the use of a hierarchy of atmospheric fluid 
dynamics models. Complex models that properly repre- 
sent the full range of physical processes may be required 
for detailed predictions or comparisons with observations, 
but their very complexity can obscure the specific physi- 
cal mechanisms causing a given phenomenon. In contrast, 
simpler models contain less physics, but they are easier to 
diagnose and can often lead to a better understanding of 
cause-and-effect in an idealized setting. Whether a given 
model contains sufficient physics to explain a given phe- 
nomenon is a question that can only be answered by ex- 
ploring a hierarchy of models with a range of complexity. 
Exploring a hierarchy of models is therefore invaluable be- 
cause it allows one to determine the minimal set of physi- 
cal ingredients that are needed to generate a specific atmo- 
spheric behavior — insight that typically cannot be obtained 
from one type of model alone. 

The equations governing atmospheric behavior derive 
from conservation of momentum, mass, and energy for a 
fluid, which we here assume to be an electrically neutral 
continuum. For three-dimensional models, where momen- 
tum is a three-dimensional vector, this implies five govern- 
ing equations, which are generally represented as five cou- 
pled partial differential equations for the three-dimensional 
velocity, density, and internal energy per mass (with other 
thermodynamic state variables determined from density and 
internal energy by the equation of state). The Navier-Stokes 
equations, described in § 2.1, constitute the canonical exam- 
ple and provide a complete representation of a continuum, 
electrically neutral, viscous fluid in three dimensions. 

So-called reduced models simplify the dynamics in one 
or more ways, for example by reducing specified equations 
to their leading-order balances. For example, because most 
atmospheres have large aspect ratios (with characteristic 
horizontal length scales for the global circulation typically 
10-100 times the characteristic vertical scales), the vertical 
momentum balance is typically close to a local hydrostatic 
balance, with the local weight of fluid parcels balancing 
the local vertical pressure gradient [see, e.g., Holton (2004, 
pp. 41-42) or VaUis (2006, pp. 80-84) for a derivation]. The 
primitive equations, described in § 2.2, formalize this fact 
by replacing the full vertical momentum equation with lo- 
cal vertical hydrostatic balance. Although the system is still 
governed by five equations, this alteration simplifies the dy- 
namics by removing vertically propagating sound waves, 
which are unimportant for most meteorological phenom- 
ena. It also leads to mathematical simplification, making 
it easier to obtain analytic and numerical solutions. This is 
the equation set that forms the basis for most cutting-edge 
global-scale climate models used for studying atmospheres 
of Solar-System planets, although some global-scale high- 



resolution models now include non-hydrostatic effects. 

A further common reduction is to simplify the dynam- 
ics to a one-layer model representing (for example) the ver- 
tically averaged flow. The most important example is the 
shallow-water model, described in § 2.3, which govern the 
behavior of a thin layer of constant-density fluid of variable 
thickness. This implies three coupled equations for hor- 
izontal momentum and mass conservation (governing the 
evolution of the two horizontal velocity components and 
the layer thickness) as a function of longitude, latitude, and 
time. Although highly idealized, the shallow-water model 
has proven surprisingly successful at capturing a wide range 
of atmospheric phenomena and has become a time-honored 
process model in atmospheric dynamics (see, for example, 
Pedlosky 1987, Chapter 3). 

A further reduction results from assuming the fluid- 
layer thickness is constant in the shallow-water model. 
Given that density is also constant, the mass conservation 
equation then becomes a statement that horizontal conver- 
gence/divergence is zero. This constraint, which has the 
effect of removing gravity (buoyancy) waves from the sys- 
tem, allows the horizontal velocity components to be rep- 
resented using a streamfunction, leading finally to a single 
governing partial differential equation for the streamfunc- 
tion as a function of longitude, latitude, and time. § 2.4 
describes this two-dimensional, non-divergent model. The 
impressive reduction from five coupled equations in five 
dependent variables (as for the Navier-Stokes or primitive 
equations) to one equation in one variable leads to great 
mathematical simplification, enabling analytic solutions in 
cases when they are otherwise difficult to obtain. Moreover, 
the exlusion of buoyancy effects, gravity waves, and verti- 
cal structure leads to a conceptual simplification, allowing 
the exploration of (for example) vortex and jet formation in 
the most idealized possible setting. 

A comparison of results from the full range of models 
described here provides a path toward identifying the rel- 
ative roles of acoustic waves, vertical structure, buoyancy 
effects, and gravity waves in affecting any given meteoro- 
logical phenomenon of interest. We now present the equa- 
tions associated with each of these models. 

2.1. Navier-Stokes Equations 

Let u = u(x, t) be velocity at position x and time t, 
where x, u G M'^. If the frictional force per unit area of the 
fluid is hnearly proportional to shear in the fluid, then it is 
a Newtonian fluid (e.g., Batchelor 1967). Such fluids are 
described by the Navier-Stokes equations: 
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is the material derivative (i.e., the deriative following the 
motion of a fluid element). Here, p is density, p is pres- 
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sure, ff, represents various body forces per mass (e.g., grav- 
ity and Coriolis), /i is molecular dynamic viscosity, and 
e = i [(Vu) + (Vu)^] and I are the strain-rate and unit ten- 
sors, respectively. In Eq. (la) the quantity inside the braces 
is the viscous stress tensor. Here, as in Eq. (3) below, the 
average normal viscous stress (bulk viscosity) has been as- 
sumed to be zero. 

Eq. (la) is closed with the following equations for mass 
(per unit volume), internal energy (per unit mass), and state: 



De 
Dt 



Dp 
Dt 



(V-u) 



= -pV-u, 



(2) 
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e:(Vu)T--(V-u)2 



+ -V • {K,VT) + Q, 
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P = Pip,T), 



(3) 



(4) 



where e = e{T, s) is specific internal energy, s is specific 
entropy, Kt- is heat conduction coefficient, T is tempera- 
ture, and Q is thermodynamic heating rate per mass. In 
Eq. (3) ":" is scalar-product (i.e., component-wise multi- 
plication) operator for two tensors. Eqs. (1-4) constitute 6 
equations for 6 independent unknowns, {u,p, p, T}. Note 
that for a homogeneous thermodynamic system, which in- 
volves a single phase, only two state variables can vary in- 
dependently; hence, there are only two thermodynamic de- 
grees of freedom for such a system. 

The neutral atmosphere is well described by Eqs. (1-4) 
when the characteristic length scale L is much larger than 
the mean free path of the constituents that make up the at- 
mosphere. Hence, the equations are valid up to heights 
where ionization is not significant and the continuum hy- 
pothesis does not break down. Under normal conditions, 
the atmosphere behaves like an ideal gas. The parameters 
jj,, Kj, and other physical properties of the fluid depend on 
T, as well as p. When appreciable temperature differences 
exists in the flow field, these properties must be regarded as 
a function of position. For large-scale atmosphere appUca- 
tions, however, the terms involving p in Eqs. (la) and (3) 
are small and can be neglected in most cases. The typical 
boimdary conditions are u • n = at the lower boundary, 
where n is the normal to the boundary, and p,p as 
z ^ oo. For local, limited area models, periodic boundary 
conditions are often used. 

2.2. The Primitive Equations 

On the large scale (to be more precisely quantified be- 
low), the motion of an atmosphere is governed by the prim- 
itive equations. They read (e.g., Salby 1996): 
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(5d) 

(5e) 



Note here that p, rather than the geometric height z, is 
used as the vertical coordinate. This coordinate, which 
simplifies the gradient term in Eq. (5a), is common in at- 
mospheric studies; it renders z = z{x,p,t) a dependent 
variable, where now x e M^. In Eq. (5) v(x, = {u,v) 
is the (eastward^, northward) velocity in a frame rotating 
with CI, where CI is the planetary rotation vector as repre- 
sented in inertial space; ^ = gz is the geopotential, where 
g is the gravitational acceleration (assumed to be constant 
and to include the centrifugal acceleration contribution; see 
Holton 2004, pp. 13-14) and z is the height above a fidu- 
cial geopotential surface; k is the local upward unit vector; 
f = 2Q sin (j) is the Coriolis parameter, the locally vertical 
component of the planetary vorticity vector 2fl; Vp is the 
horizontal gradient on a p-surface; zu = Dp/Dt is the ver- 
tical velocity; and T> represent the momentum sources 
and sinks, respectively; 6 = T{p„t/p)'^ is the potential tem- 
perature^, where p„r is a reference pressure and k = R/cp 
with R the specific gas constant and Cp the specific heat at 
constant pressure; and ^net is tiie net diabatic heating rate 
(heating minus cooling). Note that 4ici can include not only 
radiative heating/cooling but latent heating and, at low pres- 
sures where the thermal conductivity becomes large, con- 
ductive heating. The Newtonian cooling scheme, which re- 
laxes temperature toward a prescribed radiative-equilibrium 
temperature over a specified radiative time constant, is one 
simple parameterization of (jnct - 

The fundamental presumption in the use of Eqs. (5) is 
that small scale processes are parameterizable within the 
framework of large-scale dynamics. Here by "large" scales, 
it is meant typically L > a/ 10, where a is the planetary 
radius. By "small" scales, it is meant those scales that 
are not resolvable numerically by global models — typically 
< a/10. Regions of the atmosphere where small scale 
processes are important are often highly concentrated (e.g., 
fronts and convective updrafts). Their characteristic scales 
are <C a/ 10. Therefore, it is possible that the Eq. (5) set — as 
with all the other equation sets discussed in this chapter — 
leaves out some processes important for large-scale dynam- 
ics. 

To arrive at Eq. (5), one begins with Eqs. (1^) in spher- 
ical geometry (e.g., Batchelor 1967). Two approximations 
are then made. These are the "shallow atmosphere" and 



(5b) 



^Cardinal directions are defined here consistent with everyday usage; east 
is defined to be the prograde direction, that is, the direction in which the 
planet rotates. North is the direction along which SI • k becomes more 
positive. 

''The potential temperature 9 is related to the entropy s by ds = CpdlnO. 
When Cp is constant, this yields 9 = T{p^g{/p)'^. 
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the "traditional" approximations (e.g., Salby (1996)). The 
first assumes zja <C 1. The second is formally valid 
in the hmit of strong stratification, when the Prandtl ra- 
tio (TVV^^) > 1- Here, N = N{-x,z,t) is the Brunt- 
Vaisala (buoyancy) frequency, the oscillation frequency for 
an air parcel that is displaced vertically under adiabatic con- 
ditions: 



N = 



d{\ii0) 
dz 



1/2 



(6) 



These approximations allow the CorioUs terms involving 
vertical velocity to be dropped from Eq. (la) and vertical 
accelerations to be assumed small. The latter is exphcitly 
embodied in Eq. (5b), the hydrostatic balance, which we 
discuss further below. 

Hydrostatic balance renders the primitive equations vahd 
only when N^/to'^ » 1, where 27r/a; is the timescale of the 
motion under consideration. This condition, which is dis- 
tinct from the Prandtl ratio condition, restricts the vertical 
length scale of motions to be small compared to the hori- 
zontal length scale. Therefore, the hydrostatic balance ap- 
proximation breaks down in weakly stratified regions (here 
we refer to the dynamically evolving hydrostatic balance as- 
sociated with circulation-induced perturbations in pressure 
and density; the mean background density and pressure — 
i.e., those that would exist in absence of dynamics — will 
remain hydrostatically balanced even when the circulation- 
induced perturbations are not). The hydrostatic assumption 
filters vertically propagating sound waves from the equa- 
tions. 

According to Eq. (5d), when ^net = 0, individual values 
of are retained by fluid elements as they move with the 
flow. In this case, Eq. (5) also admit a dynamicaUy impor- 
tant conserved quantity, the potential vorticity: 



(C + /)k 



(7a) 



where ^ = k • V x v is the relative vorticity. This quantity 
provides the crucial connection between the primitive equa- 
tions and the physically simpler models that follow. For ex- 
ample, undulations of potential vorticity are often a direct 
manifestation of Rossby waves, which are represented in 
all the models presented in this section. The conservation 
of the potential vorticity following the flow. 



-PgpE 

Dt 



= 0, 



(7b) 



and the redistribution of implied by it, is one of the most 
important properties in atmospheric dynamics. 

2.3. ShaUow- Water Model 

For many applications, Eq. (5) is too complex and broad 
in scope. In the absence of observational information to 
properly constrain the model parameters, reduction of the 
equations is beneficial. A commonly used approach is to 
collapse the 3D primitive equations to a two-dimensional 



(2D), one-layer model. Such reduction aUows investigation 
of horizontal vortex and jet interactions in an idealized set- 
ting. 

Among the most widely used one-layer models is the 
shaUow-water model. Consider a thin layer of homoge- 
neous (i.e., constant-density) fluid, bounded above by a free 
surface and below by an impermeable boundary, so that 
its thickness is h{x,t). The dynamics of such a layer is 
governed by the following equations (e.g., Pedlosky 1987, 
chapter 3): 



Dh 
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where 
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Dt 



-gVh - /k X V 
-hV-v, 
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(8a) 
(8b) 

(8c) 



Forcing and dissipation are not included in Eq. (8), but they 
can be added in the usual way. In the absence of forcing and 
dissipation, the equations preserve the potential vorticity, 



C + f 

h ' 



(9) 



following the flow. 

If Eq. (8) is derived as the vertical mean of the flow of an 
isentropic atmosphere with a free upper boundary, h must 
be replaced by h'^ in the geopotential gradient term. If they 
are derived as a vertical mean of the flow of an isentropic 
atmosphere between rigid upper and lower boundaries, h 
must be replaced by in the geopotential gradient 

term. Note that while V • v ^ in Eq. (8b), V • u = 0, 
since the layer is homogeneous (i.e., density is constant). 
Hence, the sound speed Cg oo, and the sound waves 
are filtered out from the system. However, the system does 
retain gravity waves, which propagate at speed Cg = y/gh. 

The shallow-water equations are widely used as a pro- 
cess model in geophysical fluid dynamics. They are much 
simpler than the primitive equations, yet they still describe 
a wealth of phenomena — including vortices, jet streams, 
Rossby waves, gravity waves, and the interactions between 
them. Excluded are any processes that depend on the de- 
tails of the vertical structure — including vertically propa- 
gating waves, baroclinic instabilities (see §3.7), and depth- 
dependent flow. However, because both rotational and 
buoyancy processes are included (the latter via the variable 
layer thickness), the shallow-water model — as well as all 
the models discussed so far — includes a fundamental length 
scale called the Rossby radius of deformation (often sim- 
ply called the deformation radius), which is a natural length 
scale for a variety of phenomena that depend on both ro- 
tation and stratification. In the shallow-water system, this 
length scale is 



f 



(10) 
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2.4. TVo-Dimensional, Nondivergent Model 

This is the simplest useful one-layer model for large- 
scale dynamics. For large-scale weather systems character- 
ized by U/cg <C 1, we can apply a rigid upper boundary to 
the shallow-water model, since Cg represents external grav- 
ity wave speed in the model. Then, H is large and Eq. (8b) 
imphes V-v <C 1. Taking V • v = then gives the 2D 
nondivergent equation: 

^ = -gV/l - /k X V, (11) 

where D/Dt is same as in (8c). The V • v = restriction 
on the velocity implies that we can define a streamfunction, 
tp{x,t), such that 

dip dip^ 
dy ' dx 



(12) 



Using this definition, Eq. (11) can finally be written as 
a single governing equation for the evolution of the stream- 
function: 

-g(V2v + /)=0. (13) 

From Eq. (13), we see that = V^i/' -I- / is the materi- 
ally conserved potential vorticity for the 2D nondivergent 
model. This also results simply by letting h — » constant in 
Eq. (9). 

In addition to the nonlinear vorticity advection, this 
equation — along with all the other equation sets described 
in this section — represents the dynamical effects of latitudi- 
nally varying Coriolis parameter. This is the so-called "beta 
effect," where (3 = df /dy is the northward gradient of the 
CorioUs parameter. 

Eq. (13) describes Rossby waves, non-linear advec- 
tion, and phenomena — such as the formation of zonal jet 
streams — ^that require the interaction of all these aspects 
(see §3.6). However, it lacks a finite deformation radius 
{Ld ^ oo) and does not possess gravity wave solutions. 
Therefore, any phenomena that depend on finite deforma- 
tion radius, gravity waves, or buoyancy cannot be cap- 
tured. (These assumptions render the equation valid only 
for U/cg ^ 1 and L/Ld <C 1.) As a result, the 2D non- 
divergent model cannot serve as an accurate predictive tool 
for most apphcations; however, its very simplicity renders 
it a valuable process model for investigating jet formation 
in the simplest possible setting. For a review of examples, 
see for example Vasavada and Showman (2005) or Vallis 
(2006). 

2.5. Conserved Quantities 

Potential vorticity conservation has been emphasized 
throughout because of its central importance in atmospheric 
dynamics. There are other useful conserved quantities. For 
example, the full Navier-Stokes equation gives 

D , ^ f I 



Dt 



(rxu) 



r X 



- 2n X u - V$* 



where is effective geopotential and F represents any ad- 
ditional forces on the fluid. From this, we obtain the con- 
servation law for specific, axial angular momentum M: 



DM 

~1h 



1 dp 
pdX 



+ J^x cos (f), 



where 



A4 = {fir cos 4> + u)r cos < 



(15a) 



(15b) 



Eq. (15) relates the material change of A4 to the axial com- 
ponents of torques present. For a thin atmosphere, r can be 
replaced with a. 

Eq. (1) also gives the material conservation law for the 
specific total energy E: 



DE 

15i 



-y-{pu) + (g„et + u-F), 



(16) 



where E is the total energy including kinetic, potential, and 
internal contributions: E = ^u^ + $ + c„T with c„ the 
specific heat at constant volume and T the temperature. In 
flux form, the conservation law is: 



d_ 
dt 



(pE) + V- [{pE + p)u] = p (Qnet + u-F). (17) 



As already noted for A4, the total energy reduces in the 
appropriate way for the various simpler physical situations 
discussed in previous subsections. For example, $ and CyT 
terms do not exist for the 2D nondivergent case. An im- 
portant issue in the study of atmospheric energetics is the 
extent to which $ and CyT are available to be converted to 

2" • 



3. BASIC CONCEPTS 

The equation sets summarized in §2 describe nonhnear, 
potentially turbulent flows with many degrees of freedom. 
Unfortunately, due to the nonlinearity and complexity, an- 
alytic solutions rarely exist, and one must resort to solving 
the equations numerically on a computer. To represent the 
atmospheric circulation of a particular planet, the chosen 
equation set is solved numerically with a specified spatial 
resolution and timestep, subject to appropriate parameter 
values (e.g., composition, gravity, planetary rotation rate), 
boundary conditions, and forcing/damping (e.g., prescrip- 
tions for heating/cooling and friction). 

Such models vary greatly in complexity and numerical 
method. General Circulation Models (GCMs) in the Solar- 
System studies Uterature, for example, typically solve the 
3D primitive equations with sophisticated representations 
of radiative transfer, cloud formation, surface/atmosphere 
interactions, surface ice formation, and (if relevant) oceanic 
processes. These models are useful for exploring the inter- 
action of dynamics with surface processes, radiation, and 
cUmate and are needed for quantitative comparisons with 
observational records. 



(14) 
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However, because of their complexity, numerical simu- 
lations with full GCMs are computationally expensive, lim- 
iting such simulations to only moderate spatial resolution 
and making it difficult to broadly survey the relevant param- 
eter space. Even more problematic, because of the inherent 
complexity of nonlinear fluid dynamics and its possible in- 
teractions with radiation and surface processes, it is rarely 
obvious why a given GCM simulation produces the output it 
does. By itself, the output of a sophisticated 3D model often 
provides Uttle more fundamental understanding than the ob- 
servations of the actual atmosphere themselves. To under- 
stand how a given atmospheric circulation would vary under 
different planetary parameters, for example, an understand- 
ing of the mechanisms shaping the circulation is required. 
Although careful diagnostics of GCM results can provide 
important insights into the mechanisms that are at play, a 
deep mechanistic understanding does not always flow natu- 
rally from such simulations. 

Rather, obtaining a robust understanding requires a di- 
versity of model types, ranging from simple to complex, 
in which various processes are turned on and off and the 
results carefully diagnosed. This is called a modeling hier- 
archy and its use forms the backbone of forward progress 
in the field of atmospheric dynamics of Earth and other 
Solar-System planets (see, e.g.. Held 2005). For exam- 
ple, despite the existence of numerous full GCMs for mod- 
ern Earth climate, significant advances in our understand- 
ing of the mechanisms shaping the atmospheric circulation 
rely heavily on the usage of linear models, simplified one- 
layer non-Unear models (such as the 2D non-divergent or 
shallow-water models), and 3D models that do not include 
the sophisticated treatments of radiation and sub-gridscale 
convective processes included in full GCMs.^ Even more 
fundamentally, obtaining understanding requires the devel- 
opment of basic theory that can (at least qualitatively) ex- 
plain the results of these various models as well as obser- 
vations of actual atmospheres. One of the major goals in 
performing simplified models is to aid in the construction 
of such a theory (see, e.g., Schneider 2006). 

Exoplanet GCMs will surely be useful in the coming 
years. But, as with Solar-System planets, we expect that 
a fundamental understanding will require use of a modehng 
hierarchy as well as basic theory. In this section we survey 
key concepts in atmospheric dynamics that provide insight 
into the expected atmospheric circulation regimes. Empha- 
sis is placed on presenting a conceptual understanding and 
as such we describe not only GCM results but basic theory 
and the results of highly simplified models as well. Here 
we focus on basic aspects relevant to both gaseous and ter- 
restrial planets. Detailed presentations of issues specific to 
giant and terrestrial exoplanets are deferred to §4 and §5. 



'^To illustrate, a summary of the results of such a hierarchy for understand- 
ing Jupiter's jet streams can be found in Vasavada and Showman (2005). 



3.1. Energetics of atmospheric circulation 

Atmospheric circulations involve an energy cycle. Ab- 
sorption of starlight and emission of infrared energy to 
space creates potential energy, which is converted to kinetic 
energy and then lost via friction. Each step in the process 
involves nonlinearities, and generally the atmosphere self- 
adjusts so that, in a time mean sense, the conversion rates 
balance. 

What matters for driving the circulation is not the to- 
tal potential energy but rather the fraction of the potential 
energy that can be extracted by adiabatic atmospheric mo- 
tions. For example, a stably stratified, horizontally uniform 
atmosphere can contain vast potential energy, but none can 
be extracted — any adiabatic motions can only increase the 
potential energy of such a state. Uniformly heating the top 
layers of such an atmosphere would further increase its po- 
tential energy but would still preclude an atmospheric cir- 
culation. 

For convecting atmospheres, creating extractable poten- 
tial energy requires heating the fluid at lower altitudes than 
it is cooled. This creates buoyant air parcels (positively 
buoyant at the bottom, negatively buoyant at the top); verti- 
cal motion of these buoyant parcels releases potential en- 
ergy and drives convection.** But, most atmospheres are 
stably stratified, and in this case extractable energy — called 
available potential energy — only exists when density varies 
horizontally on isobars (Peixoto and Oort 1992, chapter 14). 
In this case, the denser regions can slide laterally and down- 
ward underneath the less-dense regions, decreasing the po- 
tential energy and creating kinetic energy (winds). Con- 
tinual generation of available potential energy (required to 
balance its continual conversion to kinetic energy and loss 
via friction) requires heating the regions of the atmosphere 
that are already hot (e.g., the tropics on Earth) and cooling 
the regions that are already cold (e.g., the poles). For Earth, 
available potential energy is generated at a global-mean rate 
of ~ 2 W m~^, which is ~ 1% of the global- mean absorbed 
and radiated flux of 240 W m~^ (Peixoto and Oort 1992, 
pp. 382-385). 

The rate of frictional dissipation can affect the mean 
state, but rigorously representing such friction in models is 
difficult. For Solar-System planets, kinetic -energy loss oc- 
curs via turbulence, waves, and friction against the surface 
(if any). Ohmic dissipation may be important in the deep 
interiors of gas giants (Kirk and Stevenson 1987; Liu et al. 
2008), as well as in the upper atmosphere where ioniza- 
tion becomes important. These processes sometimes have 
length scales much smaller (by up to several orders of mag- 
nitude) than can easily be resolved in global, 3D numer- 
ical models. In Earth GCMs, such frictional dissipation 
mechanisms are therefore often parameterized by adding to 



'To emphasize the importance of the distinction, consider a hot, isolated 
giant planet. The cooling caused by its radiation to space decreases its 
total potential energy, yet (because the cooling occurs near the top) this in- 
creases the fraction of the remaining potential energy that can be extracted 
by motions. This is what can allow convection to occur on such objects. 
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the equations quasi-empirical damping terms (e.g., a verti- 
cal diffusion to represent turbulent kinetic-energy losses by 
small-scale shear instabilities and breaking waves). A diffi- 
culty is that such prescriptions, while physically motivated, 
are often non-rigorous and the extent to which they can 
be extrapolated to other planetary environments is unclear. 
Perhaps for this reason, models of hot Jupiters published 
to date do not include such parameterizations of frictional 
processes (although they all include small-scale viscosity 
for numerical reasons).' Nevertheless, Goodman (2009) 
has highlighted the possible importance that such processes 
could play in the hot- Jupiter context, and future models of 
hot Jupiters will surely explore the possible effect that fric- 
tion may have on the mean states. 

Solar-System planets offer interesting lessons on the 
role of friction. Despite absorbing a greater solar flux 
than any other thick atmosphere in our Solar System, 
Earth's winds are relatively slow, with a mean wind speed 
of ~ 20msec~^. In contrast, Neptune absorbs a solar 
flux only 0.1% as large, but has wind speeds reaching 
400 m see" ^. Presumably, Neptune can achieve such fast 
winds despite its weak radiative forcing because its fric- 
tional damping is extremely weak. Qualitatively, this makes 
sense because Neptune lacks a surface, which is a primary 
source of frictional drag on Earth. More puzzling is the fact 
that Neptune has significantly stronger winds than Jupiter 
(Table 1) despite absorbing only 4% the solar flux absorbed 
by Jupiter. Possible explanations are that Jupiter experi- 
ences greater frictional damping than Neptune or that it has 
equilibrated to a state that has relatively slow wind speeds 
despite weak damping. This is not well understood and 
argues for humility in efforts to model the circulations of 
exoplanets. 

3.2. Timescale arguments for the coupled radiation- 
dynamics problem 

The atmospheric circulation represents a coupled radiation- 
hydrodynamics problem. The circulation advects the tem- 
perature field and thereby influences the radiation field; in 
turn, the radiation field (along with atmospheric opacities 
and surface conditions) determines the atmospheric heating 
and cooling rates that drive the circulation. Rigorously at- 
tacking this problem requires coupled treatment of both ra- 
diation and dynamics. However, crude insight into the ther- 
mal response of an atmosphere can be obtained with sim- 
ple timescale arguments. Suppose Tadvect is an advection 
time (e.g., the characteristic time for air to advect across a 
hemisphere) and Trad is the radiative time (i.e., the charac- 
teristic time for radiation to induce large fractional entropy 
changes). When Trad ^ ^advoct, we expect temperature to 
deviate only slightly from the (spatially varying) radiative 
equilibrium temperature structure. Because the radiative- 



^Note that a statistically steady (or quasi-steady) state can still occur in 
such a case; this requires the atmosphere to self-adjust so that the rates 
of generation of available potential energy and its conversion to kinetic 
energy become small. 



equilibrium temperature typically varies greatly from day- 
side to nightside (or from equator to pole), this implies that 
such a planet would exhibit large fractional temperature 
contrasts. On the other hand, when Trad ^ Tadvect, dynam- 
ical transport dominates and air will tend to homogenize its 
entropy, implying that lateral temperature contrasts should 
be modest. 

In estimating the advection time, one must distinguish 
north-south from east-west advection; east-west advection 
(relative to the pattern of steUar insolation) will often be 
dominated by the planetary rotation. For synchronously ro- 
tating planets, a characteristic horizontal advection time is 



'^advect 



a 

U' 



(18) 



where C/ is a characteristic horizontal wind speed. A simi- 
larly crude estimate of the radiative time can be obtained by 
considering a layer of pressure thickness Ap that is slightly 
out of radiative equilibrium and radiates to space as a black- 
body. If the radiative equilibrium temperature is Trad and 
die actual temperature is Trad + AT, with AT < Trad, 
then the net flux radiated to space is AaT^^AT and the 
radiative timescale is (Showman and Guillot 2002; James 
1994, pp. 65-66) 



'Trad 



Ap 



g 4(tT3' 



(19) 



In deep, optically thick atmospheres where the radiative 
transport is diffusive, a more appropriate estimate might be 
a diffusion time, crudely given by Trad ~ H'^ /D, where H 
is the vertical height of a thermal perturbation and D is the 
radiative diffusivity. 

Showman et al. (2008b) estimated advective and radia- 
tive time constants for Solar-System planets and found that, 
as expected, planets with Trad ^ ^advect generally have 
small horizontal temperature contrasts and vice versa. 

For hot Jupiters, most models suggest peak wind speeds 
of ^l-3kmscc~^ (§4.3), implying advection times of 
~10^sec based on the peak speed. Eq. (19) would then 
suggest that Trad < Tadvect at p < 1 bar whereas 
Trad 3> Tadvoct at p 1 bar Thus, one might crudely ex- 
pect large day-night temperature differences at low pressure 
and small day-night temperature difference at high pres- 
sure, with the transition occurring at ^ 0.1-1 bar. These 
estimates are generally consistent with the observational 
inference of Barman (2008) and 3D numerical simulations 
(e.g.. Showman et al. 2009; Dobbs-Dixon and Lin 2008) 
of hot Jupiters — though some uncertainties still exist with 
modehng and interpretation. 

For synchronously rotating terrestrial planets in the hab- 
itable zones of M dwarfs, a mean wind speed of 20 m sec" ^ 
(typical for terrestrial planets in our Solar System; see Ta- 
ble 1) would imply an advection time of ^3 Earth days. 
For a temperature of 300 K, Eq. (19) would then imply that 
Trad is much smaller (greater) than Tadvect when the sur- 
face pressure is much less (greater) than ^0.2 bars. This 
argument suggests that synchronously rotating terrestrial 
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exoplanets with a surface pressure much less than -^0.2 
bars should develop large day-night temperature differ- 
ences, whereas if the surface pressure greatly exceeds '~0.2 
bars, day-night temperature differences would be modest. 
As with hot Jupiters, these estimates are consistent with 3D 
GCM simulations (Joshi et al. 1997), which suggest that 
this transition should occur at ~0.1bars. These estimates 
may have relevance for whether CO2 atmospheres would 
collapse due to nightside condensation and hence whether 
such planets are habitable. 

3.3. Basic force balances: importance of rotation 

Planets rotate, and this typically constitutes a dominant 
factor in shaping the circulation. The importance of rota- 
tion can be estimated by performing a scale analysis on the 
equation of motion. Suppose the circulation has a mean 
speed U and that we are interested in flows with character- 
istic length scale L (this might approach a planetary radius 
for global-scale flows). To order-of-magnitude, the strength 
of the acceleration term is C/^/L (namely, U divided by a 
time L/U to advect fluid across a distance L), while the 
magnitude of the Coriolis term is fU. The ratio of the ac- 
celeration term to the CorioUs term can therefore be repre- 
sented by the Rossby number. 



Ro 



U 



(20) 



Whenever Ro <^ 1, the acceleration terms Dv/Dt are 
weak compared to the Coriolis force per unit mass in the 
horizontal momentum equation. Because friction is gener- 
ally weak, the only other term that can balance the horizon- 
tal Coriolis force is the pressure-gradient force, which is 
just — Vp$ in pressure coordinates. The resulting balance, 
caUed geostwphic balance, is given by 



fu- 



\dy 



(21) 



where x and y are eastward and northward distance, respec- 
tively and the derivatives are evaluated at constant pressure. 
In our Solar System, geostrophic balance holds at large 
scales in the mid- and high-latitude atmospheres of Earth, 
Mars, Jupiter, Saturn, Uranus, and Neptune. Rossby num- 
bers range from 0.01-0.1 for these rapidly rotating planets, 
but exceed unity in the stratosphere of Titan^" and reach 
^10 for Venus, implying in the latter case that the Cori- 
olis force plays a less important role in the force balance 
(Table 1). Note that, even on rapidly rotating planets, hor- 
izontal geostrophy breaks down at the equator, where the 
horizontal Coriolis forces go to zero. 

Determining Rossby numbers for exoplanets requires es- 
timates of wind speeds, which are unknown. Some models 
of hot Jupiter atmospheres have generally suggested peak 
winds of several km sec" ^ near the photosphere,'^ with 



^"Titan's Rossby number is smaller near the surface, where winds are weak. 
11 Defined here as the approximate pressure at which infrared photons can 
escape directly to space. 



mean values perhaps a factor of several smaller. To illus- 
trate the possibilities, Table 1 presents Ro values for sev- 
eral hot Jupiters assuming a range of wind speeds of 100- 
4000 m sec" ^. Generally, if mean wind speeds are fast 
(several km sec" ^), Rossby numbers approach or exceed 
unity. If mean wind speeds are hundreds of msec"^ or 
less, Rossby numbers should be much less than one, imply- 
ing that geostrophy approximately holds. One might thus 
plausibly expect a situation where the Coriolis force plays 
an important but not overwhelming role (i.e. Ro ~ 1) near 
photosphere levels, with the flow transitioning to geostro- 
phy in the interior if winds are weaker there. 

Geostrophy implies that, rather than flowing from pres- 
sure highs to lows as often occurs in a non-rotating fluid, 
the primary horizontal wind flows perpendicular to the hor- 
izontal pressure gradient. Thus, the primary flow does not 
erase the horizontal pressure contrasts; rather, the Cori- 
olis forces associated with that flow actually help pre- 
serve large-scale pressure gradients in a rotating atmo- 
sphere. Geostrophy explains why the isobar contours in- 
cluded in most mid-latitude weather maps (say of the U.S. 
or Europe) are so useful: the isobars describe not only the 
pressure field but the large-scale wind field, which flows 
along the isobar contours. 

In many cases, rotation inhibits the ability of the circu- 
lation to equalize horizontal temperature differences. On 
rapidly rotating planets hke the Earth, which is heated by 
sunlight primarily at low latitudes, the mean horizontal tem- 
perature gradients in the troposphere'^ generally point from 
the poles toward the equator. Integration of the hydrostatic 
equation (Eq. 5b) impUes that the mean pressure gradients 
also point north-south. Thus, on a rapidly rotating planet 
where geostrophy holds and the primary temperature con- 
trast is between equator and pole, the mean midlatitude 
winds will be east-west rather than north-south — thus lim- 
iting the ability of the circulation to homogenize its tem- 
perature differences in the north-south direction. We might 
thus expect that, everything else being equal, a more rapidly 
rotating planet will harbor a greater equator-to-pole temper- 
ature difference. 

In rapidly rotating atmospheres, a tight link exists be- 
tween horizontal temperature contrasts and the vertical gra- 
dients of the horizontal wind. This can be shown by taking 
the derivative with pressure of Eq. (21) and invoking the 
hydrostatic balance equation (5b) and the ideal-gas law. We 
obtain the thermal-wind equation for a shallow atmosphere 
(Holton 2004, pp. 70-75): 



/ 



du 



d{RT) 
dy 



f 



dv 
dlnp 



d{RT) 
dx 



(22) 



^^The troposphere is the bottommost, optically thick layer of an atmosphere, 
where temperature decreases with altitude and convection may play an 

important role; the tropopause defines the top of the troposphere, and the 
stratosphere refers to the stably stratified, optically thin region overlying 
the troposphere. In some cases, a stratosphere's temperature may increase 
with altitude due to absorption of sunlight by gases or aerosols; in other 
cases, however, the stratosphere's temperature can be nearly constant with 
altitude. 
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Table 1 
Planetary parameters 



Planet 


o* 


Rotation period' 


n 


gravity^ 








Ui 




Lo/a* 






(IQS km) 


(Earth days) 


(radsec"^) 


(msec-^) 


(Wm-2) 


(K) 


(km) 


(msec-^) 








Venus 


6.05 


243 


3 X IQ-'' 


8.9 


2610 


232 


5 


~ 20 


10 


70 


7 


Earth 


6.37 


1 


7.27 X 10-6 


9.82 


1370 


255 


7 


~ 20 


0.1 


0.3 


0.5 


Mars 


3.396 


1.025 


7.1 X 10-5 


3.7 


590 


210 


11 


~ 20 


0.1 


0.6 


0.6 


Titan 


2.575 


16 


4.5 X 10-s 


1.4 


15 


85 


18 


~ 20 


2 


10 


3 


Jupiter 


71.4 


0.4 


1.7 X 10-4 


23.1 


50 


124 


20 


~40 


0.02 


0.03 


0.1 


Saturn 


60.27 


0.44 


1.65 X 10-* 


8.96 


15 


95 


39 


~ 150 


0.06 


0.03 


0.3 


Uranus 


25.56 


0.72 


9.7 X 10-5 


8.7 


3.7 


59 


25 


~ 100 


0.1 


0.1 


0.4 


Neptune 


24.76 


0.67 


1.09 X 10-"' 


11.1 


1.5 


59 


20 


~ 200 


0.1 


0.1 


0.6 


WASP- 12b 


128 


1.09 


6.7 X 10-5 


11.5 


8.8 X 10"^ 


2500 


800 




0.01-0.3 


0.1 


0.2-1.5 


HD 189733b 


81 


2.2 


3.3 X 10-5 


22.7 


4.7 X 105 


1200 


200 




0.03-1 


0.3 


0.4-3 


HD 149026b 


47 


2.9 


2.5 X 10-5 


21.9 


1.8 X 10'^ 


1680 


280 




0.06-2 


0.8 


0.6-4 


HD 209458b 


94 


3.5 


2.1 X 10-5 


10.2 


1.0 X 10*^ 


1450 


520 




0.04-1 


0.4 


0.5-3 


TrES-2 


87 


2.4 


2.9 X 10-5 


21 


1.1 X 10*^ 


1475 


260 




0.03-1 


0.3 


0.4-3 


TrES-4 


120 


3.5 


2.0 X 10-5 


7.8 


2.5 X 10'^ 


1825 


870 




0.03-1 


0.4 


0.4-3 


HAT-P-7b 


97 


2.2 


3.3 X 10-5 


25 


4.7 X 10"^ 


2130 


320 




0.02-1 


0.3 


0.4-3 


GJ 436b 


31 


2.6 


2.8 X 10-5 


9.8 


4.3 X 10" 


660 


250 




0.1-3 


0.7 


0.8-5 


HAT-P-2b 


68 


5.6 


1.3 X 10-5 


248 


9.5 X 105 


1400 


21 




0.1-3 


1 


0.8-5 


Corot-Exo-4b 


85 


9.2 


7.9 X 10-s 


13.2 


3.0 X 105 


1080 


300 




0.1^ 


1 


0.9-5 



Note. — *Equatorial planetary radius. "Assumes synchronous rotation for exoplanets. Equatorial gravity at the surface. '-'Mean incident stellar flux. *Global- 
average blackbody emission temperature, which for exoplanets is calculated from Eq. (51) assuming zero albedo. ^Pressure scale height, evaluated at temperature Te. 
■'•Rough estimates of characteristic horizontal wind speed. Estimates for Venus and Titan are in the high-altitude superrotating jet; both planets have weaker winds (few 
msec-^) in the bottom scale height. In all cases, peak winds exceed the listed values by factors of two or more. ^Rossby number, evaluated in mid-latitudes using 
wind values listed in Table and L ~ 2000 km for Earth, Mars, and Titan, 6000 km for Venus, and 10* km for Jupiter, Saturn, Uranus, and Neptune. For exoplanets, we 
present a range of possible values evaluated with L = a and winds from 100 to 4000 msec- ^. *Ratio of Rossby deformation radius to planetary radius, evaluted in 
mid-latitudes with H equal to the pressure scale height and appropriate for a vertically isothermal temperature profile. * Ratio of Rhines length (Eq. 35) to planetary 
radius, calculated using the equatorial value of /? and the wind speeds listed in the Table. 
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where R is the specific gas constant (i.e., the universal gas 
constant divided by the molar mass). The equation states 
that, in a geostrophically balanced atmosphere, north-south 
temperature gradients must be associated with a vertical 
gradient in the zonal (east-west) wind, whereas east-west 
temperature gradients must be associated with a vertical 
gradient in the meridional (north-south) wind. Given the 
primarily equatorward pointing midlatitude horizontal tem- 
perature gradient in the tropospheres of Earth and Mars, 
for example, and given the weak winds at the surface of a 
terrestrial planet (a result of surface friction), this equation 
correctly demonstrates that the mean mid-latitude winds in 
the upper troposphere must flow to the east — as observed 
for the mid-latitude tropospheric jet streams on Earth and 
Mars. 

Geostrophy relates the 3D structure of the winds and 
temperatures at a given time but says nothing about the 
flow's time evolution. In a rapidly rotating atmosphere, 
both the temperatures and winds often evolve together, 
maintaining approximate geostrophic balance as they do 
so. This time evolution depends on the ageostrophic com- 
ponent of the circulation, which tends to be of order Ro 
smaller than the geostrophic component. The fact that a 
time evolving flow maintains approximate geostrophic bal- 
ance imphes that, in a rapidly rotating atmosphere, adjust- 
ment mechanisms exist that tend to re-estabhsh geostrophic 
balance when departures from it occur. 

What is the mechanism for establishing and maintaining 
geostrophy? If a rapidly rotating atmosphere deviates from 
geostrophic balance, it imphes that the horizontal pressure- 
gradient and Coriolis forces only imperfectly cancel, leav- 
ing an unbalanced residual force. This force generates a 
component of ageostrophic motion between pressure highs 
and lows. The Coriolis force on this ageostrophic motion, 
and the alteration of the pressure contrasts caused by the 
ageostrophic wind, act to re-establish geostrophy. 

To give a concrete example, imagine an atmosphere with 
zero winds and a localized circular region of high surface 
pressure surrounded on all sides by lower surface pres- 
sure. This state has an unbalanced pressure-gradient force, 
which would induce a horizontal acceleration of fluid ra- 
dially away from the high-pressure region. The horizontal 
Coriohs force on this outward motion causes a lateral de- 
flection (to the right in the northern hemisphere and left in 
the southern hemisphere^^), leading to a vortex surrounding 
the high-pressure region. The Coriohs force on this vor- 
tical motion points radially toward the high-pressure cen- 
ter, resisting its lateral expansion. This process continues 
until the inward-pointing Coriolis force balances the out- 
ward pressure-gradient force — hence establishing geostro- 
phy and inhibiting further expansion. Although this is an 
extreme example, radiative heating/coohng, friction, and 
other forcings graduaUy push the atmosphere away from 



i^Northem and southern hemispheres are here defined as those where fl • 
k > and < 0, respectively, where CI and k are the planetary rotation 
vector and local vertical (upward) unit vector. 



geostrophy, and the process described above re-establishes 
it. 

This process, called geostrophic adjustment, tends to oc- 
cur with a natural length scale comparable to the Rossby 
radius of deformation, given in a 3D system by 



NH 



f 



(23) 



where is Brunt- Vaisala frequency (Eq. 6), H is the verti- 
cal scale of the flow, and / is the Coriohs parameter. Thus, 
geostrophic adjustment naturally generates large-scale at- 
mospheric flow structures with horizontal sizes compara- 
ble to the deformation radius (see Holton (2004) or Valhs 
(2006) for a detailed treatment). 

Equation (22), as written, applies to atmospheres that 
are vertically thin compared to their horizontal dimensions 
(i.e., shallow atmospheres). However, a non- shallow analog 
of Eq. (22) can be obtained by considering the 3D vortic- 
ity equation. In a geostrophically balanced fluid where the 
friction force is weak, the vorticity balance is given by (see 
Pedlosky 1987, p. 43) 



2(fl- V)u-2rj(V-u) 



(24) 



where n is the planetary rotation vector and u is the 3D 
wind velocity. The term on the right, called the baroclinic 
term, is nonzero when density varies on constant-pressure 
surfaces. In the interior of a giant planet, however, convec- 
tion tends to homogenize the entropy, in which case frac- 
tional density variations on isobars are extremely small. 
Such a fluid, where surfaces of constant p and p align, 
is called a barotropic fluid. In this case the right side of 
Eq. (24) can be neglected, leading to the compressible-fluid 
generalization of the Taylor-Proudman theorem: 



2(J2 • V)u - 2J2(V • u) = 



(25) 



Consider a Cartesian coordinate system (x«, y*, z») with 
the axis parallel to and the .x* and axes lying in the 
equatorial plane. Eq. (25) can then be expressed in compo- 
nent form as 

= (26) 



dxt, dy^ 







(27) 



where and are the wind components along the 
and axes, respectively. These equations state that the 
wind components parallel to the equatorial plane are in- 
dependent of direction along the axis perpendicular to the 
equatorial plane, and moreover that the divergence of the 
winds in the equatorial plane must be zero. The flow thus 
exhibits a structure with winds constant on columns, called 
Taylor columns, that are parallel to the rotation axis. In 
the context of a giant planet, motion of Taylor columns to- 
ward or away from the rotation axis is disallowed because 
the planetary geometry would force the columns to stretch 
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or contract, causing a nonzero divergence in the equatorial 
plane — violating Eq. (27). Rather, the columns are free 
to move only along latitude circles. The flow then takes 
the form of concentric cylinders, centered about the rota- 
tion axis, which can move in the east-west direction. This 
columnar flow provides one model for the structure of the 
winds in the deep interiors of Jupiter, Saturn, Uranus, Nep- 
time, and giant exoplanets (see §4). 

3.4. Other force balances: The case of slowly rotating 
planets 

The pressure-gradient force can be balanced by forces 
other than the Coriolis force. An important example is cy- 
clostrophic balance, which is a balance between horizon- 
tal centrifugal and pressure-gradient force, expressed in its 
simplest form as: 



(28) 



ui _ldp 

r p dr 

where we are considering circular flow around a central 
point. Here, Ut is the tangential speed of the circular flow 
and r is the radius of curvature of the flow. The left-hand 
side of the equation is the centrifugal force per unit mass 
and the right-hand side is the radial pressure-gradient force 
per unit mass. Cyclostrophic balance is the force balance 
that occurs within dust devils and tornadoes, for example. 

In the context of a global-scale planetary circulation, if 
the primary flow is east-west, the centrifugal force mani- 
fests as the curvature term v? tamp /a, where a is planetary 
radius and cf) is latitude (see Holton 2004, pp. 31-38 for a 
discussion of curvature terms). Cyclostrophic balance can 
then be written 



(29) 



This is the force balance relevant on Venus, for exam- 
ple, where a strong zonal jet with peak speeds reach- 
ing ~100msec~^ dominates the stratospheric circulation 
(Gierasch et al. 1997). Eq. (29) can be differentiated with 
respect to pressure to yield a cyclostrophic version of the 
thermal-wind equation: 

d\np i&ncj) dy ' 

where hydrostatic balance and the ideal-gas law have been 
invoked. 

Thus, on a slowly rotating, cyclostrophically balanced 
planet like Venus, Eq. (30) would indicate that varia- 
tion of the zonal winds with height requires latitudinal 
temperature gradients as occur with geostrophic balance. 
Note, however, that for zonal winds whose strength in- 
creases with height, the latitudinal temperature gradients 
for cyclostrophic balance point equatorward regardless of 
whether the zonal winds are east or west. This differs from 
geostrophy, where latitudinal temperature gradients would 
point equatorward for a zonal wind that becomes more 
eastward with height but poleward for a zonal wind that 
becomes more westward with height. 



3.5. What controls vertical velocities? 

In atmospheres, mean vertical velocities associated with 
the large-scale circulation tend to be much smaller than hor- 
izontal velocities. This results from several factors: 

Large aspect ratio: Atmospheres generally have horizon- 
tal dimensions greatly exceeding their vertical dimensions, 
which leads to a strong geometric constraint on vertical mo- 
tions. To illustrate, suppose the continuity equation can 
be approximated in 3D by the incompressibility condition 
V ■ u = du / dx + dv / dy + dw/dz = 0. To order of mag- 
nitude, the horizontal terms can be approximated as U/L 
and the vertical term can be approximated as W/ H, where 
L and H are the horizontal and vertical scales and U and 
W are the characteristic magnitudes of the horizontal and 
vertical wind. This then suggests a vertical velocity of ap- 
proximately 

W -^^U (31) 

On a typical terrestrial planet, U ~ lOmsec"^, L ~ 
1000 km, and H ~ 10 km, which would suggest W 
O.lmsec"^, two orders of magnitude smaUer than hori- 
zontal velocities. However, this estimate of W is an upper 
limit, since partial cancellation can occur between du/dx 
and dv/dy, and indeed for most atmospheres Eq. (31) 
gready overestimates their mean vertical velocities. On 
Earth, for example, mean midlatitude vertical velocities on 
large scales (^lO'' km) are actually ~10~^ msec~^ in the 



troposphere and ^10 



in the stratosphere, much 



smaller than suggested by Eq. (31). 

Suppression of vertical motion by rotation: The geostrophic 
velocity flows perpendicular to the horizontal pressure gra- 
dient, which implies that a large fraction of this flow can- 
not cause horizontal divergence or convergence (as neces- 
sary to allow vertical motions). However, the ageostrophic 
component of the flow, which is 0{Ro) smaller than 
the geostrophic component, can cause horizontal conver- 
gence/divergence. Thus, to order-of-magnitude one might 
expect W ~ RoU/L. Considering the geostrophic flow 
and taking the curl of Eq. (21), we can show that the hori- 
zontal divergence of the horizontal geostrophic wind is 



Vp-VG 



P _ V 



(32) 



where in the rightmost expression has been evalu- 
ated assuming spherical geometry. To order of magnitude, 
V ^ U. Eq. (32) then impUes that to order-of-magnitude 
the horizontal flow divergence is not U/L but rather U/a. 
On many pkinets, the dominant flow structures are smaller 
than the planetary radius (e.g., L/a ^ 0.1-0.2 for Earth, 
Jupiter, and Saturn). This implies that the horizontal di- 
vergence of the geostrophic flow is ^L/a smaller than the 
estimate in Eq. (31). Taken together, these constraints sug- 
gest that rapid rotation can suppress vertical velocities by 
close to an order of magnitude. For Earth parameters, the 
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estimates imply W 0.01msec ^, similar to the mean 
vertical velocities in Earth's troposphere. 

Suppression of vertical motion by stable stratification: 
Most atmospheres are stably stratified, implying that en- 
tropy (potential temperature) increases with height. Adia- 
batic expansion/contraction in ascending (descending) air 
would cause temperature at a given height to decrease (in- 
crease) over time. In the absence of radiation, such steady 
flow patterns are unsustainable because they induce density 
variations that resist the motion — ascending air becomes 
denser and descending air becomes less dense than the 
surroundings. Thus, in stably stratified atmospheres, the 
radiative heating/cooling rate exerts major control over the 
rate of vertical ascent/decent. Steady vertical motion can 
only occur as fast as radiation can remove the temperature 
variations caused by the adiabatic ascent/descent (Show- 
man et al. 2008a), when conduction is negligible. The idea 
can be quantified by rewriting the thermodynamic energy 
equation (Eq. 5d) in terms of the Brunt- Vaisala frequency: 



dT 



■ V • VpT - u) 



Rp 



Qnet 



(33) 



where Hp is the pressure scale height. In situations where 
the flow is approximately steady and vertical thermal advec- 
tion is more important than (or comparable to) horizontal 
thermal advection, an estimate of the magnitude of the ver- 
tical velocity in an overturning circulation can be obtained 
by equating the right side to the last term on the left side. 
Converting to velocity expressed in height units, we obtain 
(Showman et al. 2008a; Showman and GuiUot 2002)" 



W 



Qnet R 



(34) 



For Earth's midlatitude troposphere. Hp ^ 10 km, N ^ 
0.01 sec~^, and Qnet/cp ~ 3 x 10~^Ksec^^, implying 
W ~ 10~^ m sec~^. In the stratosphere, however, the heat- 
ing rate is lower, and the stable stratification is greater, lead- 
ing to values of W 



10 ^msec ^. 



For a canonical hot Jupiter with a 3-day rotation pe- 
riod, horizontal wind speeds of km sec" ^ imply Rossby 
numbers of ^1, and for scale heights of ^300 km (Ta- 
ble 1) and horizontal length scales of a planetary ra- 
dius, one thus estimates RoUH/L ~ lOmsec"^. Like- 



wise, for R = 3700Jkg-iK-\ Hp 



300 km, N 



0.005 sec~^ (appropriate for an isothermal layer on a hot 
Jupiter with surface gravity of 20msec~^), and Qnet /cp ~ 
10~^Ksec~^ (perhaps appropriate for photospheres of a 
typical hot Jupiter), Eq. (34) suggests mean vertical speeds 
of W ~ lOmscc^i (Showman and Guillot 2002; Show- 
man et al. 2008a). Thus, despite the fact that the circula- 
tion on hot Jupiters occurs in a radiative zone expected to 
be stably stratified over most of the planet, vertical over- 
turning can occur near the photosphere over timescales of 



Hp/W ~ 10^ sec. Because the interior of a multi-Gyr- 
old hot Jupiter transports an intrinsic flux of only ~ 10- 
100 Wm~^ (as compared with typically ~ 10^ Wm~^ in 
the layer where starlight is absorbed and infrared radia- 
tion escapes to space), the net heating gnet should decrease 
rapidly with depth. The mean vertical velocities should thus 
plummet with depth, leading to very long overturning times 
in the bottom part of the radiative zone. 

3.6. Effect of eddies: jet streams and banding 

Although turbulence is a challenging problem, several 
decades of work in the fields of fluid mechanics and at- 
mosphere/ocean dynamics have led to a basic understand- 
ing of turbulent flows and how they interact with plane- 
tary rotation. In a 3D turbulent flow, vorticity stretching 
and straining drives a fluid's kinetic energy to smaller and 
smaller length scales, where it can finally be removed by 
viscous dissipation. This process, the so-called /orwara? en- 
ergy cascade, occurs when turbulence has scales smaller 
than a fraction of a pressure scale height (~10 km for Earth, 
~200 km for hot Jupiters). On global scales, however, at- 
mospheres are quasi-two-dimensional, with horizontal flow 
dimensions exceeding vertical flow dimensions by typically 
factors of ~100. In this quasi- 2D regime, vortex stretch- 
ing is inhibited, and other nonlinear processes, such as vor- 
tex merging, assume a more prominent role, forcing energy 
to undergo an inverse cascade from small to large scales 
(VaUis 2006, pp. 349-361). This process has been well- 
documented in idealized laboratory and numerical exper- 
iments (e.g., Tabeling 2002) and helps explain the emer- 
gence of large-scale jets and vortices from small-scale tur- 
bulent forcing in planetary atmospheres. 

Where the Coriolis parameter / is approximately con- 
stant, this process generally produces turbulence that is hor- 
izontally isotropic, that is, turbulence without any preferred 
directionahty (north-south versus east-west). This could 
explain the existence of some quasi-circular vortices on 
Jupiter and Saturn, for example; however, it fails to explain 
the existence of numerous jet streams on Jupiter, Saturn, 
Uranus, Neptune, and the terrestrial planets. On the other 
hand, Rhines (1975) realized that the variation of / with lat- 
itude leads to anisotropy, causing elongation of structures 
in the east-west direction relative to the north-south direc- 
tion. This anisotropy can cause the energy to reorganize 
into east- west oriented jet streams with a characteristic lat- 
itudinal length scale 



1/2 



(35) 



^■^Oscillatory vertical motions associated with (for example) fast waves can 
occur adiabatically and are not governed by this constraint. 



where Uisa characteristic wind speed and /? = df/dy is the 
gradient of the Coriolis parameter with northward distance 
y. This length is called the Rhines scale. 

Numerous idealized studies of 2D turbulent flow forced 
by injection of small-scale turbulence have demonstrated 
jet formation by this mechanism [e.g., Williams (1978), 
Cho and Polvani (1996a), Cho and Polvani (1996b), Huang 
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Fig. 3. — Relative vorticity for three 2D non-divergent simulations on a sphere initialized from small-scale isotropic turbulence. The 
three simulations are identical except for the planetary rotation rate, which is zero on the left, intermediate in the middle, and fast on the 
right. Forcing and large-scale friction are zero, but the simulations contain a hyperviscosity to maintain numerical stability. The final 
state consists of isotropic turbulence in the non-rotating case but banded flow in the rotating case, a result of the Rhines effect. The 
spheres are viewed from the equatorial plane, with the rotation axis oriented upward on the page. After Yoden et al. (1999), Hayashi 
et al. (2000), Ishioka et al. (1999), and Hayashi et al. (2007). 



and Robinson (1998), Marcus et al. (2000); see Vasavada 
and Showman (2005) for a review]. Figure 3 illustrates 
an example. In 2D non-divergent numerical simulations 
of flow energized by small-scale turbulence, the flow re- 
mains isotropic at low rotation rates (Fig. 3, left) but devel- 
ops banded structure at high rotation rates (Fig. 3, right). 
As shown in Fig. 4, the number of jets in such simulations 
increases as the wind speed decreases, qualitatively consis- 
tent with Eq. (35). 

The Rhines scale can be interpreted as a transition scale 
between the regimes of turbulence and Rossby waves, 
which are a large-scale wave solution to the dynamical 
equations in the presence of non-zero (3. Considering the 
simplest possible case of an unforced 2D non-divergent 
fluid (§2.4), the vorticity equation reads 



5C 
dt 



■ V • VC + w/? = 0, 



(36) 



where ( is the relative vorticity, v is the horizontal veloc- 
ity vector, and v is the northward velocity component. The 
relative vorticity has characteristic scale U/L, and so the 
nonlinear term has characteristic scale U"^ /L'^. The /3 term 
has characteristic scale (3U . Comparison between these two 
terms shows that, for length scales smaller than the Rhines 
scale, the nonlinear term dominates, implying the domi- 
nance of nonlinear vorticity advection. For scales exceed- 
ing the Rhines scale, the linear v(3 term dominates over the 
nonlinear term, and Rossby waves are the primary solutions 
to the equations.'^ 



Adopting a streamfunction 4> defined by m = —dip/dy and v = dip/dx, 
the 2D non-divergent vorticity equation can be written d\7'^ip/dt + 
fSdi/j/dx = when the nonlinear term is neglected. In Cartesian ge- 
ometry with constant f3, this equation has wave solutions with a dispersion 
relationship u> = —I3k^/{k'^ + l'^) where uj is the oscillation frequency 
and k and I are the wavenumbers (2tt over the wavelength) in the eastward 
and northward directions, respectively. These are Rossby waves, which 



How does jet formation occur at the Rhines scale? In 
a 2D fluid where turbulence is injected at small scales'^ 
an upscale energy cascade occurs, but when the turbulent 
structures reach sizes comparable to the Rhines scale, the 
transition from nonlinear to quasi-linear dynamics prevents 
the turbulence from easily cascading to scales larger than 
Lp. At the Rhines scale, the characteristic Rossby-wave 
frequency for a typically oriented Rossby wave, ~f3/K, 
roughly matches the turbulent frequency Uk, where k is 
the wavenumber magnitude. (The Rossby wave frequency 
results from the dispersion relation, ignoring the distinction 
between eastward and northward wavenumbers; the turbu- 
lent frequency is one over an advective timescale, L/U.) At 
larger length scales (smaller wavenumbers) than the Rhines 
scale, the Rossby-wave oscillation frequency is faster than 
the turbulence frequency and so wave/turbulent interactions 
are inefficient at transporting energy into the Rossby-wave 
regime. This tends to cause a pile-up of energy at a wave- 
length L^. As a result, flow structures with a wavelength 
of Lp often contain more energy than flow structures at any 
other wavelength. 

But why are these dominant structures generally banded 
jet streams rather than (say) quasi-circular vortices? The 
reason relates to the anisotropy of Rossby waves. In a 
2D non-divergent fluid, the Rossby-wave dispersion rela- 
tion can be written 



(3k 



(3 cos a 



(37) 



where k and I are the wavenumbers in the eastward and 



have a westward phase speed and a frequency that depends on the wave 
orientation. 

^''In a 2D model, turbulence behaves in a 2D manner at all scales, but in 
the atmosphere, flows only tend to behave 2D when the ratio of widths to 
heights 2> 1 (implying horizontal scales exceeding hundreds of km for 
Earth). 
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northward directions, k = \/k'^ + P is the total wavenum- 
ber magnitude, and a. is the angle between the Rossby-wave 
propagation direction and east. Equating this frequency to 
the turbulence frequency U k leads to an anisotropic Rhines 
wavenumber 

fc^ = ^1 cosa|. (38) 

When plotted on the k-l wavenumber plane, this relation- 
ship traces out a dumbbell pattern; kfj approaches f3/U 
when a w or 180° and approaches zero when a « ±90°. 
The Rhines scale Lp is thus ^7r([///3)^/^ for wavevectors 
with a ~ or 180°, but it tends to infinity for wavevec- 
tors with a ~ ±90°. Because of this anisotropic barrier, 
the inverse cascade can successfully drive energy to smaller 
wavenumbers (larger length scales) along the a = ±90° 
axis than along any other direction. Wavevectors with 
a w ±90° correspond to cases with k I, that is, 
flow structures whose east-west lengths become unbounded 
compared to their north-south lengths. The result is usu- 
ally an east-west banded pattern and jet streams whose lat- 
itudinal width is ~7r(C///3)^/^. Although these consider- 
ations are heuristic, detailed numerical simulations gener- 
ally confirm that in 2D flows forced by small-scale isotropic 
turbulence, the tmbulence reorganizes to produce jets with 
widths close to (Fig. 3; for a review see Vasavada and 
Showman 2005). 

In practice, when small-scale forcing is injected into the 
flow, jet formation does not generally occur as the end prod- 
uct of successive mergers of continually larger and larger 
vortices but rather by a feedback whereby small eddies in- 
teract directly with the large-scale flow and pump up jets 
(e.g., Nozawa and Yoden 1997; Huang and Robinson 1998). 
The large-scale background shear associated with jets dis- 
torts the eddies, and the interaction pumps momentum up- 
gradient into the jet cores, helping to maintain the jets 
against friction or other forces. 

The above arguments were derived in the context of 
a 2D non-divergent model, which is the simplest model 
where Rossby waves can interact with turbulence to form 
jets. Such a model lacks gravity waves and has an infi- 
nite deformation radius Lb (Eqs. 10 and 23). It is possible 
to extend the above ideas to one-layer models exhibiting 
gravity waves and a finite deformation radius, such as the 
shallow- water model. Doing so suggests that, as long as 
P is strong, jets dominate when the deformation radius is 
large (as in the 2D non-divergent model); however, the flow 
becomes more isotropic (dominated by vortices rather than 
jets) when the deformation radius is sufficiently smaller 
than (;7//3) V2 (e.g., Okuno andMasuda 2003; Smitii 2004; 
Showman 2007). 

On a spherical planet, (3 = 2Qa~^ cos (p, which implies 
that a planet should have a number of jet streams roughly 
given by 

«,,~(^)'^ 

These considerations do a reasonable job of predicting 



the latitudinal separations and total number of jets that ex- 
ist on planets in our Solar System. Given the known wind 
speeds (Table 1) and values of /?, the Rhines scale predicts 
~ 20 jet streams on Jupiter/Satum, 4 jet streams on Uranus 
and Neptune, and 7 jet streams on Earth. This is similar 
to the observed numbers (~ 20 on Jupiter and Saturn, 3 
on Uranus and Neptune, and 3 to 7 on Earth depending 
on how they are defined'^). Moreover, the scale-dependent 
anisotropy — with quasi-isotropic turbulence at small scales 
and banded flow at large scales — ^is readily apparent on 
Solar-System planets. Fig. 5 illustrates an example for Sat- 
urn: small-scale structures (e.g., the cloud-covered vortices 
that manifest as small dark spots in the figure) are relatively 
circular, whereas the large-scale structures are banded. 

These considerations yield insight into the degree of 
banding and size of dominant flow structures to be ex- 
pected on exoplanets. Equation (39) suggests that the num- 
ber of bands scales as Q^/^, so rapidly rotating exoplan- 
ets will tend to exhibit numerous strongly banded jets, 
whereas slowly rotating planets will exhibit fewer jets with 
weaker banding. Moreover, for a given rotation rate, planets 
with slower winds should exhibit more bands than planets 
with faster winds. Hot Jupiters are expected to be tidally 
locked, implying rotation rates of typically a few days. 
Although wind speeds on hot Jupiters are unknown, vari- 
ous numerical models have obtained (or assumed) speeds 
of ~0.5-3 kmsec"^ (Showman and Guillot 2002; Cooper 
and Showman 2005; Langton and Laughlin 2007; Dobbs- 
Dixon and Lin 2008; Menou and Rauscher 2009; Cho et al. 
2003, 2008). Inserting these values into Eq. (39) implies 
A/jet ~ 1-2. Consistent with these arguments, simulations 
of hot-Jupiter atmospheres forced at small scales have ob- 
tained typicaUy ~l-3 broad jets (Cho et al. 2003, 2008; 
Langton and Laughlin 2008a). Thus, due to a combi- 
nation of their slower rotation and presumed faster wind 
speeds, hot Jupiters should have only a small number of 
broad jets — ^unUke Jupiter and Saturn. If the wind speeds 
are much slower, and/or rotation rates much greater than 
assumed here, the number of bands would be larger. 

3.7. Role of eddies in 3D atmospheres 

Although atmospheres behave in a quasi-2D manner on 
large scales due to large horizontakvertical aspect ratios, 

^^In Earth's troposphere, an instantaneous snapshot typically shows distinct 
subtropical and high-latitude jets in each hemisphere, with weaker flow be- 
tween these jets and westward flow at the equator. The jet latitudes exhibit 
large excursions in longitude and time, however, and when a time- and 
longitude average is performed, only a single local maximum in eastward 
flow exists in each hemisphere, with westward flow at the equator. 

^® Published 3D models of hot Jupiters forced by day-night heating contrasts 
also exhibit only a few broad jets (e.g.. Cooper and Showman 2005; Show- 
man et al. 2008a, 2009; Dobbs-Dixon and Lin 2008; Menou and Rauscher 
2009) but for a different reason. In this case, the jet width may be con- 
trolled by the Rossby deformation radius, which is close to the planetary 
radius for many hot Jupiters. As a result, the day-night heating contrast in- 
jects significant energy into the flow directly at the planetary scale, and the 
/3 effect can reorganize this global-scale energy into a banded jet pattern. 
Because the forcing and jet scales are comparable, such a flow could lack 
an inverse cascade. 
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Fig. 4. — Zonal-mean zonal winds versus latitude for a series of global, 2D non-divergent simulations performed on a sphere showing 
the correlation between jet speed and jet width. The simulations are forced by small-scale turbulence and damped by a linear drag. Each 
simulation uses different forcing and friction parameters and equilibrates to a different mean wind speed, ranging from fast on the left to 
slow on the right. Simulations with faster jets have fewer jets and vice versa; the jet widths approximately scale with the Rhines length. 
From Huang and Robinson (1998). 




Fig. 5. — Saturn's north polar region as imaged by the Cassini Visible and Infrared Mapping Spectrometer at 5 fim wavelengths. At 
this wavelength, scattered sunlight is negligible; bright regions are cloud-free regions where thermal emission escapes from the deep 
(~ 3 — 5 bar) atmosphere, whereas dark region are covered by thick clouds that block this thermal radiation. Note the scale-dependent 
anisotropy: small-scale features tend to be circular, whereas the large-scale features are banded. This phenomenon results directly from 
the Rhines effect (see text). [Photo credit: NASAA'IMS/Bob Brown/Kevin Baines.] 
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small Rossby numbers, and stable stratification, they are of 
course three dimensional, and as a result they can experi- 
ence both upscale and downscale energy cascades depend- 
ing on the length scales of the turbulence and other factors. 
Nevertheless, the basic mechanism discussed above for the 
interaction of turbulence and the (3 effect still appUes and 
suggests that even 3D atmospheres can generally exhibit a 
banded structure with a characteristic length scale close to 
Lfj. Consistent with this idea, numerical simulations show 
that banding can indeed occur in 3D even when /3 is the only 
source of horizontal anisotropy (e.g., Sayanagi et al. 2008; 
Lian and Showman 2009). Nevertheless, on real planets, 
banding can also result from anisotropic forcing — such as 
the fact that solar heating is primarily a function of latitude 
rather than longitude on rapidly rotating planets like Earth 
and Mars. 

A variety of studies show that, in stably stratified, rapidly 
rotating atmospheres, the characteristic vertical length scale 
of the flow is approximately f /N times the characteristic 
horizontal scale (e.g., Charney 1971; Dritschel et al. 1999; 
Reinaud et al. 2003; Haynes 2005). For typical horizontal 
dimensions of large-scale flows, this often implies vertical 
dimensions of one to several scale heights. 

Several processes can generate turbulent eddies that sig- 
nificantly affect the large-scale flow. Convection is particu- 
larly important in giant planet interiors and near the surface 
of terrestrial planets. Shear instabilities can occur when the 
wind shear is sufficiently large; they transfer energy from 
the mean flow into turbulence and reduce the shear of the 
mean flow. A particularly important turbulence-generating 
process on rapidly rotating planets is baroclinic instabil- 
ity, which is a dynamical instability driven by the extrac- 
tion of potential energy from a latitudinal temperature con- 
trast. On non- (or slowly) rotating planets, the presence of 
a latitudinal temperature contrast would simply cause a di- 
rect Hadley-type overturning circulation, which efficiently 
mutes the thermal contrasts. However, on a rapidly rotat- 
ing planet, Hadley circulations cannot penetrate to high lat- 
itudes (§5.3), and, in the absence of instabilities and waves, 
the atmosphere poleward of the Hadley cell would ap- 
proach a radiative-equilibrium temperature structure, with 
strong latitudinal temperature gradients and strong zonal 
winds in thermal-wind balance with the temperature gra- 
dients. In 3D, this radiative -equilibrium structure can ex- 
perience baroclinic instabilities, which develop into three- 
dimensional eddies that push cold polar air equatorward and 
down and push warm low-latitude air poleward and up. This 
process lowers the center of mass of the fluid, thereby con- 
verting potential energy into kinetic energy and transporting 
thermal energy between the equator and poles. The fastest 
instabihty growth rates occur at length scales comparable to 
the Rossby deformation radius (Eq. 23), and the resulting 



^^When two adjacent, stably stratified air columns have differing temper- 
ature profiles, potential energy is released when the colder column slides 
imdemeath the hotter column. For this reason, baroclinic instability, which 
draws on this energy source, is sometimes called "slantwise convection." 



eddies act as a major driver for the mid-latitude jet streams 
on Earth, Mars, and perhaps Jupiter, Satum, Uranus, and 
Neptune. 

A formalism for describing the effect of eddies on the 
mean flow can be achieved by decomposing the flow into 
zonal-mean components and deviations from the zonal 
mean. Here, "eddies" are defined as deviations from the 
zonal-mean flow and can represent the effects of turbulence, 
waves, and instabiUties. Denoting zonal means with over- 
bars and deviations therefrom with primes, we can write 
u = u + u',v = v + v',9 — 9 + 6', and lo = iD + lo' . Sub- 
stituting these expressions into the 3D primitive equations 
(Eq. 5) and zonally averaging the resulting equations yields 
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and similarly for zonal averages of other quadratic terms. 
Thus, the terms involving derivatives on the right sides of 
Eqs. (40)-(41) represent the effect of transport by advection 
of the mean flow (e.g., uv) and the effect of transport by ed- 
dies (e.g., u'v'). For example, u'v' can be thought of as the 
northward transport of eastward momentum by eddies. If 
u'v' has the same sign as the background jet shear du/dy, 
then eddies pump momentum up-gradient into the jet cores. 
If they have opposite signs, the eddies transport momen- 
tum down-gradient out of the jet cores. (See Cho 2008, 
for a brief discussion in the exoplanet context). Molecu- 
lar diffusion is a down-gradient process, but in many cases 
eddies can transport momentum up-gradient, strengthening 
the jets. 

To illustrate, consider linear Rossby waves, whose dis- 
persion relation is given by Eq. (37) in the case of a 2D 
non-divergent flow. For such waves, it can be shown that 
the product u'v' is opposite in sign to the north-south com- 
ponent of their group velocity (see Vallis 2006, pp. 489- 
490). Since the group velocity generally points away from 
the wave source, this implies that Rossby waves cause a flux 
of eastward momentum toward the wave source. The gen- 
eral result is an eastward eddy acceleration at the latitudes 
of the wave source and a westward eddy acceleration at the 
latitudes where the waves dissipate or break. This process is 
a major mechanism by which jet streams can form in plan- 
etary atmospheres. 

The existence of eddy acceleration and heat transport 
triggers a mean circulation that has a back-influence on 
both u and 6. On a rapidly rotating planet, for example, 
the eddy accelerations and heat transports tend to push the 
atmosphere away from geostrophic balance, leading to a 
mismatch in the mean north-south pressure gradient and 
Coriolis forces. This unbalanced force drives north-south 
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motion, which, through the mass continuity equation, ac- 
companies vertical motions. This circulation induces ad- 
vection of momentum and entropy, and moreover causes 
an east-west CorioUs acceleration fv, all of which tend to 
restore the atmosphere toward geostrophic balance [for de- 
scriptions of this, see James (1994, pp. 100-107) or Holton 
(2004, pp. 313-323)]. Such mean circulations tend to have 
broad horizontal and vertical extents, providing a mecha- 
nism for a localized eddy perturbation to trigger a non-local 
response. 

Clearly, jets on planets do not continually accelerate, so 
any jet acceleration caused by eddies must (on average) be 
resisted by other terms in the equation. Near the equator (or 
on slowly rotating planets where Coriolis forces are small), 
this jet acceleration is generally balanced by a combination 
of friction (e.g., as represented by V in Eq. (41)) and large- 
scale advection. However, in the mid- and high-latitudes 
of rapidly rotating planets, Coriolis forces are strong, and 
away from the surface the eddy accelerations are often bal- 
anced by the Coriolis force on the mean flow fv. 

When averaged vertically, the east-west accelerations 
due to the Coriolis forces tend to cancel out^", and the ver- 
tically integrated force balance (on a terrestrial planet, for 
example) is between the vertically integrated eddy accel- 
erations and surface friction (e.g., Peixoto and Oort 1992, 
Chapter 1 1). This means that regions without active eddy 
accelerations must have weak time-averaged surface winds, 
although the time-averaged winds may still be strong aloft. 
Conversely, the existence of strong time-averaged winds at 
the surface can only occur at latitudes where eddy accelera- 
tions are active. A prime example is the eastward mean sur- 
face winds in Earth's mid-latitudes (associated with the jet 
stream), which are enabled by the eddy accelerations from 
mid-latitude barocUnic instabilities. 

4. CIRCULATION REGIMES: GIANT PLANETS 

4.1. Thermal structure of giant planets: general con- 
siderations 

Because of the enormous energy released during plan- 
etary accretion, the interiors of giant planets are hot, and 
due to the expected high opacities in their interiors, this en- 
ergy is thought to be transported through their interiors into 
their atmospheres primarily by convection rather than radia- 
tion (Stevenson 1991; Chabrier and Baraffe 2000; Burrows 
et al. 2001; Guillot et al. 2004; Guillot 2005; Fortney et al. 
2009, this volume). Jupiter, Saturn, and Neptune emit sig- 
nificantly more energy than they absorb from the Sun (by 
factors of 1.7, 1.8, and 2.6, respectively, implying interior 
heat fluxes of ^1-lOWm^^) showing that these planets 
are still cooling off, presumably by convection, even 4.6 
Gyr after their formation. Even for giant exoplanets that 
lie close to their parent stars, the stellar insolation generally 



■^"This occurs because mass continuity requires that, in a time average, tliere 
be no net north-south mass transport, hence no net east-west Coriolis ac- 
celeration. 



cannot stop the inexorable cooling of the planetary interiors 
(e.g., Guillot et al. 1996; Saumon et al. 1996; Burrows et al. 
2000; Chabrier et al. 2004; Fortney et al. 2007). Thus, we 
expect that in general the interiors of giant planets will be 
convective, and hence the interior entropy will be nearly ho- 
mogenized and the interior temperature and density profiles 
will lie close to an adiabat.^' 

The interiors will not precisely follow an adiabat because 
convective heat loss generates descending plumes that are 
colder than the background fluid. We can estimate the devi- 
ation from an adiabat as follows. Solar-System giant planets 
are rotationally dominated; as a result, convective velocities 
are expected to scale as (Stevenson 1979; Fernando et al. 
1991): 

/ Fag \ 

\UpCp J 

where w is the characteristic magnitude of the vertical ve- 
locity, F is the heat flux transported by convection, and a 
and Cp are the thermal expansivity and specific heat at con- 
stant pressure. To order-of-magnitude, the convected heat 
flux should be ^ pwcp ST, where ST is the characteristic 
magnitude of the temperature difference between a convec- 
tive plume and the environment. These equations yield 
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Inserting values for Jupiter's interior (F 
Q = 1.74xl0-^sec-\/9~ 1000kgm-3,a ~ IQ-^K-^, 
Cp w 1.3 X 10''Jkg-iK-\ and g w 20mscc~2), we 
obtain ST ~ 10"'^ K, suggesting fractional density pertur- 
bations of aST ~ 10~^. Thus, deviations from an adi- 
abat in the interior should be extremely small. Even ob- 
jects with large interior heat fluxes, such as brown dwarfs 
or young giant planets, will have only modest deviations 
from an adiabat; for example, the above equations suggest 
that an isolated 1000-K Jupiter-like planet with a heat flux 
of 6 X 10*Wm~^ would experience deviations from an 
adiabat in its interior of only ^0.1 K. 

At sufficiently low pressures in the atmosphere, the 
opacities become small enough that the outward energy 
transport transitions from convection to radiation, leading to 
a temperature profile that in radiative equilibrium is stably 
stratified and hence suppresses convection. At a minimum, 
this transition (called the radiative-convective boundary) 
will occur when the gas becomes optically thin to escaping 
infrared radiation (at pressures less than ~0.01-1 bar de- 
pending on the opacities). For Jupiter, Saturn, Uranus, and 
Neptune, this transition occurs at pressures somewhat less 
than 1 bar. However, in the presence of intense stellar irradi- 
ation, the absorbed stellar energy greatly exceeds the energy 
loss from the interior; thus, the mean photospheric temper- 
ature only sUghtly exceeds the temperature that would exist 



^ Compositional gradients that force the interior density structures to deviate 
from a uniform adiabat could exist in Uranus and some giant exoplanets 
(Podolak et al. 1991; Chabrier and Baraffe 2007). 
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in thermal equilibrium with the star. In this case, cooling of 
the interior adiabat can only continue by development of a 
thick, stably stratified layer that penetrates downward from 
the surface and deepens with increasing age. Thus, for old, 
heavily irradiated planets, the radiative-convective bound- 
ary can instead occur at large optical depth, at pressures of 
^100-1000 bars depending on age (see §4.3 and Fortney 
et al. 2009, this volume). 

To what extent can large horizontal temperature differ- 
ences develop in the observable atmosphere? This de- 
pends greatly on the extent to which the infrared photo- 
sphere^^ and radiative-convective boundaries differ in pres- 
sure. On Jupiter, Saturn, Uranus, and Neptune, the in- 
frared photospheres lie at ^300 mbar, close to the pressure 
of the radiative-convective boundary. Thus, the convective 
interiors — with their exceedingly small lateral temperature 
contrasts — effectively outcrop into the layer where radia- 
tion streams into space. On such a planet, lateral tempera- 
ture contrasts in the observable atmosphere should be small 
(Ingersoll and Porco 1978), and indeed the latitudinal tem- 
perature contrasts are typically only ~2-5 K in the upper 
tropospheres and lower stratospheres of our Solar-System 
giant planets (Ingersoll 1990). 

If instead the radiative-convective boundary lies far 
deeper than the infrared photosphere — as occurs on hot 
Jupiters — then large lateral temperature contrasts at the 
photosphere can potentially develop. Above the radiative- 
convective boundary, convective mixing is inhibited and 
thus cannot force vertical air columns in different regions 
to he along a single temperature profile. Laterally vary- 
ing absorption of light from the parent star (or other pro- 
cesses) can thus generate different vertical temperature gra- 
dients in different regions, potentially leading to significant 
horizontal temperature contrasts. In reality, such lateral 
temperature contrasts lead to horizontal pressure gradients, 
which generate horizontal winds that attempt to reduce the 
temperature contrasts. The resulting temperature contrasts 
thus depend on a competition between radiation and dy- 
namics. These processes will be particularly important 
when strongly uneven external irradiation occurs, as on 
hot Jupiters. 

Several additional processes can produce meteorologi- 
cally significant temperature perturbations on giant plan- 
ets. In particular, most giant planets contain trace con- 
stituents that are gaseous in the high pressure/temperature 
conditions of the interiors but condense in the outermost 
layers. This condensation releases latent heat and changes 
the mean molecular weight of the air, both of which cause 
significant density perturbations. On Jupiter and Satum, 
ammonia, ammonium hydrosulfide, and water condense at 
pressures of ~0.5, 2, and 6 bars, causing temperature in- 
creases of ~0.2, 0.1, and 2K, respectively, for solar abun- 
dances. For Jupiter, this corresponds to fractional density 
changes of -IQ-^ for NH3 and NH4SH and -IQ-^ for 



'We define the infrared photosphere as the pressure at which the bulk of the 
planet's infrared radiation escapes to space. 



H2O. At deep levels (~2000K, which occurs at pressures 
of ~5000bars on Jupiter), iron, silicates, and various metal 
oxides and hydrides condense, with a total latent heating 
of ~1K and a fractional density perturbation of ~10~^. 
For Uranus/Neptune, the condensation structure is similar 
but also includes methane near the tropopause, with a la- 
tent heating of ~0.2 K and a fractional density change of 
~3 X lO^-^ at solar abundance. Conversion between ortho 
and para forms of the H2 molecule is also important at tem- 
peratures < 200 K and produces a fractional density change 
of ~10~^ (Gierasch and Conrath 1985). Note that methane, 
and plausibly water and other trace constituents, are en- 
hanced over solar abundances by factors of ~3, 7, and 20- 
40 for Jupiter, Saturn, and Uranus/Neptune, respectively, so 
the actual density alterations for these planets likely exceed 
those listed above by these factors. 

The key point is that these fractional density perturba- 
tions exceed those associated with dry convection in the 
interior by orders of magnitude; condensation could thus 
dominate the meteorology in the region where it occurs. 
Such moist convection could have several effects. First, 
if the density contrasts become organized on large scales, 
then significant vertical shear of the horizontal wind can 
occur (via the thermal-wind equation) that would otherwise 
not exist in the outermost part of the convection zone. Sec- 
ond, it can generate a background temperature profile that 
is stably stratified (e.g. Stevens 2005), allowing wave prop- 
agation and various other phenomena. Third, the buoy- 
ancy produced by the condensation could act as a powerful 
driver of the circulation. For example, on Jupiter and Sat- 
um, condensation leads to powerful thunderstorms (Little 
et al. 1999; Gierasch et al. 2000) that are a leading can- 
didate for driving the global-scale jet streams (§4.2). The 
cloud formation associated with moist convection can also 
significantly alter the temperature structure due to its influ- 
ence in scattering/absorbing radiation. 

For giant planets hotter than those in our Solar System, 
the condensation sequence shifts to lower pressures and the 
more volatile species disappear from the condensation se- 
quence. Important breakpoints occur for giant exoplanets 
with effective temperatures exceeding ^150 and -^300- 
500 K, above which ammonia and water, respectively, no 
longer condense, thus removing the effects of their conden- 
sation from the meteorology. The latter breakpoint will be 
particularly significant due to the overriding dominance of 
the fractional density change associated with water conden- 
sation. Nevertheless, condensation of iron, silicates, and 
metal oxides and hydrides will collectively constitute an 
important source of buoyancy production even for hot giant 
exoplanets, as long as such condensation occurs within the 
convection zone, where its buoyancy dominates over that 
associated with dry convection (Eq. 44). For hot Jupiters, 
however, evolution models predict the existence of a quasi- 
isothermal radiative zone extending down to pressures of 
~100-1000 bars (GuiUot 2005; Fortney et al. 2009). Latent 
heating occurring within this zone is less likely to be impor- 
tant, simply because the ~1 K potential temperature change 
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associated with silicate condensation is small compared to 
the probable vertical and horizontal potential temperature 
variations associated with the global circulation (§4.3). For 
example, in the vertically isothermal radiative zone of a 
hot Jupiter, a vertical displacement of air by only ~0.5km 
would produce a temperature perturbation (relative to the 
background) of ~1 K, similar to the magnitude of thermal 
variation caused by condensation. 

4.2. Circulation of giant planets: general considera- 
tions 

We now turn to the global-scale atmospheric circulation. 
As yet, few observational constraints exist regarding the cir- 
culation of exoplanets. Infrared light curves of several hot 
Jupiters indirectly suggest the presence of fast winds able to 
advect the temperature pattem (Knutson et al. 2007; Cowan 
et al. 2007; Knutson et al. 2009b), and searches for variabil- 
ity in brown dwarfs are being made (e.g.. Morales -Calderon 
et al. 2006; Goldman et al. 2008). But, at present, our lo- 
cal Solar-System giant planets represent the best proxies to 
guide our thinking about the circulation on a wide class of 
rapidly rotating giant exoplanets. We consider basic issues 
here and take up specific models of hot Jupiters in §4.3. 

In the observable atmosphere, the global-scale circula- 
tion on Jupiter, Saturn, Uranus, and Neptune is dominated 
by numerous east-west jet streams (Fig. 6). Jupiter and 
Saturn each exhibit ~ 20 jet streams whereas Uranus and 
Neptune exhibit three jets each. Peak speeds range from 
~150mscc~^ for Jupiter to over 400mscc^^ for Saturn 
and Neptune. The winds are determined by tracking the 
motion of small cloud features, which occur at pressures of 
~0.5-4 bars, over periods of hours. The jet streams modu- 
late the patterns of ascent and descent, leading to a banded 
cloud pattem that can be seen in Fig. 6. 

Constraints from balance arguments: For giant exoplan- 
ets and Solar-System giant planets alike, the dynamical 
link between temperatures and winds encapsulated by the 
thermal-wind equation allows plausible scenarios for the 
vertical structure of the winds to be identified. The Uke- 
lihood that giant planets lose their internal energy by con- 
vection suggests that the interior entropy is nearly homog- 
enized and hence that the interiors are close to a barotropic 
state (§3.3). This would lead to the Taylor-Proudman theo- 
rem (Eqs. 25-27), which states that the winds are constant 
on cyUnders parallel to the planetary rotation axis. This 
scenario, first suggested by Busse (1976) for Jupiter, Sat- 
urn, Uranus, and Neptune, postulates that the jets observed 
in the cloud layer would simply represent the intersection 
with the surface of eastward- and westward-moving Taylor 
columns that penetrate throughout the molecular envelope. 

How closely will the molecular interior of a giant planet 
follow the Taylor-Proudman theorem? It might not apply if 
geostrophy does not hold on the length scale of the jets (i.e. 
Rossby number > 1), but this is unUkely, at least for the 
rapidly rotating giant planets of our Solar System (Liu et al. 
2008). However, for giant exoplanets, geostrophy could 



break down if the interior wind speeds are sufficiently fast, 
planetary rotation rate is sufficiently low, the interior heat 
flux is sufficiently large (so that convective buoyancy forces 
exceed Coriolis forces Q.u associated with the jets, for ex- 
ample), or magnetic effects dominate. Altemately, even if 
the interior exhibits geostrophic columnar behavior, the in- 
terior could exhibit a shear of the wind along the Taylor 
columns if density variations on isobars in the interior are 
sufficiently large. To illustrate how this might work, we 
cast Eq. (24) in a cylindrical coordinate system and take the 
longitudinal component, which gives 
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where u is the zonal wind, is the coordinate parallel to 
the rotation axis (see §3.3), and A is the unit vector in the 
longitudinal direction. Suppose that horizontal density vari- 
ations (on isobars) occur only in latitude. A geometrical 
argument shows that the magnitude of Vp x Vp is then 
\Vp\{dp/dy)p, where {dp/dy)p is the partial derivative of 
density with northward distance, y, on a surface of constant 
pressure. The magnitude of the pressure gradient is domi- 
nated by the hydrostatic component, —pg, and thus we can 
write 

du g [ dp^ 
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which, to order-of-magnitude, implies that the magnitude 
of zonal-wind variation along a Taylor column can be ex- 
pressed as 

where L is the width of the jets in latitude, {dp/p)p is the 
fractional density difference (on isobars) across a jet, and 
Az* is the distance along the Taylor column over which 
this fractional density difference occurs. In the interior, 
convection suppUes the fractional density contrasts, and as 
shown in §4. 1 these plausibly have extremely small values 
of 10~^-10~* below the region where condensation can oc- 
cur. We can obtain a crude order-of-magnitude estimate 
for the thermal-wind shear along a Taylor column by as- 
suming that these fractional density differences are orga- 
nized on the jet scale and coherently extend across the full 
^IC* km vertical length of a Taylor column. Inserting these 
values into Eq. (47), along with parameters relevant to a 
Jupiter-like planet (g « 20msec~^, L « 5000 km, Q = 
1.74 X IQ-^scc-i) yields Am - 0.001-0.1 mscc-i. If 
this estimate is correct, then Jupiter-like planets would ex- 
hibit only small thermal- wind shear along Taylor columns, 
and the Taylor-Proudman theorem would hold to good ap- 
proximation in the molecular interior. 

However, the above estimate is crude; if the jet-scale 
horizontal density contrasts were much larger than the den- 
sity contrasts associated with convective plumes, for exam- 
ple, then the wind shear along Taylor colunms could ex- 
ceed the values estimated above. The quantitative validity 
of Eq. (44), on which the wind-shear estimates are based. 
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Fig. 6. — Top: Jupiter, Saturn, Uranus, and Neptune shown to scale in visible-wavelength images from the Cassini and Voyager 
spacecraft. Bottom: Zonal-mean zonal wind profiles obtained from tracking small-scale cloud features over a period of hours. Jupiter 
and Saturn have an alternating pattern of ~ 20 east- west jets, whereas Uranus and Neptune have 3 broad jets. 



also remains uncertain. To attack the problem more rig- 
orously, Kaspi et al. (2009) performed fuU 3D numerical 
simulations of convection in giant-planet interiors, suggest- 
ing that the compressibility may play an important role in 
allowing the generation of thermal-wind shear within giant 
planets, especially in the outermost layers of their convec- 
tion zones. Such wind shear is particularly likely to be im- 
portant for young and/or massive planets with large interior 
heat fluxes. 

Although fast winds could exist throughout the molec- 
ular envelope, Lorentz forces probably act to brake the 
zonal flows in the underlying metallic region at pressures 
exceeding ^l-3Mbar (Kirk and Stevenson 1987; Grote 
et al. 2000; Busse 2002). The winds in the deep, metal- 
lic interior are thus often assumed to be weak. The transi- 
tion between molecular and metallic occurs gradually, and 
there exists a wide semi-conducting transition zone over 
which the Taylor-Proudman theorem approximately holds 
yet the electrical conductivity becomes important. Liu et al. 
(2008) argued that, if the observed jets on Jupiter and Sat- 
urn penetrated downward as Taylor columns, the Ohmic 
dissipation would exceed the luminosity of Jupiter (Sat- 
urn) if the jets extended deeper than 95% (85%) of the 
planetary radius. The observed jets are thus probably shal- 
lower than these depths. However, if the jets extended part- 
way through the molecular envelope, terminating within the 
semi-conducting zone where the Ohmic dissipation occurs, 
the shear at the base of the jets must coexist with lateral 
density contrasts on isobars via the thermal wind relation- 



ship (Eq. 46). This is problematic, because, as previously 
discussed, sufficiently large sources of density contrasts are 
lacking in the deep interior Liu et al. (2008) therefore 
argued that the jets must be weak throughout the molecu- 
lar envelope up to shallow levels where alternate buoyancy 
sources become available (e.g., latent heating and/or transi- 
tion to a radiative zone). 

In the outermost layers of a giant planet, density varia- 
tions associated with latent heating and/or a transition to a 
radiative zone allow significant thermal wind-shear to de- 
velop. On Jupiter and Saturn, this so-called "baroclinic" 
layer begins with condensation of iron, silicates, and vari- 
ous metal oxides and hydrides at pressures of ^^10"' bars and 
continues with condensation of water, NH4SH, and ammo- 
nia (at '-^lO, 2, and 0.5 bars, respectively), finally transition- 
ing to the stably stratified stratosphere at ~0.2 bars. Early 
models showed that the observed jet speeds can plausibly 
result from thermal-wind shear within this layer assuming 
the winds in the interior are zero (Ingersoll and Cuzzi 1969). 
For example, water condensation could cause a fractional 
density contrast on isobars of up to 10^^, and if these den- 
sity differences extend vertically over ^100 km (roughly 
the altitude difference between the water-condensation level 
and the observed cloud deck), then Eq. (47) implies Au ^ 
30 msec~^, similar to the observed speed of Jupiter's mid- 
latitude jets (Fig. 6). 

To summarize, current understanding suggests that giant 
planets should exhibit winds within the convective molecu- 
lar envelope that exhibit columnar structure parallel to the 
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rotation axis, ttansitioning in the outennost layers (pres- 
sures less than thousands of bars) to a baroclinic zone where 
horizontal density contrasts associated with latent heating 
and/or a radiative zone allow the zonal winds to vary with 
altitude. 

Constraints from observations: Few observations yet ex- 
ist that constrain the deep wind structure in Jupiter, Saturn, 
Uranus, and Neptune. On Neptune, gravity data from the 
1989 Voyager 2 flyby imply that the fast jets observed in 
the cloud layers are confined to the outermost few percent 
or less of the planet's mass, corresponding to pressures less 
than ~10^ bars (Hubbard et al. 1991). By comparison, only 
weak constraints currently exist for Jupiter and Saturn. The 
Galileo Probe, which entered Jupiter's atmosphere at a lat- 
itude of ~7°N in 1995, showed that the equatorial jet pen- 
etrates to at least 22 bars, roughly 150 km below the vis- 
ible cloud deck (Atkinson et al. 1997), and indirect infer- 
ences suggest that the jets penetrate to at least ~5-10 bars at 
other latitudes (e.g. Dowling and IngersoII 1989; Sanchez- 
Lavega et al. 2008). In 2016, however, NASA's Juno mis- 
sion will measure Jupiter's gravity field, finally determin- 
ing whether Jupiter's jets penetrate deeply or are confined 
to pressures as shallow as thousands of bars or less. These 
results will have important implications for understanding 
the deep- wind structure in giant exoplanets generally. 

Jet-pumping mechanisms: Two scenarios exist for the 
mechanisms to pump the zonal jets on the giant planets 
(for a review see Vasavada and Showman 2005). In the 
"shallow forcing" scenario, the jets are hypothesized to 
result from moist convection (e.g., thunderstorms), baro- 
clinic instabilities, or other turbulence-generating processes 
in the baroclinic layer in the outermost region of the planet. 
Two-dimensional and shallow-water models of this process, 
where random turbulence is injected as an imposed forc- 
ing (or alternately added as a turbulent initial condition), 
show success in producing altemating zonal jets through the 
Rhines effect (§3.6). For appropriately chosen forcing and 
damping parameters (or for appropriate initial velocities 
when the turbulence is added as an initial condition), such 
models typically exhibit jets of approximately the observed 
speed and spacing (WilUams 1978; Cho and Polvani 1996a; 
Huang and Robinson 1998; Scott and Polvani 2007; Show- 
man 2007; Sukoriansky et al. 2007). Most models of this 
type produce an equatorial jet that flows westward rather 
than eastward as on Jupiter and Saturn, though robust east- 
ward equatorial flow has been obtained in a recent shallow- 
water study (Scott and Polvani 2008). Three-dimensional 
models of this process can explicitly resolve the turbulent 
energy generation and therefore need not inject turbulence 
by hand. These models likewise exhibit multiple zonal jets 
via the Rhines effect, in some cases spontaneously pro- 
ducing an eastward equatorial jet as on Jupiter and Saturn 
(e.g. WilUams 1979, 2003; Lian and Showman 2008, 2009; 
Schneider and Liu 2009). Figure 7, for example, illustrates 
3D simulations from Lian and Showman (2009) where the 



circulation is driven by latent heating associated with con- 
densation of water vapor; the condensation produces tur- 
bulent eddies that interact with the planetary rotation (i.e., 
non-zero /?) to generate jets. Globes show the zonal wind 
for a Jupiter case (top row), Saturn case (middle), and a case 
representing Uranus/Neptune (bottom row). The Jupiter 
and Saturn cases develop ~ 20 jets, including an eastward 
equatorial jet, whereas the Uranus/Neptune case develops 
a 3-jet structure with a broad westward equatorial flow — 
qualitatively similar to the observed jet patterns on these 
planets (Fig. 6). Note that shallow /orc/ng need not imply 
shallow jets; because the atmosphere's response to a pertur- 
bation is nonlocal, deep jets (penetrating many scale heights 
below the clouds) can result from jet pumping confined to 
the cloud layer if the frictional damping in the interior is 
sufficiently weak (Showman et al. 2006; Lian and Show- 
man 2008). 

In the "deep forcing" scenario, convection throughout 
the molecular envelope is hypothesized to drive the jets. 
To date, most studies of this process involve 3D numeri- 
cal simulations of convection in a self -gravitating spherical 
shell excluding the effects of magnetohydrodynamics (e.g., 
Aurnou and Olson 2001; Christensen 2001, 2002; Heimpel 
et al. 2005). The boundaries are generally taken as free-slip 
impermeable spherical surfaces with an outer boundary at 
the planetary surface and an inner boundary at a radius of 
0.5 to 0.9 planetary radii. Most of these studies assume the 
mean density, thermal expansitivity, and other fluid proper- 
ties are constant with radius. These studies generally pro- 
duce jets penetrating through the shell as Taylor columns. 
When the shell is thick (inner radius < 0.7 of the outer 
radius), only ^ 3-5 jets form (Aurnou and Olson 2001; 
Christensen 2001, 2002; Aurnou and Heimpel 2004), in- 
consistent with Jupiter and Saturn. Only when the shell has 
a thickness ^ 10% or less of the planetary radius do such 
simulations produce a Jupiter- or Satum-hke profile with 
~ 20 jets (Heimpel et al. 2005). However, on giant planets, 
the density and thermal expansivity each vary by several or- 
ders of magnitude from the photosphere to the deep interior, 
and this might have a major effect on the dynamics. Only 
recently has this effect been included in detailed numeri- 
cal models of the convection (Glatzmaier et al. 2009; Kaspi 
et al. 2009), and these studies show that the compressibil- 
ity can exert a significant influence on the jet structure. A 
difficulty in interpreting all the convective simulations de- 
scribed here is that, to integrate in reasonable times, they 
must adopt viscosities and heat fluxes orders of magnitude 
too large. 

The fast interior jets in these convection simulations im- 
ply the existence of strong horizontal pressure contrasts, 
which in these models are supported by the impermeable 
upper and lower boundaries. On a giant planet, however, 
there is no solid surface to support such pressure variations. 
If these fast winds transitioned to weak winds within the 
metalUc region, significant lateral density variations would 
be required via the thermal-wind relationship. As discussed 
previously, however, there exists no obvious source of such 
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Fig. 7. — Illustration of the effect of latent heating on the circulation of a rapidly rotating giant planet. Shows zonal (east-west) wind 
from 3D atmospheric simulations of Jupiter (top row), Saturn (middle row), and a case representing Uranus/Neptune (bottom row) where 
the circulation is driven by condensation of water vapor, with an assumed deep abundance of 3, 5, and 30 times solar for the Jupiter, 
Saturn, and Uranus/Neptune cases, respectively. Left column shows an oblique view and right column shows a view looking down 
over the north pole. Note the development of numerous jets (including superrotating equatorial jets on Jupiter and Saturn), which bear 
qualitative resemblance to the observed jets in Fig. 6. From Lian and Showman (2009). 
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density variations (on isobars) within the deep interior. Re- 
solving this issue will require 3D convection simulations 
that incorporate magnetohydrodynamics and simulate the 
entire molecular+metaUic interior. Several groups are cur- 
rently making such efforts, so the next decade should see 
significant advances in this area. 

4.3. Hot Jupiters and Neptunes 

Because of their likelihood of transiting their stars, hot 
Jupiters (giant planets within 0. 1 AU of their primary star) 
remain the best characterized exoplanets and thus far have 
been the focus of most work on exoplanet atmospheric 
circulation. Dayside photometry and/or infrared spectra 
now exist for a variety of hot Jupiters, constraining the 
dayside temperature structure. For some planets, such as 
HD 189733b, this suggests a temperature profile that de- 
creases with altitude from ~ 0.01-1 bar (Charbonneau et al. 
2008; Barman 2008; Swain et al. 2009), whereas other hot 
Jupiters, such as HD 209458b, TrES-2, TrES-4, and XO-lb 
(Knutson et al. 2008, 2009a; Machalek et al. 2008), appear 
to exhibit a thermal inversion layer (i.e., a hot stratosphere) 
where temperatures rise above 2000 K. 

Infrared light curves at 8 and/or 24 ^m now exist for HD 
189733b, Ups And b, HD 209458b, and several other hot 
Jupiters (Knutson et al. 2007, 2009b; Harrington et al. 2006; 
Cowan et al. 2007), placing constraints on the day-night 
temperature distribution of these planets. Light curves for 
HD 189733b and HD 209458b exhibit nighlside brightness 
temperatures only modestly (~ 20-30%) cooler than the 
dayside brightness temperatures, suggesting efficient redis- 
tribution of the thermal energy from dayside to nightside. 
On the other hand, Ups And b and HD 179949b exhibit 
larger day-night phase variations that suggests large (per- 
haps > 500 K) day-night temperature differences. 

These and other observations provide sufficient con- 
straints on the dayside temperature structure and day-night 
temperature distributions to make comparison with detailed 
atmospheric circulation models a useful exercise. Here we 
survey the basic dynamical regime and recent efforts to 
model these objects (see also Showman et al. 2008b; Cho 
2008). 

The dynamical regime on hot Jupiters differs from that 
on Jupiter and Saturn in several important ways. First, 
because of the short spindown times due to their proxim- 
ity to their stars, hot Jupiters are expected to rotate nearly 
synchronously^^ with their orbital periods, implying rota- 
tion periods of 1-5 Earth days for most hot Jupiters dis- 
covered to date. This is 2-12 times slower than Jupiter's 
10-hour rotation period, implying that, compared to Jupiter, 
hot Jupiters experience significantly weaker Corolis forces 
for a given wind speed. Nevertheless, Coriohs forces can 
still play a key role: for global-scale flows (length scales 
L ^ 10^ m) and wind speeds of 2kinscc~-'^ (see below), 
the Rossby number for a hot Jupiter with a 3-day period is 



or pseudo-synchronously in the case of hot Jupiters on highly eccentric 
orbits 



~ 1, implying a three-way force balance between Coriolis, 
pressure-gradient, and inertial (i.e. advective) terms in the 
horizontal equation of motion. Slower winds would imply 
smaller Rossby numbers, implying greater rotational domi- 
nance. 

Second, hot Jupiters receive enormous energy fluxes 
(~ 10^-10^ W m~^ on a global average) from their parent 
stars, in contrast to the ~ 10Wm~^ received by Jupiter 
and 1 Wm~^ for Neptune. The resulting high tempera- 
tures lead to short radiative time constants of days or less 
at pressures < 1 bar (Showman and Guillot 2002; Iro et al. 
2005; Showman et al. 2008a, and Eq. (19)), in contrast to 
Jupiter where radiative time constants are years. Even in the 
presence of fast winds, these short radiative time constants 
allow the possibility of large fractional day-night tempera- 
ture differences, particularly on the most heavily irradiated 
planets — in contrast to Jupiter where temperatures vary hor- 
izontally by only a few percent. 

Third, over several Gyr of evolution, the intense irradi- 
ation received by hot Jupiters leads to the development of 
a nearly isothermal radiative zone extending to pressures of 
~ 100-1000 bars (GuiUot et al. 1996; GuiUot 2005; Fort- 
ney et al. 2009). The observable weather in hot Jupiters thus 
occurs not within the convection zone but within the stably 
stratified radiative zone. Moreover, as described in §4.1, 
this separation between the photosphere and the radiative- 
convective boundary may allow large horizontal tempera- 
ture differences to develop and support significant thermal- 
wind shear For example, assuming that a lateral tempera- 
ture contrast of AT ~ 300 K extending over several scale 
heights occurs over a lateral distance of 10^ m, the result- 
ing thermal-wind shear is Aw, ^ 2 kniscc^^. Thus, balance 
arguments suggest that the existence of significant horizon- 
tal temperature contrasts would naturally require fast wind 
speeds. 

The modest rotation rates, large stable stratifications, and 
possible fast wind speeds suggest that the Rossby defor- 
mation radius and Rhines scale are large on hot Jupiters. 
On Jupiter, the deformation radius is ~ 2000 km and the 
Rhines scale is ~ 10^ km (much smaller than the planetary 
radius of 71, 400 km), which helps explain why Jupiter and 
Saturn have an abundant population of small vortices and 
numerous zonal jets (Figs. 5-6) (WilUams 1978; Cho and 
Polvani 1996a; Vasavada and Showman 2005). In contrast, 
a hot Jupiter with a rotation period of 3 Earth days and wind 
speed of 2 km sec ~^ has a deformation radius and Rhines 
scale comparable to the planetary radius. Thus, dynamical 
structures such as vortices and jets should be more global in 
scale than occurs on Jupiter and Saturn (Menou et al. 2003; 
Showman and Guillot 2002; Cho et al. 2003). This scal- 
ing argument is supported by detailed nonlinear numerical 
simulations of the circulation, which generally show the de- 
velopment of only ^ 1-3 broad jets (Showman and Guillot 
2002; Cho et al. 2003, 2008; Dobbs-Dixon and Lin 2008; 
Langton and Laughlin 2007, 2008a; Showman et al. 2008a, 
2009). 

A variety of 2D and 3D models have been used to inves- 
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tigate the atmospheric circulation of hot Jupiters. Cho et al. 
(2003, 2008) performed global numerical simulations of hot 
Jupiters on circular orbits using the one-layer equivalent 
barotropic equations (a one-layer model that is mathemat- 
ically similar to the shallow-water equations described in 
§ 2.3). They initialized the simulations with small-scale bal- 
anced turbulence and forced them with a large-scale deflec- 
tion of the surface to provide a crude representation of the 
pressure effects of the day-night heating gradient. However, 
no expUcit heating/cooling was included. Together, these 
effects led to the production of several broad, meandering 
jets and drifting polar vortices (Fig. 8). The mean wind 
speed in the final state is to a large degree determined by 
the mean speed of the initial turbulence, which ranged from 
100-800 m sec~^ in their models. At the large (~planetary- 
scale) deformation radius relevant to hot Jupiters, the equa- 
torial jet in the final state can flow either eastward or west- 
ward depending on the initial condition and other details. 
The simulations exhibit significant time variabiUty that, if 
present on hot Jupiters, would lead to detectable orbit-to- 
orbit variabihty in light curves and secondary eclipse depths 
(Cho et al. 2003; Rauscher et al. 2007, 2008). 

Langton and Laughlin (2007) performed global, 2D sim- 
ulations of hot Jupiters on circular orbits using the shallow- 
water equations with a mass source on the dayside and 
mass sink on the nightside to parameterize the effects of 
dayside heating and nightside cooling. When the obhq- 
uity is assumed zero and the mass sources/sinks are suffi- 
ciently large, their forced flows quickly reach a steady state 
with wind speeds of ~ 1 km sec" ^ and order-unity lateral 
variations in the thickness of the shallow-water layer. On 
the other hand, Langton and Laughlin (2008a) numerically 
solved the 2D fully compressible equations for the horizon- 
tal velocity and temperature of hot Jupiters on eccentric or- 
bits and obtained very turbulent, time-varying flows. Ap- 
parently, in this case, the large-scale heating patterns pro- 
duced hemisphere-scale vortices that were dynamically un- 
stable, leading to the breakdown of these eddies into smaU- 
scale turbulence. 

Several authors have also performed 3D numerical sim- 
ulations of hot Jupiters. Showman and Guillot (2002), 
Cooper and Showman (2005, 2006), Showman et al. 
(2008a), and Menou and Rauscher (2009) performed global 
simulations with the 3D primitive equations where the day- 
side heating and nightside coohng was parameterized with a 
Newtonian heating/cooling scheme, which relaxes the tem- 
peratures toward a prescribed radiative-equilibrium tem- 
perature profile (hot on the dayside, cold on the nightside) 
over a prescribed radiative timescale. Dobbs-Dixon and Lin 
(2008) performed simulations with the fully compressible 
equations in a limited-area domain, consisting of the equa- 
torial and mid-latitudes but with the poles cut off. Like the 
studies listed above, they also adopted a simplified method 
for forcing their flow, in this case using a radiative diffusion 
scheme. Showman et al. (2009) coupled their global 3D 
dynamical solver to a state-of-the-art, non-gray, cloud-free 
radiative transfer scheme with opacities calculated assum- 



ing local chemical equilibrium (Figs. 9 and 10). The above 
models models all generally obtain wind structures with 
~ 1-3 broad jets with speeds of ^ l-4kmsec~^. 

Despite the diversity in modehng approaches, the studies 
described above agree in several key areas. 

• First, most of the above studies generally produce 
peak wind speeds similar to within a factor of 2- 
3, in the range of one to several kmsec^^.^"* This 
similarity is not a coincidence but results from the 
force balances that occur for a global-scale flow in 
the presence of large fractional temperature differ- 
ences. Consider the longitudinal force balance at the 
equator, for example, and suppose the day-night heat- 
ing gradient produces a day-night temperature dif- 
ference AThoriz that extends vertically over a range 
of log-pressures A In p. This temperature difference 
causes a day-night horizontal pressure-gradient ac- 
celeration which, to order-of-magnitude, can be writ- 
ten i?ATiiorizA Inp/a, where a is the planetary ra- 
dius. At high latitudes, this could be balanced by the 
Coriolis force arising from a north-south flow, but the 
horizontal Coriolis force is zero at the equator. At the 
equator, such a force instead tends to cause accelera- 
tion of the flow in the east-west direction. Balancing 
the pressure-gradient acceleration by v • Vv, which to 
order-of-magnitude is U'^/a for a global-scale flow, 
we have 

U - V^AThorizAlnp (48) 

This should be interpreted as the characteristic vari- 
ation in zonal wind speed along the equator. For 
R ^ 3700Jkg-iK-\ AThoriz ~ 400 K, and 
Alnp ~ 3 (appropriate to a temperature difference 
extending verticaUy over three scale heights), this 
yields U 2 km sec" ^. 

Likewise, consider the latitudinal force balance in 
the mid-latitudes. To order-of-magnitude, the lati- 
tudinal pressure-gradient acceleration is again given 
by i^AThorizA Inp/a, where here AThoriz is the lat- 
itudinal temperature contrast (e.g., from equator to 
pole) that extends vertically over A In p. If Ro^ 1, 
this is balanced by the advective acceleration U"^ /a, 
whereas if i?o <C 1, it would instead be balanced by 
the Coriolis acceleration fU, where / is the Coriohs 
parameter (§3.3). The former case recovers Eq. (48), 
whereas the latter case yields 



U 



i?AThorizAlnp 
fa 



(49) 



Here, U is properly interpreted as the characteristic 
difference in horizontal wind speed vertically across 
Aln^^. Inserting R = 3700 J kg"! K-\ AThoriz ~ 

■^''in intercomparing studies, one must be careful to distinguish mean versus 
peak speeds and, in the case of 3D models, the pressure level at which 
those speeds are quoted; such quantities can differ by a factor of several in 
a single model. 
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Fig. 8. — One-layer, global equivalent barotropic simulation of the hot Jupiter HD 209458b from Cho et al. (2003). Globes show 
equatorial (left) and polar (right) views of the potential vorticity (see §2); note the development of large polar vortices and small-scale 
turbulent structure, which result from a combination of the turbulent initial condition and large-scale "topographic" forcing intended to 
qualitatively represent the day-night heating contrast. Plot shows zonal-mean zonal wind versus latitude, illustrating the 3-jet structure 
that develops. 




Fig. 9. — Results from a 3D simulation of the hot Jupiter HD 189733b from Showman et al. (2009). The dynamics are coupled to 
a realistic representation of cloud-free, non-gray radiative transfer assuming solar abundances. Left: Temperature (colorscale, in K) 
and winds (arrows) at the 30 mbar pressure, which is the approximate level of the mid-IR photosphere. A strong eastward equatorial 
jet develops that displaces the hottest regions eastward from the substellar point (which lies at 0° longitude, 0° latitude). Right: Light 
curves in Spitzer bandpasses calculated from the simulation (curves; labels show wavelength in /im) in comparison to observations 
(points) from Knutson et al. (2007, 2009b), Charbonneau et al. (2008), and Deming et al. (2006). 
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Fig. 10. — 3D simulation of HD 209458b, including realistic cloud-free, non-gray radiative transfer, illustrating development of a 
dayside stratosphere when opacity from visible-absorbing species (in this case TiO and VO) are included. Such stratospheres are 
relevant to explaining Spitzer data of HD 209458b, TrES-2, TrES-4, XO-lb, and other hot Jupiters. Here, the stratosphere (shown in 
red) begins at pressures of ~ 10 mbar and widens with altitude until it covers most of the dayside at pressures of ~ 0.1 mbar. Panels 
show temperature (colorscale, in K) and winds (arrows) at the 0.08, 1, and 30 mbar levels from top to bottom, respectively. Substellar 
longitude and latitude are (0°, 0°). From Showman et al. (2009). 
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400K, Alnp ~ 3, / ~ 2 X IQ-^ (appropriate in 
mid-latitudes to a hot Jupiter with a 3-day period), 
and aw 10^ m, we again obtain U ^ 2 kniscc^^. 

While not minimizing the real differences in the nu- 
merical results obtained in the various studies, the 
above estimates show that the existence of wind 
speeds of a few km sec ~^ in the various numerical 
studies is a basic outcome of force balance in the 
presence of lateral temperature contrasts of hundreds 
ofK. 

• Second, the various numerical studies all produce a 
small number (~ 1^) of broad jets. This was first 
pointed out by Showman and Guillot (2002), Menou 
et al. (2003), and Cho et al. (2003) and, as described 
previously, results from the fact that the Rossby de- 
formation radius and Rhines scale are comparable to 
a planetary radius for conditions relevant to typical 
hot Jupiters. At least two specific mechanisms seem 
to be relevant. When the deformation radius is close 
to the planetary radius, the day-night forcing tends 
to inject energy at horizontal scales comparable to 
the planetary radius; the /3 effect then anisotropizes 
this energy into zonal jets whose widths are the or- 
der of a planetary radius. Alternately, if the energy 
injection occurs primarily at small horizontal scales, 
an inverse cascade can reorganize the energy into 
planetary-scale jets if the planetary rotation is suf- 
ficiently slow and friction sufficiently weak (allow- 
ing fast winds and a large Rhines scale). Follow-up 
studies are consistent with these general expectations 
(Cooper and Showman 2005; Cho et al. 2008; Lang- 
ton and Laughhn 2007, 2008a; Dobbs-Dixon and Lin 
2008; Showman et al. 2008a, 2009), but additional 
work to clarify the relative roles of L]j and in set- 
ting the jet widths would be beneficial. 

Despite these similarities, there remain some key dif- 
ferences in the numerical results obtained by the various 
groups so far. These include the following: 

• As yet, little agreement exists on whether hot Jupiters 

should exhibit significant temporal variability on the 
global scale. The 2D simulations of Cho et al. (2003, 
2008) and Langton and Laughhn (2008a) produce 
highly turbulent, time-variable flows, whereas the 
shallow-water study of Langton and Laughhn (2007) 
and most of the 3D studies (Showman and Guillot 
2002; Cooper and Showman 2005; Dobbs-Dixon and 
Lin 2008; Showman et al. 2008a, 2009) exhibit rel- 
atively steady flow patterns with only modest time 
variability for zero-obliquity, zero-eccentricity hot 



"The similarity of the numerical estimates from Eq. (48) and (49) results 
from the fact that we assumed the same horizontal temperature differences 
in longitude and latitude and that Ro ~ 1 for the parameter regime ex- 
plored here. 



Jupiters. The Cho et al. simulations essentially guar- 
antee a time-variable outcome because of the tur- 
bulent initial condition, although breaking Rossby 
waves and dynamical instabilities during the simula- 
tions also play a large role; the Langton and Laugh- 
lin (2008a) simulations were initialized from rest and 
naturally develop turbulence via dynamical instabil- 
ity throughout the course of the simulation. 

Time variability can result from several mechanisms, 
including dynamical (e.g., barotropic or baroclinic) 
instabilities, large-scale oscillations, and modifica- 
tion of the jet pattern by atmospheric waves. If 
the jet pattem is dynamically stable (or, more rig- 
orously, if relevant instability growth times are sig- 
nificantly longer than the characteristic atmospheric 
heating/cooling times), then this mechanism for pro- 
ducing variabihty would be inhibited; on the other 
hand, if the jet pattern is unstable (with instability 
growth rates comparable to or less than the heat- 
ing/coohng times), the jets will naturally break up 
and produce a wide range of turbulent eddies. In this 
regard, it is crucial to carefully distinguish between 
2D and 3D models, because they exhibit very dif- 
ferent stability criteria (e.g., Dowling 1995). For ex- 
ample, in a 2D, nondivergent model, theoretical and 
numerical work shows that jets are guaranteed stable 
only if 

o^<^ <»' 

a relationship called the barotropic stability criterion. 
Numerical simulations show that 2D, nondivergent 
fluids initiahzed with jets violating Eq. (50) gener- 
ally develop instabilities that rob energy from the jets 
until the jets no longer violate the equation. However, 
jets in a 3D fluid can in some cases strongly violate 
Eq. (50) while remaining stable; Jupiter and Saturn, 
for example, have several jets where d'^u/ dy^ exceed 
/3 by a factor of 2-3 (IngersoU et al. 1981; Sanchez- 
Lavega et al. 2000), and yet the jet pattern has re- 
mained almost unchanged over multi-decade time 
scales. This may simply imply that the barotropic 
stability criterion is irrelevant for a 3D fluid (Dowl- 
ing 1995; Vasavada and Showman 2005). Thus, it is 
a priori imsurprising that 2D and 3D models would 
make different predictions for time variability. 

However, a recent study by Menou and Rauscher 
(2009) shows that the issue is not only one of 2D 
versus 3D. Making straightforward extensions to an 
Earth model, they performed forced 3D shallow sim- 
ulations initialized from rest that developed highly 
turbulent, time variable flows on the global scale, 
showing that time variability from the emergence of 
horizontal shear instability may occur in 3D under 
some conditions relevant to hot Jupiters. Like several 
previous studies, they forced their flows with a New- 
tonian heating/cooling scheme that generates a day- 



29 



night temperature difference. However, their forc- 
ing set-up differs significantly from those of previ- 
ous studies. They place an impermeable surface at 1 
bar and apply their maximum day/night forcing im- 
mediately above the surface; the day-night radiative- 
equilibrium temperature difference decreases with al- 
titude and reaches zero near the top of their model. 
This is the reverse of the set-up used in Cooper and 
Showman (2005, 2006) and Showman et al. (2008a, 
2009) where the forcing amphtude peaks near the 
top of the domain and decreases with depth, and 
where the surface is placed significantly deeper than 
the region of atmospheric heating/cooling to mini- 
mize its interaction with the circulation in the ob- 
servable atmosphere. The differences in the forcing 
schemes and position of the surface presumably de- 
termine whether or not a given 3D simulation devel- 
ops global-scale instability and time variabihty. Fur- 
ther work will be necessary to identify the specific 
forcing conditions that lead to time-variable or quasi- 
steady conditions at the global scale. 

It is also worth remembering that hot Jupiters them- 
selves exhibit huge diversity; among the known tran- 
siting planets, for example, incident stellar flux varies 
by a factor of ~ 150, gravities range over a factor of 
70, and expected rotation rates vary over nearly 
a factor of 10 (e.g.. Table 1). Given this diversity, 
it is reasonable to expect that some hot Jupiters will 
exhibit time- variable conditions while others will ex- 
hibit relatively steady flow patterns. 

• Another area of difference between models regards 
the direction of the equatorial jet. The shallow- water 
and equivalent barotropic models can produce ei- 
ther eastward or westward equatorial jets, depend- 
ing on the initial conditions and other modeling de- 
tails (Cho et al. 2003, 2008; Langton and LaughUn 
2007, 2008a). In contrast, all of the 3D models pub- 
hshed to date for synchronously rotating hot Jupiters 
with zero obliquity and eccentricity have produced 
robust eastward equatorial jets (Showman and Guil- 
lot 2002; Cooper and Showman 2005, 2006; Show- 
man et al. 2008a, 2009; Dobbs-Dixon and Lin 2008; 
Menou and Rauscher 2009). The mechanisms re- 
sponsible for generating the eastward jets in these 3D 
models remain to be diagnosed in detail but presum- 
ably involve the global-scale eddies generated by the 
day-night heating contrast. 

To date, most work on hot- Jupiter atmospheric circula- 
tion has emphasized planets on circular orbits. However, 
several transiting hot Jupiters and Neptunes have eccentric 
orbits, including GJ 436b, HAT-P-2b, HD 17156b, and HD 
80606b, whose orbital eccentricities are 0.15, 0.5, 0.67, and 
0.93, respectively. When the eccentricity is large, the inci- 
dent stellar flux varies significantly throughout the orbit — 
HAT-P-2b, for example, receives 9 times more flux at pe- 



riapse than at apoapse, while for HD 80606b, the maximum 
incident flux exceeds the minimum incident flux by a fac- 
tor of over 800! These variations dwarf those experienced 
in our Solar System and constitute an uncharted regime of 
dynamical forcing for a planetary atmosphere. An 8-/xm 
fight curve of HD 80606b has recently been obtained with 
the Spitzer Space Telescope during a 30-hour interval sur- 
rounding periapse passage (Laughlin et al. 2009), which 
shows the increase in planetary flux that presumably ac- 
companies the flash heating as the planet passes by its star. 
2D and 3D numerical simulations suggest that these plan- 
ets may exhibit dynamic, time- variable flows and show that 
the planet's emitted infrared flux can peak hours to days af- 
ter periapse passage (Langton and Laughlin 2008a,b; Lewis 
et al. 2009). 

4.4. Chemistry as a probe of the meteorology 

Chemistry can provide important constraints on the me- 
teorology of giant planets. On Jupiter, for example, gaseous 
CO, PH3, GeH4, and ASH3 have been detected in the upper 
troposphere (p < 5 bars) with mole fractions of ~ 0.8, 0.2, 
0.6, and 0.5-1 ppb, respectively. None of these compounds 
are thermochemically stable in the upper troposphere — the 
chemical equiUbrium mole fractions are < 10"^'' for ASH3 
and < 10~^° for the other three species (Fegley and Lod- 
ders 1994). However, the equilibrium abundances of all 
four species rise with depth, exceeding 1 ppb at pressures 
greater tiian 500, 400, 160, and 20 bars for CO, PH3, GeH4, 
and ASH3, respectively. Thus, the high abundances of these 
species in the upper troposphere appear to result from rapid 
convective mixing from the deep atmosphere.^^ Kinetic re- 
action times are rapid at depth but plummet with decreasing 
temperature and pressure, allowing these species to persist 
in disequiUbrium (i.e., "quench") in the low-pressure and 
temperature conditions of the upper troposphere. Knowl- 
edge of the chemical kinetic rate constants allows quantita- 
tive estimates to be made of the vertical mixing rate needed 
to explain the observed abundances. For all four of the 
species listed above, the required vertical mixing rates yield 
order-of-magnitude matches with the mixing rate expected 
from thermal convection at Jupiter's known heat flux (e.g. 
Prinn and Barshay 1977; Fegley and Lodders 1994; Bezard 
et al. 2002). Thus, the existence of these disequiUbrium 
species provides evidence that the deep atmosphere is in- 
deed convective. 

Given the difficulty of characterizing the meteorology of 
exoplanets with spatially resolved observations, these types 
of chemical constraints will hkely play an even more im- 
portant role in our understanding of exoplanets than is the 
case for Jupiter. Several brown dwarfs, including Gl 229b, 
Gl 570d, and 2MASS J0415-0935, show evidence for dise- 
quilibrium abundances of CO and/or NH3 (e.g., Noll et al. 
1997; Saumon et al. 2000, 2006, 2007). These observations 
can also plausibly be explained by vertical mixing. For 



'Jupiter also has a high stratospheric CO abundance that is inferred to result 
from exogenous sources (Bezard et al. 2002). 
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these objects, NH3 is quenched in the convection zone, but 
CO is quenched in the overlying radiative zone where ver- 
tical mixing must result from some combination of small- 
scale turbulence (due, for example, to the breaking of verti- 
cally propagating gravity waves) and vertical advection due 
to the large-scale circulation (§3.7). Current models sug- 
gest that such mixing corresponds to an eddy diffusivity 
in the radiative zone of ^ 1-10 scc^^ for Gl 229b and 
10-100 m^sec-i for 2MASS J0415-0935 (Saumon et al. 
2007). These inferences on vertical mixing rate provide 
constraints on the atmospheric circulation that would be dif- 
ficult to obtain in any other way. 

Several hot Jupiters likewise reside in temperature 
regimes where CH4 is the stable form in the observable 
atmosphere and CO is the stable form at depth, or where 
CH4 is the stable form on the nightside and CO is the stable 
form on the dayside. The abundance and spatial distribution 
of CO and CH4 on these objects could thus provide impor- 
tant clues about the wind speeds and other aspects of the 
meteorology. CH4 has been detected via transmission spec- 
troscopy on HD 189733b (Swain et al. 2008) and CO has 
been indirectly inferred from the shape of the dayside emis- 
sion spectrum (Barman 2008), although further observa- 
tions and radiative-transfer models are needed to determine 
their abundances and spatial variability (if any). Cooper 
and Showman (2006) coupled a simple model for CO/CH4 
interconversion kinetics to their 3D dynamical model of 
hot Jupiters to demonstrate that CO and CH4 should be 
quenched in the observable atmosphere, where thermo- 
chemical CO«-»CH4 interconversion times are orders of 
magnitude longer than plausible dynamical times. If so, one 
would expect CO and CH4 to have similar abundances on 
the dayside and nightside of hot Jupiters, in contrast to the 
chemical-equilibrium situation where CH4 would be much 
more prevalent on the nightside (note, however, that their 
models neglect photochemistry, which could further modify 
the abundances of both species). For temperatures relevant 
to objects like HD 209458b and HD 189733b, Cooper and 
Showman (2006) found that the quenching occurs at pres- 
sures of ~ 1-10 bars, and that the quenched abundance 
depends primarily on the temperatures and vertical veloci- 
ties in this pressure range. 

Thus, future observational determination of CO and CH4 
abundances and comparison with detailed 3D models may 
allow constraints on the temperatures and vertical veloci- 
ties in the deep (1-10 bar) atmospheres of hot Jupiters to 
be obtained. Such constraints probe significantly deeper 
than the expected visible and infrared photospheres of hot 
Jupiters and hence would nicely complement the charac- 
terization of hot- Jupiter thermal structure via light curves 
and secondary-eclipse spectra. Note also that an analogous 
story involving N2 and NH3 interconversion may occur on 
even cooler objects (Tgff 700 K "warm" Jupiters). 



5. CIRCULATION REGIMES: TERRESTRIAL PLAN- 
ETS 

With current groundbased searches and spacecraft such 
as NASA's Kepler mission, terrestrial planets are likely to 
be discovered around main-sequence stars within the next 
five years, and basic characterization of their atmospheric 
composition and temperature structure should follow over 
the subsequent decade. At base, we wish to understand not 
only whether such planets exist but whether they have at- 
mospheres, what their atmospheric composition, structure 
and climate may be, and whether they can support life. Ad- 
dressing these issues will require a consideration of rele- 
vant physical and chemical cUmate feedbacks as well as the 
plausible circulation regimes in these planets' atmospheres. 

The study of terrestrial exoplanet atmospheres is just be- 
ginning, and only a handful of papers have specifically in- 
vestigated the possible circulation regimes on these objects. 
However, a vast literature has developed to understand the 
climate and circulation of Venus, Titan, Mars, and espe- 
cially Earth, and many of the insights developed for under- 
standing these planets may be generalizable to exoplanets. 
Our goal here is to provide conceptual and theoretical guid- 
ance on the types of climate and circulation processes that 
exist in terrestrial planet atmospheres and discuss how those 
processes vary under diverse planetary conditions. In §5.1 
we survey basic issues in climate. The following sections, 
§5.2-§5.5, address basic circulation regimes on terrestrial 
planets, emphasizing those aspects (e.g., processes that de- 
termine the horizontal temperature contrasts) relevant to 
future exoplanet observations. §5.6 discusses regimes of 
exotic forcing associated with synchronous rotation, large 
obliquities, and large orbital eccentricities. 

5.1. Climate 

Climate can be defined as the mean condition of a 
planet's atmosphere/ocean system — the temperatures, pres- 
sures, winds, humidities, and cloud properties — averaged 
over time intervals longer than the timescale of typical 
weather events. The climate on the terrestrial planets re- 
sults from a wealth of interacting physical, chemical, ge- 
ological, and (when relevant) biological effects, and even 
on Earth, understanding the past and present climate has 
required a multi-decade interdisciplinary research effort. In 
this section, we provide only a brief sampling of the sub- 
ject, touching only on the most basic physical processes 
that help to determine a planet's global-mean conditions. 

The mean temperatures at which a planet radiates to 
space depend primarily on the incident stellar flux and 
the planetary albedo. Equating emitted infrared energy 
(assumed blackbody) with absorbed stellar flux yields 
a planet's global-mean effective temperature at radiative 
equilibrium: 



^eff 



4cr 



1/4 



(51) 



where is the incident stellar flux, is the planet's 
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global-mean bond albedo^^, and a is the Stefan-Boltzmann 
constant. The factor of 4 results from the fact that the planet 
intercepts a stellar beam of area i^o? but radiates infrared 
from its full surface area of 47ra^. 

Given an incident stellar flux, a variety of atmospheric 
processes act to determine the planet's surface temperature. 
First, the circulation and various atmospheric feedbacks can 
play a key role in determining the planetary albedo. Bare 
rock is fairly dark, and for terrestrial planets the albedo is 
largely determined by the distribution of clouds and surface 
ice. The Moon, for example, has bond albedo of 0.1 1, yield- 
ing an effective temperature of 274 K. In contrast, Venus 
has a bond albedo of 0.75 because of its global cloud cover, 
yielding an effective temperature of 232 K. With partial 
cloud and ice cover. Earth is an intermediate case, with a 
bond albedo of 0.31 and an effective temperature of 255 K. 
These examples illustrate that the circulation, via its effect 
on clouds and surface ice, can have a major influence on 
mean conditions — ^Venus' effective temperature is less than 
that of the Earth and Moon despite receiving nearly double 
the Solar flux! 

Second, Teg is not the surface temperature but the 
mean blackbody temperature at which the planet radiates 
to space. The surface temperatures can significantly exceed 
Teff through the greenhouse effect. Planets orbiting Sun- 
like stars receive most of their energy in the visible, but 
they tend to radiate their energy to space in the infrared. 
Most gases tend to be relatively transparent in the visible, 
but H2O, CO2, CH4, and other molecules absorb signif- 
icantly in the infrared. On a planet hke Earth or Venus, 
a substantial fraction of the sunUght therefore reaches the 
planetary surface, but infrared emission from the surface 
is mostly absorbed in the atmosphere and cannot radiate 
directly to space. In turn, the atmosphere radiates both up 
(to space) and down (to the surface). The surface therefore 
receives a double whammy of radiation from both the Sun 
and the atmosphere; to achieve energy balance, the surface 
temperature becomes elevated relative to the effective tem- 
perature. This is the greenhouse effect. Contrary to popular 
descriptions, the greenhouse effect should not be thought 
of as a situation where the heat is "trapped" or "cannot 
escape." Indeed, the planet as a whole resides in a near- 
balance where infrared emission to space (mostly from the 
atmosphere when the greenhouse effect is strong) almost 
equals absorbed sunlight. 

To have a significant greenhouse effect, a planet must 
have a massive-enough atmosphere to experience pressure 
broadening of the spectral lines. Mars, for example, has 
~ 15 times more CO2 per area than Earth, yet its green- 
house effect is significantly weaker because its surface pres- 
sure is only 6 mbar. 

The cUmate is influenced by numerous feedbacks that 
can affect the mean state and how the atmosphere responds 
to perburbations. Some of the more important feedbacks 



^'^The bond albedo is the fraction of light incident upon a planet that is scat- 
tered back to space, integrated over all wavelengths and directions. 



are as follows: 

• Thermal feedback: Increases or decreases in the tem- 
perature at the infrared photosphere lead to enhanced 
or reduced radiation to space, respectively. This is a 
negative feeback that allows planets to reach a stable 
equiUbrium with absorbed starUght. 

• Ice-albedo feedback: Surface ice can form on plan- 
ets exhibiting trace gases that can condense to solid 
form. Because of the brightness of ice and snow, an 
increase in snow/ice coverage decreases the absorbed 
starlight, promoting colder conditions and growth of 
even more ice. Conversely, melting of surface ice 
increases the absorbed starlight, promoting warmer 
conditions and continued melting. This is a positive 
feedback. 

Simple models of the ice-albedo feedback show that, 
over a range of values (depending on the strength 
of the greenhouse effect), an Earth-like planet can ex- 
hibit two stable equilibrium states for a given value 
of F*: a warm, ice-poor state with a low albedo and 
a cold, ice-covered state with a high albedo (e.g.. 
North et al. 1981). Geologic evidence suggests that 
~ 0.6-2.4 Gyr ago Earth experienced several multi- 
Myr-long glaciations with global or near-global ice 
cover (dubbed "snowball Earth" events; Hoffman and 
Schrag 2002), suggesting that Earth has flipped back 
and forth between these equilibria. The susceptibility 
of a planet to entering such a snowball state is highly 
sensitive to the strength of latitudinal heat transport 
(Spiegel et al. 2008), thereby Unking this feedback to 
the global circulation. 

• Condensable-greenhouse-gas feedback: When an at- 
mospheric constituent is a greenhouse gas that also 
exists in condensed form on the surface, the atmo- 
sphere can experience a positive feedback that af- 
fects the temperature and the distribution of this con- 
stituent between the atmosphere and surface. An 
increase in surface temperature increases the con- 
stituent's saturation vapor pressure, hence the atmo- 
spheric abundance and therefore the greenhouse ef- 
fect. A decrease in surface temperature decreases the 
saturation vapor pressure, reducing the atmospheric 
abundance and therefore the greenhouse effect. Ther- 
mal perturbations are therefore amplified. For Earth, 
this process occurs with water vapor (Earth's most 
important greenhouse gas) and is called the "water- 
vapor feedback" in the Earth climate literature (Held 
and Soden 2000). The water- vapor feedback will 
play a major role in determining how Earth responds 
to anthropogenic increases in carbon dioxide over the 
next century. 

In some circumstances, this positive feedback is 
so strong that it can trigger a runaway that shifts 
the atmosphere into a drastically different state. In 
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the case of warming, this would constitute a run- 
away greenhouse that leads to the complete evapora- 
tion/subUmation of the condensable constituent from 
the surface. If early Venus had oceans, for example, 
they might have experienced runaway evaporation, 
leading to a monstrous early water vapor atmosphere 
(Ingersoll 1969; Kasting 1988) that would have had 
major effects on subsequent planetary evolution. 

In the case of cooling, such a runaway would remove 
most of the condensable constituent from the atmo- 
sphere, and if the relevant gas dominates the atmo- 
sphere, this could lead to atmospheric collapse. For 
example, if a Venus-like planet (with its ~ 500-K 
greenhouse effect) were moved sufficiently far from 
the Sun, CO2 condensation would initiate, potentially 
collapsing the atmosphere to a Mars-hke state with 
most of the CO2 condensed on the ground, a cold, 
thin CO2 atmosphere in vapor-pressure equilibrium 
with the smface ice, and minimal greenhouse effect. 
The CO2 cloud formation that precedes such a col- 
lapse is often used to define the outer edge of the 
classical habitable zone for terrestrial planets whose 
greenhouse effect comes primarily from CO2 (Kast- 
ing etal. 1993). 

Because CO2 condensation naturally initiates in 
the coldest regions (the poles for a rapidly rotating 
planet; the nightside for a synchronous rotator), the 
conditions under which atmospheric coUapse initi- 
ates depend on the atmospheric circulation. Weak 
equator-to-pole (or day-night) heat transport leads to 
colder polar (or nightside) temperatures for a given 
solar flux, promoting atmospheric coUapse. These 
dynamical effects have yet to be fuUy included in 
climate models. 

• Clouds: Cloud coverage increases the albedo, lessen- 
ing the absorption of starlight and promoting cooler 
conditions. On the other hand, because the tro- 
pospheric temperatures generally decrease with alti- 
tude, cloud tops are typically cooler than the ground 
and therefore radiate less infrared energy to space, 
promoting warmer conditions (an effect that depends 
sensitively on cloud altitude and latitude). These ef- 
fects compete with each other. By determining the 
detailed properties of clouds (e.g., fractional cover- 
age, latitudes, altitudes, and size-particle distribu- 
tions), the atmospheric circulation therefore plays a 
major role in determining the mean surface tempera- 
ture. 

To the extent that cloud properties depend on atmo- 
spheric temperature, clouds can act as a feedback that 
affects the mean state and ampUfies or reduces a ther- 
mal perturbation to the climate system. However, be- 
cause the sensitivity of cloud properties to the global 
circulation and mean climate is extremely difficult to 
predict, the net sign of this feedback (positive or neg- 



ative) remains unknown even for Earth. For exoplan- 
ets, clouds could plausibly act as a positive feedback 
in some cases and negative feedback in others. 

• Long-term atmosphere/geology feedbacks: On geo- 
logical timescales, terrestrial planets can experience 
significant exchange of material between the interior 
and atmosphere. An example that may be particu- 
larly relevant for Earth-like (i.e. ocean and continent- 
bearing) exoplanets is the carbonate-silicate cycle, 
which can potentially buffer atmospheric CO2 in a 
temperature-dependent way that tends to stabilize the 
atmospheric temperature against variations in solar 
flux (see, e.g. Kasting and Catling 2003). Such feed- 
backs, while beyond the scope of this chapter, may 
be critical in determining the mean composition and 
temperature of terrestrial exoplanets. 

In addition to the major feedbacks outlined above, there 
exist dozens of additional interacting physical, chemical, 
and biological feedbacks that can influence the mean cli- 
mate and its sensitivity to perturbations. Hartmann and 
coauthors (2003) provide a thorough assessment for the 
Earth's current climate; although the details wiU be differ- 
ent for other planets, this gives a flavor for the complexity 
that can be expected. 

5.2. Global circulation regimes 

For atmospheres with radiative time constants greatly 
exceeding the planet's solar day, the atmosphere cannot 
respond rapidly to day-night variations in stellar heating 
and instead responds primarily to the daily-mean insolation, 
which is a function of latitude. This situation applies to 
most planets in the Solar System, including the (lower) at- 
mospheres of Venus, Earth, Titan, Jupiter, Satum, Uranus, 
and Neptune.^^ The resulting patterns of temperature and 
winds vary little with longitude in comparision to their vari- 
ation in latitude and height. In such an atmosphere, the pri- 
mary task of the atmospheric circulation is to transport ther- 
mal energy not from day to night but between the equator 
and the poles. On the other hand, when the atmosphere's 
radiative time constant is much shorter than the solar day, 
the day-night heating gradient will be paramount and day- 
night temperature variations at the equator could potentially 
rival temperature variations between the equator and poles. 

In the next several subsections we survey our under- 
standing of global atmospheric circulations for these two 
regimes. Because the first situation applies to most Solar- 
System atmospheres, it has received the vast majority of 
work and remains better understood. This regime, which 



'Mars is a transitional case, with a radiative time constant ~ 2-3 times its 
24.6-hour day. 

'At low obliquity (< 51°), yearly averaged starlight is absorbed primarily 
at low latitudes, so the circulation transports thermal energy poleward, but 
at high obliquity (> 54°), yearly averaged starlight is absorbed primarily 
at the poles (Ward 1974), and the yearly averaged energy transport by the 
circulation is equatorward. At high obliquity, strong seasonal cycles will 
also occur. 
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we discuss in §§5.3-5.5, will apply to planets whose or- 
bits are sufficiently far from their stars that the planets have 
not despun into a synchronously rotating state; it can also 
apply to atmosphere-bearing moons of hot Jupiters even at 
small distances from their stars. The latter regime, sur- 
veyed in §5.6, prevails for planets exhibiting small enough 
orbital eccentricities and semi-major axes to become syn- 
chronously locked to their stars. This situation probably 
apphes to most currently known transiting giant exoplan- 
ets and may also apply to terrestrial planets in the habitable 
zones of M dwarfs, which are the subject of current obser- 
vational searches. Despite its current relevance, this novel 
forcing regime has come under investigation only in the past 
decade and remains incompletely understood due to a lack 
of Solar-System analogs. 

5.3. Axisymmetric flows: Hadley cells 

Perhaps the simplest possible idealization of a circula- 
tion that transports heat from equator to poles is an ax- 
isymmetric circulation — that is, a circulation that is inde- 
pendent of longitude — where hot air rises at the equator, 
moves poleward aloft, cools, sinks at high latitudes, and re- 
turns equatorward at depth (near the surface on a terrestrial 
planet). Such a circulation is termed a //a<i/e>' cell, and was 
first envisioned by Hadley in 1735 to explain Earth's trade 
winds. Most planetary atmospheres in our Solar System, 
including those of Venus, Earth, Mars, Titan, and possibly 
the giant planets, exhibit Hadley circulations. 

Hadley circulations on real planets are of course not truly 
axisymmetric; on the terrestrial planets, longitudinal varia- 
tions in topography and thermal properties (e.g., associated 
with continent-ocean contrasts) induce asymmetry in longi- 
tude. Nevertheless, the fundamental idea is that the longitu- 
dinal variations are not crucial for driving the circulation. 
This differs from the circulation in midlatitudes, whose 
longitudinally averaged properties are fundamentally con- 
trolled by the existence of non-axisymmetric baroclinic ed- 
dies that are inherently three-dimensional (see §3.7). 

Planetary rotation generally prevents Hadley circulations 
from extending all the way to the poles. Because of plan- 
etary rotation, equatorial air contains considerable angular 
momentum about the planetary rotation axis; to conserve 
angular momentum, equatorial air would accelerate to un- 
realistically high speeds as it approached the pole, a phe- 
nomenon which is dynamically inhibited. To illustrate, the 
specific angular momentum about the rotation axis on a 
spherical planet is M = (fiacos^ + u)a cos cj), where the 
first and second terms represent angular momentum due to 
planetary rotation and winds, respectively (recall that fl, a, 
and (p are planetary rotation rate, planetary radius, and lat- 
itude). If w = at the equator, then M = Qa^, and an 
angular-momentum conserving circulation would then ex- 
hibit winds of 

u = na'^ (52) 

cos 

Given Earth's radius and rotation rate, this equation im- 



plies zonal-wind speeds of 134 msec at 30° latitude, 
700 m sec" 1 at 60° latitude, and 2.7 km sec" ^ at 80° lat- 
itude. Such high-latitude wind speeds are unrealistically 
high and would furthermore be violently imstable to 3D in- 
stabilities. On Earth, the actual Hadley circulations extend 
to ~ 30° latitude. 

The Hadley circulation exerts strong control over the 
wind structure, latitudinal temperature contrast, and cli- 
mate. Hadley circulations transport thermal energy by the 
most efficient means possible, namely straightforward ad- 
vection of air from one latitude to another. As a result, the 
latitudinal temperature contrast across a Hadley circulation 
tends to be modest; the equator-to-pole temperature con- 
trast on a planet will therefore depend strongly on the width 
of the Hadley cell. Moreover, on planets with condensable 
gases, Hadley cells exert control over the patterns of cloudi- 
ness and rainfall. On Earth, the rising branch of the Hadley 
circulation leads to cloud formation and abundant rainfall 
near the equator, helping to explain for example the preva- 
lence of tropical rainforests in Southeast Asia/Indonesia, 
Brazil, and central Africa.^'^ On the other hand, because 
condensation and rainout dehydrates the rising air, the de- 
scending branch of the Hadley cell is relatively dry, which 
explains the abundances of arid chmates on Earth at 20-30° 
latitude, including the deserts of the African Sahara, South 
Africa, Australia, central Asia, and the southwestern United 
States. The Hadley cell can also influence the mean cloudi- 
ness, hence albedo and thereby the mean surface tempera- 
ture. Venus's slow rotation rate leads to a global equator- 
to-pole Hadley cell, with a broad ascending branch in low 
latitudes. Coupled with the presence of trace condensable 
gases, this widespread ascent contributes to a near-global 
cloud layer that helps generate Venus' high bond albedo 
of 0.75. Different Hadley cell pattems would presumably 
cause different cloudiness patterns, different albedos, and 
therefore different global-mean surface temperatures. On 
exoplanets, Hadley circulations wiU also likewise help con- 
trol latitudinal heat fluxes, equator-to-pole temperature con- 
trasts, and chmate. 

A variety of studies have been carried out using fully 
nonlinear, global 3D numerical circulation models to de- 
termine the sensitivity of the Hadley cell to the plane- 
tary rotation rate and other parameters (e.g. Hunt 1979; 
Wifliams and Hofloway 1982; Williams 1988a,b; Del Genio 
and Suozzo 1987; Navarra and Boccaletti 2002; Walker and 
Schneider 2005, 2006). These studies show that as the rota- 
tion rate is decreased the width of the Hadley cell increases, 
the equator-to-pole heat flux increases, and the equator-to- 
pole temperature contrast decreases. Figs. 11-12 illustrates 
examples from Navarra and Boccaletti (2002) and Del Ge- 
nio and Suozzo (1987). For Earth parameters, the circula- 
tion exhibits mid-latitude eastward jet streams that peak in 
the upper troposphere (^ 200 mbar pressure), with weaker 
wind at the equator (Fig. 1 1). The Hadley cells extend from 
the equator to the equatorward flanks of the mid-latitude 



'Regional circulations, such as monsoons, also contribute. 



34 




Fig. 1 1 . — Left: Zonal-mean circulation versus latitude in degrees (abscissa) and pressure in mbar (ordinate) in a series of Earth-based 
GCM experiments from Navarra and Boccaletti (2002) where the rotation period is varied from 18 hours (top) to 360 hours (bottom). 
Greyscale and thin grey contours depict zonal-mean zonal wind, u, and thick black contours denote streamfunction of the circulation in 
the latitude-height plane (the so-called "meridional circulation"), with solid being clockwise and dashed being coimterclockwise. The 
meridional circulation flows parallel to the streamfunction contours, with greater mass flux when contours are more closely spaced. The 
two cells closest to the equator correspond to the Hadloy cell. As rotation period increases, the jets move poleward and the Hadley cell 
widens, becoming nearly global at the longest rotation periods. 
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jets. As the rotation rate decreases, the Hadley cells widen 
and the jets shift poleward. At first, the jet speeds increase 
with decreasing rotation rate, which results from the fact 
that as the Hadley cells extend poleward (i.e., closer to the 
rotation axis) the air can spin-up faster (cf Eq. 52). Even- 
tually, once the Hadley cells extend almost to the pole (at 
rotation periods exceeding ~ 5-10 days for Earth radius, 
gravity, and vertical thermal structure), further decreases in 
rotation rate reduce the mid-latitude jet speed. 

Perhaps more interestingly for exoplanet observations, 
these changes in the Hadley cell significantly influence 
the planetary temperature structure. This is illustrated in 
Fig. 12 from a series of simulations by Del Genio and 
Suozzo (1987). Because the Hadley cells transport heat 
extremely efficiently, the temperature remains fairly con- 
stant across the width of the Hadley cells. Poleward of the 
Hadley cells, however, the heat is transported in latitude 
by barochnic instabilities (§§3.7 and 5.4), which are less 
efficient, so a large latitudinal temperature gradient exists 
within this so-called "baroclinic zone." The equator-to-pole 
temperature contrast depends strongly on the width of the 
Hadley cell. 

Despite the value of the 3D circulation models described 
above, the complexity of these models tends to obscure the 
physical mechanisms goveming the Hadley circulation's 
strength and latitudinal extent and cannot easily be extrapo- 
lated to different planetary parameters. A conceptual theory 
for the Hadley cell, due to Held and Hou (1980), provides 
considerable insight into Hadley cell dynamics and allows 
estimates of how, for example, the width of the Hadley cell 
should scale with planetary size and rotation rate [see re- 
views in James (1994, pp. 80-92), ValHs (2006, pp. 457- 
466), and Schneider (2006)]. Stripped to its basics, the 
scheme envisions an axisymmetric two-layer model, where 
the lower layer represents the equatorward flow near the 
surface and the upper layer represents the poleward flow 
in the upper troposphere. For simplicity. Held and Hou 
(1980) adopted a basic-state density that is constant with 
altitude. Absorption of sunlight and loss of heat to space 
generate a latitudinal temperature contrast that drives the 
circulation; for concreteness, let us parameterize the radi- 
ation as a relaxation toward a radiative-equilibrium poten- 
tial temperature profile that varies with latitude as ^rad = 
9o — A(?i.ad sin^ (p, where is the radiative-equilibrium po- 
tential temperature at the equator and A^rad is the equator- 
to-pole difference in radiative-equilibrium potential temper- 
ature. If we make the small-angle approximation for sim- 
phcity (valid for a Hadley cell that is confined to low lati- 
tudes), we can express this as ^rad = — ^Orad'P'^- 

In the lower layer, we assume that friction against the 
ground keeps the wind speeds low; in the upper layer, as- 
sumed to occur at an altitude H, the flow conserves angu- 
lar momentum. The upper layer flow is then specified by 
Eq. (52), which is just u = U.acj)^ in the small-angle hmit. 
We expect that the upper-layer wind will be in thermal-wind 



balance with the latitudinal temperature contrast^': 
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where du/dz is simply given by u/H in this two-layer 
model. Inserting u = Qa(j)^ into Eq. (53), approximat- 
ing the Coriolis parameter as / = /3y (where /3 is treated 
as constant), and integrating, we obtain a temperature that 
varies with latitude as 



e = e, 



equator 



(54) 



where ^equator is a constant to be determined. 

At this point, we introduce two constraints. First, Held 
and Hou (1980) assumed the circulation is energetically 
closed, i.e. that no net exchange of mass or thermal energy 
occurs between the Hadley cell and higher latitude circula- 
tions. Given an energy equation with radiation parameter- 
ized using Newtonian cooling, dO/dt = (^rad — ^)/Trad. 
where Tiad is a radiative time constant, the assumption that 
the circulation is steady and closed requires that 
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where we are integrating from the equator to the poleward 
edge of the Hadley cell, at latitude (pn- Second, tempera- 
ture must be continuous with latitude at the poleward edge 
of the Hadley cell. In the axisymmetric model, baroclinic 
instabilities are suppressed, and the regions poleward of the 
Hadley cells reside in a state of radiative equihbrium. Thus, 
9 must equal ^rad at the poleward edge of the cell. Insert- 
ing our expressions for 9 and 6'rad into these two constraints 
yields a system of two equations for (pH and ^equator- The 
solution yields a Hadley cell with a latitudinal half-width of 
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in radians. This solution suggests that the width of the 
Hadley cell scales as the square root of the fractional 
equator-to-pole radiative-equilibrium temperature differ- 
ence, the square root of the gravity, the square root of 
the height of the cell, and inversely with the rotation 
rate. Inserting Earth annual-mean values (A^^ad ~ 70 K, 
61 = 260K, 5 = 9.8msec-2, H = 15km, a = 6400km, 
and n = 7.2x 10'^ scc'^) yields - 30°. 

Redoing this analysis without the small-angle approx- 
imation leads to a transcendental equation for (pH (see 
Eq. (17) in Held and Hou 1980), which can be solved nu- 
merically. Figure 13 (solid curve) illustrates the solution. 
As expected, (pH ranges from 0° as f2 — > oo to 90° as 

^ 0, and, for planets of Earth radius with Hadley cir- 
culations ~ 10 km tall, bridges these extremes between ro- 
tation periods of ~ 0.5 and 20 days. Although deviations 



^^This form differs slightly from Eq. (22) because Eq. (53) adopts a constant 
basic-state density (the so-called "Boussinesq" approximation) whereas 
Eq. (22) adopts the compressible ideal-gas equation of state. 
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Fig. 12. — Zonally averaged temperature in °C versus latitude at 500 mbar (in the mid-troposphere) for a sequence of Earth-like GCM 
runs that vary the planetary rotation period between 1 and 64 Earth days (labelled in the graph). Latitude runs from 90° on the left to 
—90° on the right, with the equator at the center of the horizontal axis. The flat region near the equator in each run results from the 
Hadley cell, which transports thermal energy extremely efficiently and leads to a nearly isothermal equatorial temperature structure. The 
region of steep temperature gradients at high latitudes is the baroclinic zone, where the temperature structure and latitudinal heat transport 
are controlled by eddies resulting from baroclinic instability. Note that the width of the Hadley cell increases, and the equator-to-pole 
temperature difference decreases, as the planetary rotation period is increased. From Del Genio and Suozzo (1987). 



exist, the agreement between the simple Held-Hou model 
and the 3D GCM simulations are surprisingly good given 
the simplicity of the Held-Hou model. 

Importantly, the Held-Hou model demonstrates that lat- 
itudinal confinement of the Hadley cell occurs even in an 
axisymmetric atmosphere. Thus, the cell's latitudinal con- 
finement does not require (for example) three-dimensional 
baroclinic or barotropic instabilities associated with the jet 
at the poleward branch of the cell. Instead, the confinement 
results from energetics: the twin constraints of angular- 
momentum conservation in the upper branch and thermal- 
wind balance specify the latitudinal temperature profile in 
the Hadley circulation (Eq. 54). This generates equato- 
rial temperatures colder than (and subtropical temperatures 
warmer than) the radiative equilibrium, implying radiative 
heating at the equator and cooling in the subtropics. This 
properly allows the circulation to transport thermal energy 
poleward. If the cell extended globally on a rapidly ro- 
tating planet, however, the circulation would additionally 
produce high-latitude temperatures colder than the radia- 
tive equilibrium temperature, which in steady state would 
require radiative heating at high latitudes. This is thermo- 
dynamically impossible given the specified latitudinal de- 
pendence of ^rad- The highest latitude to which the cell can 
extend without encountering this problem is simply given 
byEq. (55). 

The model can be generalized to consider a more re- 
alistic treatment of radiation than the simplified Newto- 
nian cooling/heating scheme employed by Held and Hou 
(1980). Caballero et al. (2008) reworked the scheme using a 
two-stream, non-grey representation of the radiative trans- 
fer with parameters appropriate for Earth and Mars. This 
leads to a prediction for the width of the Hadley cell that 
differs from Eq. (56) by a numerical constant of order unity. 



Although the prediction of Held-Hou-type models for 
(f)H provides impoitant insight, several failures of these 
models exist. First, the model underpredicts the strength 
of the Earth's Hadley cell (e.g., as characterized by the 
magnitude of the north-south wind) by about an order of 
magnitude. This seems to result from the lack of turbu- 
lent eddies in axisymmetric models; several studies have 
shown that turbulent three-dimensional eddies exert stresses 
that act to strengthen the Hadley cells beyond the predic- 
tions of axisymmetric models (e.g. Kim and Lee 2001; 
Walker and Schneider 2005, 2006; Schneider 2006).^^ 
ond, the Hadley cells on Earth and probably Mars are not 
energetically closed; rather, mid-latitude baroclinic eddies 
transport thermal energy out of the Hadley cell into the 
polar regions. Third, the poleward-moving upper tropo- 
spheric branches of the Hadley cells do not conserve an- 
gular momentum — although the zonal wind does become 
eastward as one moves poleward across the cell, for Earth 
this increase is a factor of ^ 2-3 less than predicted by 
Eq. (52). Overcoming these failings requires the inclusion 
of three-dimensional eddies. 

Several studies have shown that turbulent eddies in the 
mid- to high-latitudes — which are neglected in the Held- 
Hou and other axisymmetric models — can affect the width 
of the Hadley circulation and alter the parameter depen- 
dences suggested by Eq. (56) (e.g., Del Genio and Suozzo 
1987; Walker and Schneider 2005, 2006). Turbulence can 
produce an acceleration/deceleration of the zonal-mean 



'Held and Hou (1980)'s original model neglected tlie seasonal cycle, and 
it has been suggested that generalization of the Held-Hou axisymmetric 
model to include seasonal effects could alleviate this failing (Lindzen and 
Hou 1988). Although this improves the agreement with Earth's observed 
annual-mean Hadley-cell strength, it predicts large solstice/equinox os- 
cillations in Hadley-cell strength that are lacking in the observed Hadley 
circulation (Dima and Wallace 2003). 
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Fig. 13. — Latitudinal width of the Hadley cell from a sequence of Earth-like GCM runs that vary the planetary rotation period 
(symbols) and comparison to the Held and Hou (1980) theory (solid curve), the small-angle approximation to the Held-Hou theory from 
Eq. (56) (dotted curve), and the Held (2000) theory from Eq. (57) (dashed-dotted curve). GCM results are from Del Genio and Suozzo 
(1987) (squares and triangles, depicting northern and southern hemispheres, respectively), Navarra and Boccaletti (2002) (filled circles), 
and Korty and Schneider (2008) (asterisks; only cases adopting 0rad « 70K and their parameter F = 0.7 are shown). Parameters 
adopted for the curves are 5 = 9.8 m sec"^, if = 15 km, A6»rad = 70 K, a = 6400 km, 9o = 260 K, and A6'„ = 30 K. In the GCM 
studies, different authors define the width of the Hadley cell in different ways, so some degree of scatter is inevitable. 



zonal wind, which breaks the angular-momentum conserva- 
tion constraint in the upper-level wind, causing u to deviate 
from Eq. (52). With a different u{(j)) profile, the latitudi- 
nal dependence of temperature will change (via Eq. 53), 
and hence so will the latitudinal extent of the Hadley cell 
required to satisfy Eq. (55). Indeed, within the context of 
axisymmetric models, the addition of strong drag into the 
upper-layer flow (parameterizing turbulent mixing with the 
slower-moving surface air, for example) can lead Eq. (55) 
to predict that the Hadley cell should extend to the poles 
even for Earth's rotation rate (e.g., Farrell 1990). 

It could thus be the case that the width of the Hadley 
ceU is strongly controlled by eddies. For example, in the 
midlatitudes of Earth and Mars, baroclinic eddies generally 
accelerate the zonal flow eastward in the upper troposphere; 
in steady state, this is generally counteracted by a westward 
CorioUs acceleration, which requires an equatorward up- 
per tropospheric flow — ^backwards from the flow direction 
in the Hadley ceU. Such eddy effects can thereby terminate 
the Hadley cell, forcing its confinement to low latitudes. 
Based on this idea. Held (2000) suggested that the Hadley 
cell width is determined by the latitude beyond which the 
troposphere first becomes baroclinically unstable (requir- 
ing isentrope slopes to exceed a latitude-dependent critical 
value). Adopting the horizontal thermal gradient implied 
by the angular-momentum conserving wind (Eq. 54), mak- 
ing the small-angle approximation, and utilizing a common 
two-layer model of baroclinic instabiUty, this yields (Held 
2000) 

'^"'^[w^J ^^^^ 



in radians, where A9y is the vertical difference in potential 
temperature from the surface to the top of the Hadley cell. 
Note that the predicted dependence of (I>h on planetary ra- 
dius, gravity, rotation rate, and height of the Hadley cell is 
weaker than predicted by the Held-Hou model. Earth-based 
GCM simulations suggest that Eq. (57) may provide a better 
representation of the parameter dependences (Frierson et al. 
2007; Lu et al. 2007; Korty and Schneider 2008). Neverthe- 
less, even discrepancies with Eq. (57) are expected since the 
actual zonal wind does not follow the angular-momentum 
conserving profile (implying that the actual thermal gradi- 
ent will deviate from Eq. 54). Substantially more work is 
needed to generalize these ideas to the full range of condi- 
tions relevant for exoplanets. 

5.4. High-latitude circulations: the baroclinic zones 

For terrestrial planets heated at the equator and cooled at 
the poles, several studies suggest that the equator-to-pole 
heat engine can reside in either of two regimes depend- 
ing on rotation rate and other parameters (Del Genio and 
Suozzo 1987). When the rotation period is long, the Hadley 
cells extend nearly to the poles, dominate the equator-to- 
pole heat transport, and thereby determine the equator-to- 
pole temperature gradient. Baroclinic instabilities are sup- 
pressed because of the small latitudinal thermal gradient 
and large Rossby deformation radius (at which baroclinic 
instabilities have maximal growth rates), which exceeds the 
planetary size for slow rotators such as Titan and Venus. 
On the other hand, for rapidly rotating planets like Earth 
and Mars, the Hadley cells are confined to low latitudes; 
in the absence of eddies the mid- and high-latitudes would 
relax into a radiative-equilibrium state, leading to minimal 
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equator-to-pole heat transport and a large equator-to-pole 
temperature difference. This structure is baroclinically un- 
stable, however, and the resulting baroclinic eddies provide 
the dominant mechanism for transporting thermal energy 
from the poleward edges of the Hadley cells (~ 30° lati- 
tude for Earth) to the poles. This baroclinic heat transport 
significantly reduces the equator-to-pole temperature con- 
trast in the mid- and high-latitudes. For gravity, planetary 
radius, and heating rates relevant for Earth, the breakpoint 
between these regimes occurs at a rotation period of ~ 5-10 
days (Del Genio and Suozzo 1987). 

On rapidly rotating terrestrial planets hke Earth and 
Mars, then, the equatorial and high-latitude circulations 
fundamentally differ: the equatorial Hadley cells, while 
strongly affected by eddies, do not require eddies to ex- 
ist, nor to transport heat poleward. In contrast, the circula- 
tion and heat transport at high latitudes (poleward of ^ 30° 
latitude on Earth and Mars) fundamentally depends on the 
existence of eddies. At high latitudes, interactions of ed- 
dies with the mean flow controls the latitudinal tempera- 
ture gradient, latitudinal heat transport, and structure of the 
jet streams. This range of latitudes is called the baroclinic 
zone. 

The extent to which baroclinic eddies can reduce the 
equator-to-pole temperature gradient has important imph- 
cations for the mean climate. Everything else being equal, 
an Earth-hke planet with colder poles will develop more ex- 
tensive polar ice and be more susceptible to an ice-albedo 
feedback that triggers a globally glaciated "snowball" state 
(Spiegel et al. 2008). Likewise, a predominantly CO2 at- 
mosphere can become susceptible to atmospheric coUapse 
if the polar temperatures become sufficiently cold for CO2 
condensation. 

There is thus a desire to understand the extent to which 
baroclinic eddies can transport heat poleward. General cir- 
culation models (GCMs) attack this problem by spatially 
and temporally resolving the full life of every barocUnic 
eddy and their effect on the mean state, but this is com- 
putationally intensive and often sheds httle hght on the un- 
derlying sensitivity of the process to rotation rate and other 
parameters. 

Two simplified approaches have been advanced that illu- 
minate this issue. Although GCMs for terrestrial exoplanets 
will surely be needed in the future, simplified approaches 
can guide our understanding of such GCM results and pro- 
vide testable hypotheses regarding the dependence of heat 
transport on planetary parameters. They also may provide 
guidance for parameterizations of the latitudinal heat trans- 
port in simpler energy-balance climate models that do not 
explicitly attempt to resolve the full dynamics. We devote 
the rest of this section to discussing these simplified ap- 
proaches. 

The first simplified approach postulates that baroclinic 
eddies relax the midlatitude thermal structure into a state 
that is neutrally stable to baroclinic instabihties (Stone 
1978), a process called baroclinic adjustment. This idea 
is analogous to the concept of convective adjustment: when 



radiation or other processes drive the vertical temperature 
gradient steeper than an adiabat (dT/dz < —g/cp for 
an ideal gas, where z is height), convection ensues and 
drives the temperature profile toward the adiabat, which 
is the neutrally stable state for convection. If the convec- 
tive overtum timescales are much shorter than the radiative 
timescales, then convection overwhelms the ability of ra- 
diation to destabilize the environment, and the temperature 
profile then deviates only shghtly from an adiabat over a 
wide range of convective heat fluxes (see, e.g., §4.1). In 
a similar way, the concept of baroclinic adjustment postu- 
lates that when timescales for barochnic eddy growth are 
much shorter than the timescale for radiation to create a 
large equator-to-pole temperature contrast, the eddies will 
transport thermal energy poleward at just the rate needed to 
maintain a profile that is neutral to baroclinic instabihties. 

This idea was originally developed for a simplified two- 
layer system, for which baroclinic instability first initiates 
when the slope of isentropes exceeds ~ H/a, where H is 
a pressure scale height and a is the planetary radius. Com- 
parison of this critical isentrope slope with Earth observa- 
tions shows impressive agreement poleward of 30-40° lati- 
tude, where baroclinic instabilities are expected to be active 
(Stone 1978). A substantial hterature has subsequently de- 
veloped to explore the idea further (for a review see Zurita- 
Gotor and Lindzen 2006). The resulting latitudinal temper- 
ature gradient is then approximately H/a times the verti- 
cal potential temperature gradient. Interestingly, this the- 
ory suggests that otherwise identical planets with different 
Brunt- Vaisala frequencies (due to differing opacities, verti- 
cal heat transport by large-scale eddies, or role for latent 
heating, all of which can affect the stratification) would 
exhibit very different midlatitude temperature gradients — 
the planet with smaller stable stratification (i.e., smaller 
Brunt- Vaisala frequency; Eq. 6) exhibiting a smaller lati- 
tudinal temperature gradient. Moreover, although Earth's 
ocean transports substantial heat between the equator and 
poles, the theory also suggests that the ocean may only ex- 
ert a modest influence on the latitudinal temperature gradi- 
ents in the atmosphere (Lindzen and Farrell 1980a): in the 
absence of oceanic heat transport, the atmospheric eddies 
would simply take up the slack to maintain the atmosphere 
in the baroclinically neutral state. (There could of course 
be an indirect effect on latitudinal temperature gradients if 
removing/adding the oceans altered the tropospheric strat- 
ification.) Finally, the theory suggests that the latitudinal 
temperature gradient in the baroclinic zone does not depend 
on planetary rotation rate except indirectly via the influence 
of rotation rate on static stability. It is worth emphasizing, 
however, that even if the latitudinal temperature gradient 
in the baroclinic zone were constant with rotation rate, the 
total equator-to-pole temperature difference would still de- 
crease with decreasing rotation rate because, with decreas- 
ing rotation rate, the Hadley cell occupies a greater latitude 
range and the baroclinically adjusted region would be com- 
pressed toward the poles (see Fig. 12). 

Despite Stone (1978)'s encouraging results, several 
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complicating factors exist. First, the time-scale separa- 
tion between radiation and dynamics is much less obvi- 
ous for baroclinic adjustment than for convective adjust- 
ment. In an Earth or Mars-Uke context, the convective 
instability timescales are ^ 1 hour, the baroclinic insta- 
bility timescales are days, and the radiation timescale is 
~ 20 days (Earth) or ~ 2 days (Mars). The ability of baro- 
clinic eddies to adjust the environment to a baroclinically 
neutral state may thus be marginal, especially for planets 
with short radiative time constants (such as Mars). Sec- 
ond, the neutrally stable state in the two-layer model — 
corresponding to an isentrope slope of ~ H/ a — is an arte- 
fact of the vertical discretization in that model. When a 
multi-layer model is used, the critical isentrope slope for 
initiating barochnic instabihty decreases as the number of 
layers increases, and it approaches zero in a vertically con- 
tinuous model. One might then wonder why a baroclinic 
zone can support any latitudinal temperature contrast (as it 
obviously does on Earth and Mars). The reason is that the 
instability growth timescales are long at small isentropic 
slopes (Lindzen and Farrell 1980b); the instabilities only 
develop a substantial ability to affect the mean state when 
isentrope slopes become steep. In practice, then, the baro- 
clinic eddies are perhaps only able to relax the atmosphere 
to a state with isentrope slopes of ~ H/ a rather than some- 
thing significantly shallower. 

The second simplified approach to understanding equator- 
to-pole temperature contrasts on terrestrial planets seeks to 
describe the poleward heat transport by baroclinic eddies as 
a diffusive process (Held 1999). This approach is based on 
the idea that baroclinic eddies have maximal growth rates 
at scales close to the Rossby radius of deformation, Ln 
(see Eq. 23), and if the barochnically unstable zone has 
a latitudinal width substantially exceeding Ld then there 
will be scale separation between the mean flow and the 
eddies. The possible existence of a scale separation in 
the barochnic instability problem implies that representing 
the heat transport as a function of local mean-flow quan- 
tities (such as the mean latitudinal temperature gradient) 
is a reasonable prospect. Still, caution is warranted, since 
an inverse energy cascade (§3.6) could potentially transfer 
the energy to scales larger than Ld, thereby weakening the 
scale separation. 

The diffusive approach is typically cast in the context 
of a ID energy-balance model that seeks to determine the 
variation with latitude of the zonally averaged surface tem- 
perature. This approach has a long history for Earth climate 
studies (see review in North et al. 1981), Solar-System plan- 
ets (e.g. Hoffert et al. 1981), and even in preliminary stud- 
ies of the climates and habitability of terrestrial exoplanets 
(WiUiams and Kasting 1997; Williams and Pollard 2003; 



'As pointed out by Held (1999), this differs from many instability problems 
in fluid mechanics, such as shear instability in pipe flow or convection in a 
fluid driven by a heat flux between two plates, where maximal growth rates 
occur at length scales comparable to the domain size and no eddy/mean- 
flow scale separation occurs. 



Spiegel et al. 2008, 2009). In its simplest form, the govem- 
ing equation reads 

= V • (cDVT) + S{^) - I (58) 

where T is surface temperature, D (m^sec~^) is the dif- 
fusivity associated with heat transport by barochnic eddies, 
S{(f)) is the absorbed stellar flux as a function of latitude, 
and / is the emitted thermal flux, which is a function of tem- 
perature. In the simplest possible case, I is represented as 
a linear function of temperature, I = A + BT (North et al. 
1981), where A and B are positive constants. In Eq. (58), c 
is a heat capacity (with units Jm~^K~^) that represents the 
atmosphere and oceans (if any). The steady-state solution 
in die absence of transport (£> = 0) is T = - A]/B. 
On the other hand, when transport dominates {D — > oo), 
the solution yields constant T. 

The term 5(0) includes the effect of latitudinally varying 
albedo (e.g., due to ice cover), but if albedo is constant with 
latitude, then S represents the latitudinally varying insola- 
tion. To schematically illustrate the effect of heat transport 
on the equator-to-pole temperature gradient, parameterize S 
as a constant plus a term proportional to the Legendre poly- 
nomial P2(cos (f)). In this case, and if D and c are constant, 
the equation has a steady analytic solution with an equator- 
to-pole temperature difference given by (Held 1999) 

where a is the planetary radius and ATj-ad is the equator-to- 
pole difference in radiative-equihbrium temperature. The 
equation implies that atmospheric heat transport signifi- 
cantly influences the equator-to-pole temperature difference 
when the diffiisivity exceeds ~ Ba? / (6c). For Earth, B « 
2Wm-2K-i, a = 6400km, and c w lO^Jm-^K^i, 
suggesting that atmospheric transport becomes important 
whenD > lO'^m^sec-^ 

The question comes down to how to determine the diffu- 
sivity, D. In Earth models, D is typically chosen by tuning 
the models to match the current climate, and then that value 
of D is used to explore the regimes of other climates (e.g. 
North 1975). However, that approach sidesteps the underly- 
ing physics and prevents an extrapolation to other planetary 
environments. 

Motivated by an interest in understanding the feedback 
between barochnic instabihties and the mean state, substan- 
tial effort has been devoted to determining the dependence 
of the diffusivity on control parameters such as rotation 
rate and latitudinal thermal gradient (for a review see Held 
1999). This work can help guide efforts to understand the 
temperature distribution on rapidly rotating exoplanets. We 
expect that the diffusivity will scale as 

D W Weddy^eddy (60) 

where Weddy and Z/eddy are the characteristic velocities 
and horizontal sizes of the heat-transporting eddies. Be- 
cause baroclinic eddies will be in near-geostrophic balance 
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(Eq. 21) on a rapidly rotating planet, the characteristic eddy 
velocities, Uoddy , will relate to characteristic eddy potential- 
temperature perturbations via the thermal-wind rela- 
tion, giving Uoddy ^ .g^cddy'*cddy/(/icddy), whcrC ft-oddy 

is the characteristic vertical thickness of the eddies. Under 
the assumption that the ratio of vertical to horizontal scales 
is /leddy/ioddy ~ f/N (Chamcy 1971; Haynes 2005), this 
yields Weddy ~ 5^cddy/(^o-^) (where N is Brunt- Vaisala 
frequency), which simply states that eddies with larger ther- 
mal perturbations will also have larger velocity perturba- 
tions. Under the assumption that the thermal perturbations 
scale as ^eddy ~ LeddydO/dy, we then have 
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Greater thermal gradients and eddy-length scales lead to 
greater diffusivity and, importantly, the diffiisivity depends 
on the square of the eddy size. 

Several proposals for the relevant eddy size have been 
put forward, which we summarize in Table 2. Based on 
the idea that baroclinic instabilities have the greatest growth 
rates for lengths comparable to the deformation radius. 
Stone (1972) suggested that .Leddy is the deformation ra- 
dius, NH/f. On the other hand, baroclinic eddies could 
energize an inverse cascade, potentially causing the domi- 
nant heat- transporting eddies to have sizes exceeding Ljj. 
In the limit of this process (Green 1970), the eddies would 
reach the width of the barocUnic zone, .Lzone (potentially 
close to a planetary radius for a planet with a narrow Hadley 
cell). This simply leads to Eq. (61) with Lcddy = -^zonc- In 
contrast. Held and Larichev (1996) argued that the inverse 
cascade would produce an eddy scale not of order Lzone but 
instead of order the Rhines scale, (?icddy 

Finally, 

Barry et al. (2002) used heat-engine arguments to propose 
that 
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where gnet is the net radiative heating/cooling per mass that 
the eddy fluxes are balancing, a is the planetary radius, and 
e is a constant of order unity. Equation 62 is probably most 
robust for the dependences on heating rate q^^t and thermal 
gradient 89 /dy, which they varied by factors of ~ 200 and 
6, respectively; planetary radius and rotation rates were var- 
ied by only ~ 70%, so the dependences on those parameters 
should be considered tentative. 

As can be seen in Table 2, these proposals have diver- 
gent implications for the dependence of diffusivity on back- 
ground parameters. Stone (1972)'s diffusivity is propor- 
tional to the latitudinal temperature gradient and, because 
of the variation of Ld with rotation rate, inversely propor- 
tional to the square of the planetary rotation rate. Green 
(1970)'s diffusivity likewise scales with the latitudinal tem- 
perature gradient; it contains no explicit dependence on the 
planetary rotation rate, but a rotation-rate dependence could 
enter because the Hadley cell shrinks and L^onc increases 
with increasing rotation rate. Held and Larichev (1996)'s 



diffusivity has the same rotation-rate dependence as that of 
Stone (1972), but it has a much stronger depenence on grav- 
ity and temperature gradient [g^ and {dO/dy)^\. Moreover, 
it increases with decreasing static stability as N~^, un- 
like Stone's diffusivity that scales with N Finally, Barry 
et al. (2002)'s diffusivity suggests a somewhat weaker de- 
pendence, scaUng with as {dO/dy)^/^ and fi"^/^. Note 
that the rotation-rate dependences described above should 
be treated with caution, because altering the rotation rate or 
gravity could change the circulation in a way that alters the 
Brunt- Vaisala frequency, leading to additional changes in 
the diffusivity. 

Additional work is needed to determine which of these 
schemes (if any) is valid over a wide range of parameters 
relevant to terrestrial exoplanets. Most of the work per- 
formed to test the schemes of Green (1970), Stone (1972), 
Held and Larichev (1996) and others has adopted simpli- 
fied two-layer quasi-geostrophic models^^ under idealized 
assumptions such as planar geometry with constant Corio- 
lis parameter (a so-called "f-plane"). Because variation of 
/ with latitude alters the properties of barocUnic instabiUty, 
schemes developed for constant / may not translate directly 
into a planetary context. More recently, two-layer planar 
models with non-zero /? have been explored by Thompson 
and Young (2007) and Zurita-Gotor and Vallis (2009), while 
full three-dimensional GCM calculations on a sphere (with 
~ 30 layers) investigating the latitudinal heat flux have been 
performed by Barry et al. (2002) and Schneider and Walker 
(2008) under conditions relevant to Earth. Additional simu- 
lations in the spirit of Barry et al. (2002) and Schneider and 
Walker (2008), considering a wider range of planetary and 
atmospheric parameters, can clarify the true sensitivities of 
the heat transport rates in three-dimensional atmospheres. 

5.5. Slowly rotating regime 

At slow rotation rates, the GCM simulations shown in 
Figs. 11-12 develop near-global Hadley cells with high- 
latitude jets, but the wind remains weak at the equator. In 
contrast. Titan and Venus (with rotation periods of 16 and 
243 days, respectively) have robust (~100msec~^) super- 
rotating winds in their equatorial upper tropospheres. Be- 
cause the equator is the region of the planet lying farthest 
from the rotation axis, such a flow contains a local maxi- 
mum of angular momentum at the equator. Axisymmetric 
Hadley circulations (§5.3) cannot produce such superrota- 
tion; rather, up-gradient transport of momentum by eddies 



^*The Af~^ dependence in Held and Larichev (1996)'s gives the impres- 
sion that the diffusivity becomes unbounded as the atmospheric verti- 
cal temperature profile becomes neutrally stable (i.e. as JV — > 0), but 
this is misleading. In reality, one expects the latitudinal temperature 
gradient and potential energy available for driving baroclinic instabili- 
ties to decrease with decreasing N. Noting that the slope of isentropes 
is me = {d6/dy){d6/dz)~^, one can re-express Held and Larichev 
(1996)'s diffusivity as scaling with mgN^. Thus, at constant isentrope 
slope, the diffusivity properly drops to zero as the vertical temperature 
profile becomes neutrally stable. 

'•'That is, models where geostrophic balance is imposed as an external con- 
straint. 
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Table 2 

Proposed diffusivities for high-latitude heat transport on terrestial planets 



Scheme ieddy D 



Stone (1972) Lu 

Green (1970) L.one izonej^ff 



Held and Larichev (1996) (ff ) 

Barry etal. (2002) L, f) (|)'^' 



is required. This process remains poorly understood. One 
class of models suggests that this transport occurs from high 
latitudes; for example, the high- latitude jets that result from 
the Hadley circulation (which can be seen in Fig. 11) expe- 
rience a large-scale shear instability that pumps eddy mo- 
mentum toward the equator, generating the equatorial su- 
perrotation (Del Genio and Zhou 1996). Alternatively, ther- 
mal tide or wave interaction could induce momentum trans- 
ports that generate the equatorial superrotation (e.g., Fels 
and Lindzen 1974). Recently, a variety of simplified Venus 
and Titan GCMs have been developed that show encourag- 
ing progress in capturing the superrotation (e.g., Yamamoto 
and Takahashi 2006; Richardson et al. 2007; Lee et al. 2007; 
Herrnstein and Dowling 2007). 

The Venus/Titan superrotation problem may have impor- 
tant impUcations for understanding the circulation of syn- 
chronously locked exoplanets. All published 3D circulation 
models of synchronously locked hot Jupiters (§4.3), and 
even the few published studies of sychrounously locked ter- 
restrial planets (Joshi et al. 1997; Joshi 2003, see §5.6), de- 
velop robust equatorial superrotation. The day-night heat- 
ing pattern associated with synchronous locking should 
generate thermal tides and, in analogy with Venus and Ti- 
tan, these could be relevant in driving the superrotation 
in the 3D exoplanet models. SimpUfied Earth-based two- 
layer calculations that include longitudinally varying heat- 
ing, and which also develop equatorial superrotation, seem 
to support this possibility (Suarez and Duffy 1992; Sara- 
vanan 1993). 

5.6. Unusual forcing regimes 

Our discussion of circulation regimes on terrestrial plan- 
ets has so far largely focused on annual-mean forcing con- 
ditions. On Earth, seasonal variations represent relatively 
modest perturbations around the annual-mean climate due 
to Earth's 23.5° obhquity. On longer timescales, the secu- 
lar evolution of Earth's orbital elements is thought to be re- 
sponsible for paleoclimatic trends, such as ice ages, accord- 
ing to Milankovitch's interpretation (e.g. Kawamura et al. 
2007). While it may be surprising to refer to Earth's seasons 



or ice ages as minor events from a human perspective, it is 
clear that they constitute rather mild versions of the more 
diverse astronomical forcing regimes expected to occur on 
extrasolar worlds. 

In principle, terrestrial exoplanets could possess large 
obliquities {i 90deg) and eccentricities (e 1). This 
would result in forcing conditions with substantial seasonal 
variations around the annual mean, such as order-of magni- 
tude variations in the global insolation over the orbital pe- 
riod at large eccentricities. Even the annual mean climate 
can be dramatically affected under unusual forcing condi- 
tions, for instance at large obliquities when the poles receive 
more annual-mean insolation than the planetary equator (for 
i > 54°). Yet another unusual forcing regime occurs if the 
terrestrial planet is tidally locked to its parent star and thus 
possesses permanent day and night sides much Uke the hot 
Jupiters discussed in §4.3. Relatively few studies of circu- 
lation regimes under such unusual forcing conditions have 
been carried out to date. 

WilHams and Kasting (1997) and Spiegel et al. (2009) 
investigated the climate of oblique Earth-like planets with 
simple, diffiisive energy-balance models of the type de- 
scribed by Eq. (58). While seasonal variations were found 
to be severe at high obliquity, Wilhams and Kasting (1997) 
concluded that highly obUque planets could nevertheless 
remain regionally habitable, especially if they possessed 
thick C02-enriched atmospheres with large thermal iner- 
tia relative to Earth, as may indeed be expected for terres- 
trial planets in the outermost regions of a system's habit- 
able zone. Spiegel et al. (2009) confirmed these results and 
highlighted the risks that global glaciation events and partial 
atmospheric CO2 collapse constitute for highly oblique ter- 
restrial exoplanets. Williams and Pollard (2003) presented 
a much more detailed set of three-dimensional cUmate and 
circulation models for Earth-like planets at various obliq- 
uities. Despite large seasonal variations at extreme obliq- 
uities, WilUams and Pollard (2003) found no evidence for 
any runaway greenhouse effect or global glaciation event 
in their climate models and thus concluded that Earth-like 
planets should generally be hospitable to life even at high 
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obliquity. 

Williams and Pollard (2002) also studied the climate 
on eccentric versions of the Earth with a combination 
of detailed three-dimensional climate models and simpler 
energy-balance models. They showed that the strong vari- 
ations in insolation occurring at large eccentricities can be 
rather efficiently buffered by the large thermal inertia of 
the atmosphere+ocean climate system. These authors thus 
argued that the annually averaged insolation, which is an 
increasing function of e, may be the most meaningful quan- 
tity to describe in simple terms the climate of eccentric, 
thermally blanketed Earth-Uke planets. It is presumably 
the case that terrestrial planets with reduced thermal iner- 
tia from lower atmospheric/oceanic masses would be more 
strongly affected by such variations in insolation. 

It should be emphasized that, despite the existing work 
on oblique or eccentric versions of the Earth, a fundamen- 
tal understanding of circulation regimes under such unusual 
forcing conditions is still largely missing. For instance, 
variations in the Hadley circulation regime or the baro- 
clinic transport efficiency, expected as forcing conditions 
vary along the orbit, have not been thoroughly explored. 
Similarly, the combined effects of a substantial obliquity 
and eccentricity have been ignored. As the prospects for 
finding and characterizing exotic versions of the Earth im- 
prove, these issues should become the subject of increasing 
scrutiny. 

One of the best observational prospects for the next 
decade is the discovery of tidally locked terrestrial plan- 
ets around nearby M-dwarfs (e.g., Irwin et al. 2009; Sea- 
ger et al. 2008). In anticipation of such discoveries, Joshi 
et al. (1997) presented a careful investigation of the circu- 
lation regime on this class of unusually forced planets, with 
permanent day and night sides. The specific focus of their 
work was to determine the conditions under which CO2 at- 
mospheric collapse could occur on the cold night sides of 
such planets and act as a trap for this important greenhouse 
gas, even when circulation and heat transport are present. 
Using a simplified general circulation model, these authors 
explored the problem's parameter space for various global 
planetary and atmospheric attributes and included a study 
of the potential risks posed by violent flares from the stel- 
lar host. The unusual circulation regime that emerged from 
this work was composed of a direct circulation cell at sur- 
face levels (i.e. equatorial day to night transport with po- 
lar return) and a superrotating wind higher up in the atmo- 
sphere. While detrimental atmospheric collapse did occur 
under some combinations of planetary attributes (e.g., for 
thin atmospheres), the authors concluded that efficient at- 
mospheric heat transport was sufficient to prevent collapse 
under a variety of plausible scenarios. More recently, Joshi 
(2003) revisited some of these results with a much more 
detailed climate model with explicit treatments of non-gray 
radiative transfer and a hydrological cycle. It is likely that 
this interesting circulation regime will be reconsidered in 
the next few years as the discovery of such exotic worlds is 
on our horizon. 



6. RECENT HIGHLIGHTS 

Since the first discovery of giant exoplanets around Sun- 
like stars via the Doppler velocity technique (Mayor and 
Queloz 1995; Marcy and Butler 1996), exoplanet research 
has experienced a spectacular series of observational break- 
throughs. The first discovery of a transiting hot Jupiter 
(Charbonneau et al. 2000; Henry et al. 2000) opened the 
door to a wide range of clever techniques for character- 
izing such transiting planets. The drop in flux that oc- 
curs during secondary eclipse, when the planet passes be- 
hind its star, led to direct infrared detections of the thermal 
flux from the planet's dayside (Deming et al. 2005; Char- 
bormeau et al. 2005). This was quickly followed by the 
detection of day/night temperature variations (Harrington 
et al. 2006), detailed phase curve observations (e.g. Knut- 
son et al. 2007, 2009b; Cowan et al. 2007), infrared spectral 
and photometric measurements (Grillmair et al. 2007, 2008; 
Richardson et al. 2007; Charbonneau et al. 2008; Knutson 
et al. 2008) and a variety of transit spectroscopic constraints 
(Tinetti et al. 2007; Swain et al. 2008; Barman 2008). Col- 
lectively, these observations help to constrain the compo- 
sition, albedo, three-dimensional (3D) temperature struc- 
ture, and hence the atmospheric circulation regime of hot 
Jupiters. 

The reality of atmospheric circulation on these objects is 
suggested by exquisite light curves of HD 189733b showing 
modest day-night infrared brightness temperature variations 
and displacement of the longitudes of minimum/maximum 
flux from the antistellar/substellar points, presumably the 
result of an efficient circulation able to distort the tem- 
perature pattern (Knutson et al. 2007, 2009b). Dayside 
photometry also hints at the existence of an atmospheric 
circulation on this planet (Barman 2008). Nevertheless, 
other hot Jupiters, such as Ups And b and HD 179949b, 
apparently exhibit large day-night temperature variations 
with no discernable displacement of the hot regions from 
the substellar point (Harrington et al. 2006; Cowan et al. 
2007). Secondary-eclipse photometry suggests that some 
hot Jupiters have dayside temperatures that decrease with 
altitude, while other hot Jupiters appear to exhibit dayside 
thermal inversion layers where temperatures increase with 
altitude. Some authors have suggested that these two issues 
are hnked (Fortney et al. 2008), although subtantial addi- 
tional observations are needed for a robust assessment. Un- 
derstanding the dependence of infrared phase variations and 
dayside temperature profiles on planetary orbital and phys- 
ical properties remains an ongoing observational challenge 
for the coming decade. 

In parallel, this observational vanguard has triggered a 
growing body of theoretical and modeling studies to in- 
vestigate plausible atmospheric circulation patterns on hot 
Jupiters. The models agree on some issues, such as that syn- 
chronously rotating hot Jupiters on ~ 3-day orbits should 
exhibit fast winds with only a few broad zonal jets. On the 
other hand, the models disagree on other issues, such as the 
details of the 3D flow patterns and the extent to which sig- 
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nificant global-scale time variability is likely. The effect of 
stellar flux, planetary rotation rate, obliquity, atmospheric 
composition, and orbital eccentricity on the circulation have 
barely been explored. Further theoretical work and compar- 
ison with observations should help move our understanding 
onto a firmer foundation. 

7. FUTURE PROSPECTS 

Over the next decade, the observational characteriza- 
tion of hot Jupiters and Neptunes will continue apace, 
aided especially by the warm Spitzer mission and subse- 
quently JWST. Equally exciting, NASA's Kepler mission 
and groundbased surveys (e.g., MEarth; see Irwin et al. 
2009) may soon lead to the discovery of not only additional 
super-Earths but Earth-sized terrestrial exoplanets, includ- 
ing "hot Earths" as well as terrestrial planets cool enough 
to lie within the habitable zones of their stars. If sufficiently 
favorable systems are discovered, follow-up by JWST may 
allow basic characterization of their spectra and light curves 
(albeit at significant resources per planet), providing con- 
straints on their atmospheric properties including their com- 
position, chmate, circulation, and the extent to which they 
may be habitable. Finally, over the next decade, signif- 
icant opportunities exist for expanding our understanding 
of basic dynamical processes by widening the parameters 
adopted in circulation models beyond those typically ex- 
plored in the study of Solar-System planets. 

To move the observational characterization of exoplanet 
atmospheric circulation into the next generation, we rec- 
ommend full-orbit light curves for a variety of hot Jupiters 
and Neptunes (and eventually terrestrial planets). Target ob- 
jects should span a range of incident stellar fluxes, planetary 
masses, and orbital periods and eccentricities (among other 
parameters) to characterize planetary diversity. To constrain 
the planetary thermal-energy budgets, Ught curves should 
sample the blackbody peaks at several wavelengths, includ- 
ing both absorption bands and low-opacity windows as ex- 
pected from theoretical models; this will provide constraints 
on the day-night temperature contrasts — and hence day- 
night heat transport — at a range of atmospheric pressure 
levels. Dayside infrared spectra obtained from secondary- 
eclipse measurements, while not necessarily diagnostic of 
the atmospheric circulation in isolation, will provide crucial 
constraints when combined with infrared light curves. By 
providing constraints on atmospheric composition, transit 
spectroscopy will enable much better interpretation of light 
curves and secondary-ecUpse spectra. Secondary ecUpses 
(or fall light curves) that are observed repeatedly will pro- 
vide powerful constraints on global-scale atmospheric vari- 
ability. Kepler may prove particularly useful in this regard, 
with its ability to observe hundreds of secondary eclipses 
of hot Jupiters (such as HAT-P-7b; Borucki and colleagues 
2009) in its field of view. Such an extensive temporal 
dataset could provide information on the frequency spec- 
trum of atmospheric variability (if any), which would help 
to identify the mechanisms of variability as well as the 



background state of the mean circulation. Other novel tech- 
niques, such as the possibility of separately measuring tran- 
sit spectra on the leading and trailing limbs and the possibil- 
ity of detecting the Doppler shifts associated with planetary 
rotation and/or winds, should be considered. 

Specific research questions for the next decade include 
but are not Umited to the following: 

• How warm are the nightsides of hot Jupiters? How 
does this depend on the stellar flux, atmospheric com- 
position, and other parameters? 

• What are the dynamical mechanisms for shifting the 
hot and cold points away from the substellar and an- 
tistellar points, respectively? Why does HD 189733b 
appear to have both its hottest and coldest points on 
the same hemisphere? 

• What is the observational and physical relationship 
(if any) between the amplitude of the day-night tem- 
perature contrast and existence or absence of a day- 
side temperature inversion that has been inferred on 
some hot Jupiters? 

• How does the expected atmospheric circulation 
regime depend on planetary size, rotation rate, grav- 
ity, obliquity, atmospheric composition, and incident 
stellar flux? 

• What are the important dissipation mechanisms in 
the atmospheres of hot Jupiters? What are the rela- 
tive roles of turbulence, shocks, and magnetohydro- 
dynamic processes in frictionally braking the flows? 

• How coupled are the atmospheric flows on hot 
Jupiters to convection in the planetary interior? How 
are momentum, heat, and constituents transported 
across the deep radiative zone in hot Jupiters? 

• Does the atmospheric circulation influence the long- 
term evolution of hot Jupiters? Can downward trans- 
port of atmospheric energy into the interior help 
explain the radii of some anomalously large hot 
Jupiters? 

• What is the mechanism for equatorial superrotation 
that occurs in most 3D models of synchronously ro- 
tating hot Jupiters? Does such superrotation exist on 
real hot Jupiters? 

• Are the atmospheres of hot Jupiters strongly time 
variable on the global scale, and if so, what is the 
mechanism for the variabiUty? 

• How does the circulation regime of planets in highly 
eccentric orbits differ from that of planets in circular 

orbits? 

• To what extent can obliquity, rotation rate, atmo- 
spheric mass, and atmospheric circulation proper- 
ties be inferred from disk-integrated spectra or Ught 
curves? 
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• To what extent, if any, does the atmospheric circula- 
tion influence the escape of an atmosphere to space? 

• Do the deep molecular envelopes of giant planets dif- 
ferentiaUy rotate? 

• What controls the wind speeds on giant planets? Why 
are Neptune's winds faster than Jupiter's, and what 
can this teach us about strongly irradiated planets? 

• What can observations of chemical disequilibrium 
species tell us about the atmospheric circulation of 
giant planets and brown dwarfs? 

• What is the role of clouds in affecting the circulation, 
climate, evolution, and observable properties of exo- 
planets? 

• For terrestrial planets, how does the atmospheric cir- 
culation affect the boundaries of the classical habit- 
able zone? 

• What is the role of the atmospheric circulation in 
affecting the ice-albedo feedback, atmospheric col- 
lapse, runaway greenhouse, and other chmate feed- 
backs? 

• Can adequate scaling theories be developed to predict 
day-night or equator-to-pole temperature contrasts as 
a function of planetary and atmospheric parameters? 

• What is the dynamical role of a surface (i.e. the 
ground) in the atmospheric circulation? For super- 
Earths, is there a critical atmospheric mass beyond 
which the surface becomes unimportant and the cir- 
culation behaves Uke that of a giant planet? 

Ultimately, unraveling the atmospheric circulation and 
climate of exoplanets from global-scale observations will 
be a difficult yet exciting challenge. Degeneracies of in- 
terpretation will undoubtedly exist (e.g., multiple circula- 
tion patterns explaining a given light curve), and forward 
progress will require not only high-quality observations but 
a careful exploration of a hierarchy of models so that the 
nature of these degeneracies can be understood. Despite 
the challenge, the potential payoff will be the unleashing 
of planetary meteorology beyond the confines of our So- 
lar System, leading to not only an improved understand- 
ing of basic circulation mechanisms (and how they may 
vary with planetary rotation rate, gravity, incident stellar 
flux, and other parameters) but a glimpse of the actual at- 
mospheric circulations, climate, and habitability of planets 
orbiting other stars in our neighborhood of the Milky Way. 
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